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16190 THERE WERE GIANTS ON THE EARTH IN THOSE DAY S re aorta 
KEY "WORDS: _ Adobe; Archaeology; Construction methods; Drainage; we 
Geotechnical engineering; Historic sites; History; — Masonry; © Mounds pe = 
H (archaeolog ); Piles; Roads; Sociological factors 
| ABSTRACT: There are numerous examples of noteworthy pre- pre-Columbian 
Psi geotechnology in the Americas. These include earth mounds, stone temple pyramids, — 
| pile foundations, roads, retaining walls, terraces, irrigation structures, and small dams. 
' | They exhibit well-conceived designs and organized construction. Some investigators 
‘speculate that -hese sophisticated works reflect unrecorded technology from the Old 
_ World or possibly extra-terrestial intervention. However, the responsibility can be 
_ found in both indigeneous, ingeneous responses to local social and economic problems 
- as well as well-organized community efforts. The lack of a formal system of logic, 
— computation, and written communication inhibited the development and 
interaction of the embryonic technologies. Their decay was completed by the European — 
REFERENCE: , “There Were Giants the Earth i in Those Days,’ 
Journal of the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT4, 
Paper 16190, April, 1981, pp. 383- 419 
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16178 WETTING AND DRYING ON SHEAR STRENGTH 
KEY WORDS: Climatic changes; 
Strength; Soil mechanics; Strength; -Stress- curves; T Time. 


WwW 
_ ABSTRACT: Climatic Pony such as repeated wetting and drying over geological | 
ages, lead to the degeneration of the parent rock to form soil. Simultaneously these 
actions can also cause some aggregation of soil particles and the production of bonds, 
_ called desiccation bonds, which impart an intrinsic effective stress to the soil. This 
_ intrinsic effective stress can influence the shear strength behavior of soils subjected to 
_ such climatic actions. When a soil is subjected to repeated wetting and drying in the 
| laboratory, an intrinsic effective stress is imparted to it which results in greater shear 
_ strength and in a stiffer stress-strain. -Tesponse. This intrinsic effective stress may be 
| REFERENCE: Allam, Mehter Mohamed, and Sridharan, Asuri, “Effect of Wetting | 
oh and Drying on Shear Strength,” Journal of the Geotechnical Engineering Division, 
ASCE, Vol. 107, No. GT4, Proc. 16178, April, 1981, pp. 42 421- 1-438 
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16173 IN NATURAL SAND DEPOSITS 210 u 
KEY ‘WORDS: Dry density; Danes; Hardening; I Impregnation; Laboratory © 
| tests; Peel tests; Porosity; Relative density; Relative density determination; 
a . ABSTRACT: The proposed technique enables porosity measurement to be made on 
= p undisturbed samples taken from sand deposits. It is based on the impregnation i 
_ hardening of sand in the field and the preparation of small cube samples which are 
- tested in the laboratory. The method proved successful in both dry and wet sand. Its 
_ application is demonstrated in the study of porosity variations in a dune. The internal 7 
_ structure of the dune was studied by means of replicas or “peels” which reveal in 
_ detail the layering of the dune sand. A definite relationship was established between __ 
_ porosity and relative density of the sand and its internal structure and mode of _ 
_ deposition. Very high densities were found in sand deposited by accretion on the © 
dune’s back. Relative densities in accretion deposits approach and even exceed 100 


REFERENCE: Denekamp, S. Tsur-Lavie, Y., “Measurement of Porosity in 
Natural Sand Desposits,” Journal of the Geotechnical Engineering Division, ASC E, 
107, No. GT4, Proc. Paper 16173, April, 1981, pp. 439- IneA , 


= 
= 
* 
| 
| 


16179 STOCHASTIC. FEM IN SETTLEMENT PREDICTIONS RT 
= | KEY WORDS: Deformation; Finite element method; Footings; Foundation __ 
design; Foundations; Probability theory; Settlement analysis; Settlement 
measurement; Settlement (structural); Stochastic processes 
ABSTRACT: Stochastic finite element analysis is used to predict uncertainties in total ra 
and differential settlement under a large flexible footing. The results are compared with 
one-dimensional stochastic solutions already in the literature. Differences between the _ 
one- and two-dimensional analyses, particularly for differential settlement, are distinct. 
These differences seem primarily attributable to randomness in the stress field which 
cannot be included in one-dimensional models, and to mechanistic correlations by 
common dependence on the realizations of particular random variables. In principle, — 
second-moment techniques can be extended to a broad range of analyses now 
performed using fi finite element and finite difference techni ues. 
REFERENCE: Beacher, Gregory B.,. and Ingra, ‘theme S., “Stochastic FEM In ay 
- Settlement Predictions,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 
107, No. GT4, Proc. Paper 16179, April, 1981, pp. 449-463 


16174 SOIL ORGANIC-CONTENT MEASUREMENT vb 
_ KEY WORDS: Decomposing organic matter; Decomposition; Ignition tests; 
Measurement; Organic particle composites; Organic soils; Soil 
Soil mechanics; Soil microorganisms; Soil physical 
| ABSTRACT: | Problems with organic soils are often avoided by engineers because of 

uncertainties in organic-content measurement and lack of information on potential : 


__ decomposition. High temperatures, used in the ignition test for measurement of the 
organic fraction, require that a correction be made for loss of surface hydration water 
_ from the mineral solids. Model soils, prepared from kaolinite and pulp fiber, and the _ 
_ dehydration curves for these materials have permitted development of a correction 


Suggest that this correction factor may be suitable for natural organic soils. The 
| ignition test does not distinguish between undecomposed organic matter and 
‘B microorganisms formed during decomposition; hence, the degree of decomposition has 
not been available. A_ relationship—between microorganism cell content, 


- f _ factor for use with the ignition test. A review of published data and an error analysis “a 


b undecomposed organic matter, and the mineral fraction—has been developed which 

permits computation of the degree of 
‘REFERENCE: Al- Kafaji, Abdul- Amir Ww. N., and Andersland, Orlando B., “Ignition 

a ‘Test for Soil Organic-Content Measurement,” Journal of the Geotechnical Engineering : 
{ Division, ASCE, Vol. 107, No. GT4, Proc. Paper 16174, April, 1981, pp. 465- 479 £ 


16180 MODELING OF COAL WASTE EMBANKMENTS 
KEY W ORDS: Centrifugation; Coal; Embankments; Laboratory 
varying geometries 
were contuesal of coal waste material, accelerated to 100-120, and subjected to 
various reservoir depths until failure occurred. Complimentary laboratory tests were | 
performed on the waste materials to obtain strength and permeability parameters for 
analytical calculations to predict observed centrifugal model behavior. rane ; 
occurred in model embankments with 1 Vertical:1.5 Horizontal and 1 Vertical:2 — 
Horizontal slopes when the phreatic surface due to seepage exited on the downstream — 
slope and initiated sloughing at the toe. Positive seepage control, i.e., a toe drain, — 
_ successfully prevented sloughing and failures. Slope stability analyses agreed well, 
albeit conservatively, with observed model performance. These models verified existing — 
- criteria of minimum rong factor of 1.3 and the necessity of positive seepage control. >i 
REFERENCE: AL Hussaiai, Mosaid M., Goodings, Deborah J., Schofield, Andrew N., 
and Townsend, Frank C., “Centrifuge Modeling of Coal Waste Embankments,” ‘rt 
the Geotechnical Engineering Division, ASCE, Vol. 107, No. 
16160, April, 1981, pp. 481-499 
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16212 EARTHQUAKE-IN 


KEY WORDS: Case reports; Cyclic leads; Earthquakes; Guatemala; 
Liquefaction; Penetration resistance; Pumice; Sand; Sandy soils; Weight 7 
mass 

_ ABSTRACT: One of the effects of the February 4, 1976, Guatemala pene was 1 
the extensive liquefaction which occurred at the settlement of La Playa on the © 
northeast shore of Lake Amatitlan. Ground accelerations in the area resulting from the ag 
magnitude 7.5 earthquake are estimated to be of the order of 0.12 g to 0.15 8. The 
characteristics of the sand are investigated, including its standard penetration resistance 

and its cyclic loading characteristics. Despite the fact that the sand is somewhat lighter S: 
in weight than sand deposits which have liquefied in other earthquakes, its liquefaction — 
characteristics are apparently influenced by the same factors as other sand deposits; its a 
overall behavior is consistent with that exhibited by other sands. The high degree of ‘ 
liquefaction at the La Playa site was probably due in large measure to the lightweight ie. 
nature of the pumi eous sands. 

_ REFERENCE: Seed, H. Bolton, Arango, ade Chan, Clarence K., Gomez-Masso, 7 
Alberto, and Ascoli, Rebecca Grant, “Earthquake-Induced Liquefaction Near Lake 
_Amatitlan, Guatemala,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 
(107, No. GT4, Proc. Paper 16212, 2, April, 1981, PP. 501- 
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In In accordance with the October, 1970 action of the ASCE Board of Direction, which stated 
- that all publications of the Society should list all measurements in both U.S. Customary and 
SI (International System) units, the following list contains conversion factors to enable readers 
_ to compute the SI unit values of measurements. A complete guide to the SI system and its | 
use has been published by the American Society for Testing and Materials. Copies of this 
_ publication (ASTM E-380) can be purchased from ASCE at a price of $3.00 each; orders must 
All authors of Journal papers are being asked to prepare their papers in this dual-unit format. 
— provide preliminary assistance to authors, the following list of co list of conversion factors and guides 


miles (miles) kilometers (km) 
inches (sq in. square millimeters (mm 
square feet (sq ft) = Square meters (m 
square yards (sq yd) square meters (m 
square miles (sq miles) square kilometers (km 
acres (acre) hectares (ha) 
inches (cu in. cubic millimeters (m: (mm* 


cubic yards (cu yd) cubic meters (m 
pounds (Ib) mass kilograms (kg) 
tons (ton) mass kilograms (kg) 
pound force (bf) newtons (N) 
kilogram force (kgf) > newtons (N) 
pounds per square foot (p (psf)  pascals (Pa) 
pounds per square inch (psi) wugee kilopascals (kPa 
U.S. gallons (gal) mires liters(L) 
_acre-fe feet (acre-ft) nes 


| 
0.093 
0.836 
0.405 — 
16400 
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his father, George B. Sowers, F. ASCE, on inspections of tunnel, foundation, — 
- and port construction. His undergraduate education was at Case Institute of ; 
Technology; his graduate work was at Harvard where he was ome student: of 
Karl Terzaghi and Arthur Casagrande ; rise 
George worked as a hydraulic and hydrologic engineer with the United States 
TVA and served in the United States Navy during WW Il, 
; _ Since 1947 he has pursued a dual career: teacher and consultant, presently — 
Regents Professor of Civil Engineering at the Georgia Institute of Technology, 
he pioneered geotechnical engineering education in the southeast. His — 
__ Introductory Soil Mechanics and Foundations is in its fourth United States 
~ edition, plus Spanish, international, and even pirated Taiwan editions. In addition, 
he has published more than 100 technical publications, eee 
George is a cofounder of Law Engineering, one of the largest geotechnical 
environmental and materials engineering firms of the United States. He has 
a consultant on geotechnical design throughout the United ‘States and in 


George Sowers has always a civil ¢ engineer. As a child he 


=. 


__He has been President of the Georgia Section ASCE and Chairman of the 

1 Geotechnical Engineering Division, ASCE. Presently he is the Vice President 
for North Americaofthe ISSMFE. 

a He is married; his wife, Frances, has practiced —_— ay 


_ have four children, two of whom are geologists. | 
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THERE GIANTS ON THE EARTH” 


in the Americas. Some reflect well-organized technologies comparable to a? 
4 of the Old World of the same time period. Earth fills, stone and adobe structures, 


highways, and water management systems exhibit both well- cp designs — 


- q technologies. They were vulnerable to decay and were vi 
Flying high above the earth, about 370 km south- southeast of inna: Peru, 
Baar crossing a broad, arid valley 50 km from the sea, one can see patterns 
in the | valley floor below. Even the earth resources satellites, 950 km (570 miles) 
high, can discern straight lines, some several kilometers or miles long. At lower _ 
_ altitudes, the lines become more distinct, up to 50 m wide, in strange criss-cross 
: _ patterns on the 60 km (35 mile) long, 2 km ae a mile) wide braided- co 
terrace 500m (1,600 ft) above sealevelh 


"Regents Prof. of Civ. Engrg., Georgia Inst. of Tech., Atlanta, Ga. 30332; also Sr. 
Vice Pres., Law Engrg. Testing Co., 2749 Delk Road, S.E., Marietta, Ga. 30067, 0 
Note.—Discussion open until September 1, 1981. To extend the closing date one month, 
. a written request must be filed with the Manager of Technical and Professional Publications, = 
7 4 ASCE. Manuscript was submitted for review for possible publication on February 5, 
4 1980. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 


of the American Society of Civil 5 alae Vol. 107, No. GT4, April, 1981. 
-6405 /81 /0004-0385 / $01. 00. 
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7 -_ a ew investigators speculate that such sophisticated works represent unre- j 

_ corded technology transfer from the Old World; some pseudohistorians — | 

visits from outer space. However, the evidence demonstrates that indigenous, ‘ 

7 ingenious responses to local social and economic problems and organized 

\ 

\ 

| 

i 
| 


at low altitude, 600 m (2,000 ft) above the surface (Fig. 1). The Swiss science- -fiction 8 
4 author, Eric Von Daniken (3), in Chariots + of the Gods, : asks, ““What is wrong — ou 
with the idea that the lines were laid out to say to the gods, ‘Land here! Everything © 4 
; . has been prepared as you ordered!’ . . . What other purpose could they have 7 
_ served?” (Although Von Daniken’s ideas are not taken seriously by archeologists, ssid 
he has an amazingly large audience.) == 
With a closer look, the illusion of a sind field is dispelled by the criss-cross a 4 
‘pe Some lines are oriented at random while others radiate from a common. “24 
center like the points of a compass. Scattered among the straight lines are 
curves. Some are geometric, , but most are cartoon-like outlines of familiar 


creatures, such as a ‘monkey, a spider, or a hummingbird (Fig. 2). All are about | 


75 m-100 m across. 


oes their construction is seen to be simple. The gravel-sized stones - 

of desert pavement are covered with a dark brown coating, i.e., desert varnish 7 

(a dark brown, shiny coating of iron oxide and silica brought to the exposed 

surface of bed rock and gravel by capillarity and deposited by. eveperstien — 

in a severe desert environment). When the stones are moved aside or ame, 

the unexposed surface or light-colored protected sand is exposed. The patierns © 

are merely strips of desert from which the | pavement has been moved aside 

forming little brown ridges that outline the light colored lines and figures. A 

brief experiment suggests that a person can clear | m m’ in 10 min. Even one 

a of the larger lines, 30 m wide and 300 m long, could be constructed in 100 

-man-days-200 man-days. Archaeologists are still debating the meaning of the 

ines (12,17). Most are oriented with | Prominent stars, the moon, or the ‘sun 


at special times of the y year. All havea 


P 
J 


the time that j is required for desert varnish 1 to a the ue . 
— 500 yr-1,500 yr old. This is compatible with archaeological evidence of the 
pre-Inca civilization that inhabited this valley. It agrees with carbon-14 dating — 

_ Of wood stakes at the ends of the lines that were used to establish alignment 

» These lines and figures are not unique. There are figures outlined in the 
- southern California desert that are remarkably similar. High on a mountain 

ridge in northern Wyoming, there is outlined a lopsided wheel about 22 m (75 _ a 

_ ft) diam made of light colored limestone cobbles and small boulders that contrast . 
> with the green grass. Because some of the spokes are aligned with sunrise 
if at the ene al equinox, it also appears to have an astronomical pation 


FIG. 2. —Nazca | 100 m m (330 ft) Wise 


‘The Indians of a century ago knew that it was sacred but were ignorant of 


meaning or of who built built it. Be 


or Western RN (EuROPEAN) TECHNOLOGY 


‘The: of technology ie recsonably well 
oe First, there was the empirical or trial-and-error approach that marked 
the craftsmen of the Middle East and Egypt. Later, the Greeks added geometry ; 
and a formal logic for the behavior of physical objects. The Romans combined _ 
‘fms with their genius for organizing and communicating their empire. 
The Arabs added innovations in mathematics. 
a the Renaissance and Industrial Revolution resurrected these concepts from their 
_ suppression Caring the Dark Ages to establish the framework for engineering | 


|| 
| 
4 7 
i! 


ore 
— ae the of Civil Engineers are products 
that tradition. With this background of empiricism, keen observation, , logic, : 
and mathematics, such a genius as Karl Terzaghi was able to conceive imaginative 
and innovative solutions to problems in soil and rock that had confounded — 


engineers for milleniums. The engineering structures of this century are a 


= ago. It was conceived and constructed without two of the tools of modern | 
Logic and formal mathematics. One religious writer 
“suggests that it was built by Gods from other worlds—the Egyptians could ne 
not have undertaken such a feat. His idea is not new. The ancient Greeks _ 
attributed the large polygonal stones of the walls of Mycenae to the Cyclops— ia a 
mythical one-eyed giants. Our term ‘ ‘cyclopean masonry,”’ walls of very large es. 
_ boulders, is derived from their tradition. More than 3,000 years ago, Hebrew 


herdsmen monuments of giant stones that they found int desert 


Spanish co found flourishing civilizations in Mexico and 
& destroyed many of the native works in a passion for gold and religion 

z 4, 10). Remains of some older, dying, or dead civilizations were gutted for 
building materials or looted by souvenir seekers.§ 

_ The construction relics are still being uncovered; lost cities are coming to 

light in the jungles of Yucatan, Central America and the Andes of South America. * &g 

Their scale is staggering. What was probably the largest structure in the world, _ ‘a 
i Pyramid of Tepanapa, 1.5 times the volume of the Pyramid of Cheops, 

is yet to be fully excavated. Countless other structural works on both eiiiieiens 
only now being analyzed by engineers.§ 
- Itis the purpose of this paper to examine pre- -Columbian structures and public — 

_ works from the civil engineering. point of view. It considers their features, a 
their geographic and chronologic development, and presents hypotheses for their 
conception and construction. It explores whether they were conceived by occult 7» 

= forces, such as visitors from outer space (as suggested by Von Deniken) @, 
eer by lost mariners and early explorers from | Europe and ‘Africa who 

; never reported back to the old world, or were the products o of independent 
_ innovation. Most important, it relates the development and decline of | pre- Colum- 

bian technology to the evolution of engineering today, 


“‘elemaas or burial of the dead. Both practices are common aus not universal) 
_ thoughout the world. The need may have been sanitation, or possibly the religious © 

_ idea that both excrement and corpses have magic, malevolent power that must 
@ be kept from view. On the other hand, the Mayan farmers of Meso-America — 
7 and the Creek Indians of the Southeastern ‘United States buried their dead 
below the floors of their houses so that their spirits could rise and guide the © 


_ Where | hard soil made excavation difficult, heaps of stones substituted for 


ie, 
| 
| 
| 
| 


burial. An n eskimo grave in Greenland is a crude cairn of small boulders over 
Both shallow pits and heaps of stones were used for shelter. The early Indians 
_ of the Mesa Verde area (Southwest Colorado, about 500 A.D.) found the ground — : 
‘was warmer than the outside air in winter and cooler in summer. They aaa ; 
elongated pits about 1 m deep and 3 m-6 m diam; where the soil was loose, — 
the holes were lined with stones. Above ground, low rubble stone or laced 
7 twigs formed walls; poles, twigs, and grass formed a roof. Similar pit houses — 
are found along the north coast of Peru, dated about 3000 B.C. ‘beste 
_ The more advanced Indian civilization of the American Southwest retain ned a 
the pit house in the kiva, Fig. 3. This is a stonelined pit, about 3 m-5 m 
_ diam and 3 m-4 m deep with a roof of logs and earth. The kiva served | ‘ 


5 


a fraternal-clan retreat for the men, isolated from the bustle of daily work. 
Although the details vary, most included a stone-lined fresh air duct at the >» 
base, a draft deflector at the inlet, a raised platform soph a fire, and | a small 
=. the floor that served as a ee connection 


FT = 


li of Small Kiva, Mesa Mesa Verde, Southwes orado a. 


ivas, some below ground, but others above (where the soil is s shallow), 
a ‘The Creek Indians also built earth covered ceremonial or clan lodges with 
_ provisions for a ‘fire and a small floor hole. These date between 500 A.D. — 
_ These excavations possibly reflect the instinctive feeling of security that a 
child « exhibits when he retreats to a corner - when under duress, or the dugouts 
“made by older children to escape the adult world. The pit is not limited a 
_ Primitive societies or children; the United States defense headquarters is far ; 
4 underground i inthe Colorado mountains. 


hae heaps of earth and —_ of stones were constructed ape men in all 


; 


e Trash : ee are one form. Despite the romantic iat of historians, some 
American Indian tribes merely dumped their trash, garbage, broken pottery, 


oy from about 300 B.C. are found along the sea coast in a Peru. Similar shell prece 
* occur along the South Atlantic and Florida Coasts in the United States dating — 
- from about 2000 B.C. Some became mounds as high as 6 m (20 ft). Although ra 
ag they are now rich sources for archeologists, they must have bred flies a ae 
rodents, and exuded rich odors when in use. 
a Burial mounds, , constructed between 800 B.C. and about 150 A.D. by Woodland — : 
: Pe Indians, are widely distributed in the Eastern United States. Many are from > 
es 3 m-10 m (10 ft-33 ft) high with side slopes typically 30°-35° from the horizontal — 
and have round tops. Most appear to have been built during a relatively short 


(a). —Etowsh Mound, Cartersville, Ga. 


vid time period, representing only a generation or two, based on a comparison 
of possessions th that accompanied | the burials and on very limited carbon- 14 dating, = 
_ Long- -term use of a site produced a succession of mounds instead of a ‘single, © a / 
* growing burial. Burial mounds are also found in the coastal settlements of South i 
America, dating from about SOOB.C.(11), cont 
Platform Mounds.—Platform mounds for worship ceremonies or bases of 
_ temples appeared i in a the Eastern United States about 700 A. D. Similar structures ~e 
2a mounds, built by ‘filling over t old pit and wall houses with soil and 
_ debris are found on the north and central coast of Peru (11). The debris associated 


> 
with them dates about 2UUU B.C. by Carbon-14. Ihe mounds were constructe 

materials with uniformly sloping sides and flat tops. Most 


ru 
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in stages over a ‘ile period of time, eventually becoming the largest structures 


__ The Etowah group, near Cartersville, Ga., 35 miles north- northwest of Atlanta. ‘ 
a were built on a low terrace on the right bank of the Etowah River. It consists 
of three mounds, a village area, and a 1.3 m (4 ft) deep perimeter ditch enclosing 
a circular area about 450 m (1,500 ft) diam. The smallest mound has been 
completely excavated . It was constructed in stages between about 700 A.D. a 
and 1200 A.D. and consists of a series of concentric truncated square pyramids 
with side slopes of 33°-37° with the horizontal. The top was flat and was similar 
in size for each stage. On each top was a simple wood temple with walls of _ 
vertical poles set in sockets in the ground. The bark-stained sockets remain 
and mark each platform level. At the death of a chief or priest, marking the 


— i end of a period of unified leadership, the wood structure was destroyed, and 


the leader was | buried in the i¢ mound with weapons, jewels, and vestments. A 
new: platform and temple | were then constructed 1 m- 3m 3 -ft-10 ft) above | 


| 


the earlier one. The stages in construction are reflected in in different soils and | ; 
_ The largest mound of the Etowah group, Fig. - 4(a), contains about 110,000 

om or 150,000 cu 1 yd of soil. It has not been excavated. There are remains 

of a wood temple or council house on top. It appears to have been built in 
the same stage construction as the smallest mound. 

_ Monks Mound of the Cahokia group, adjacent to ‘a small creek about — 

miles east of St. Louis is the largest i in the eastern United States. It is apie 


i it was built in 1 about 14 stages, commencing about 900 A. D. and ending 
about 1150 A.D. (16). The stages range in thickness between 1 m (3 ft) and 
5 m (17 ft). The total fill volume is about 750,000 m’ or 1,000,000 cu yd. 


urs Protective Mounds. oe ew mounds: are in the f orm of continuous embankments _ 


| 
| 
- 
| 
near Cincinnati, Ohio, a perimeter mound 6 km (3.5 miles) long 


~ throughout the eastern United States. Most are low, 1 m-3 m (3 ft-10 ft) — 
and cover areas 10 m to more than 100 m (33 ft to more than 330 ft) diam. ee, 
_ Animals, such as the panther, bear, deer, buffalo, turtle, beaver, and dog — 
"predominate, but birds are also found. Most are on bluffs above rivers and 
are associated with small burial mounds. Those excavated appear to have been a 
built in a 
burial ‘they appear to be clan symbols | or ceremonial structures. 
_ The Serpent Mound of southeastern Ohio, Fig. 4b and 4c is the largest and 


best known of the effigy mounds in the United States. It lies on a bluff Ce a 

a small river. It is made of gravelly clay, following a pattern laid out with a 
- arow of stones and ashes after removing the topsoil. The approximate volume 

is 2,000 or 2,500 cuyd. bs te 
Mees Construction.—Those mounds that have been excavated or core drilled _ 7 

all exhibit evidence of similar construction. They were placed in level layers 

2 cm-8 cm or | in.-3 in. thick. In some cases the edges of the layers were 
_ confined by logs secured by stakes: simple forms. These forms were removed _ 


and reset as the work progressed, but occasionally left to rot in place as shown © 


= by remains of bark and rotted wood. Side slopes were typically 33°-37° with 
horizontal with no obvious correlation height, soil texture, and 


40 ha (100 acre) segment ound varies from 
— - 3 m-5 m (10 ft-16 ft) high with side slopes of 30°-40°. Other fort-like or stockade 7 
5 mounds have been found in the upper Mississippi and Ohio River area. 
: _ Effigy Mounds.—A number of mounds are shaped like animals. Most are __ 
4 
Bet... 

| | 
| \ 


eroded gully was sometimes filled with stones. A. 4) 

The tools for construction have been pieced together from fr ragments etsitialin, 

in the mounds. The soil was loosened with a pointed digging stick, hardened | 
=z fire. The same sticks were also used for planting corn. I have seen them 


q use today in remote rural Mexico, Honduras, Peru, and Bolivia. Excavation . 


was by hoes made from shoulder blades of deer lashed to poles. ‘The earth 7 
sh carried in baskets or leather bags” on the backs of the workers, about — 
kg or 45 Ib at a time. wove 
_ Experiments by various investigators, cited by Erasmus (6) suggest that a 
team of one digger and three basket carriers could excavate 3 m’* (4 cu yd) 
and move and place the soil 200 m (650 ft) from the source in one day. At 
this rate, the largest Etowah mound required about 100,000 ‘man-days and the © 
largest Cahokia mound 1,000,000 man-days. The Incas of South America were 
required to devote about 90 days / yr to community labor. Reed et al. (16) —¥ 


_ that the North American Indian may have contributed 40 days. At this rate, ‘ 
_Etowah required 2500 man-years and Cahokia 25,000 man-years. Both were j 
built over T periods of about 250 years. Thus the Fequired average annual labor 


however, ¥ was built in 14 stages averaging 18 years (probably 14 generations), = 

thus Cahokia required 1,800 man-years/ generation. On the same basis Etowah = 

Peo only 200. The large volumes are not unreasonable for the modest labor 

forces available, when | the reli and leaders the 


vwabes from the burial mound. In Meso-America and Peru the temple inne 
_ developed independently. The earliest were flat platforms of earth and boulders 
us faced with walls of rubble ‘masonry and topped with flat stones (11). On top 
‘ were stone monoliths and remains of small buildings that probably were worship _ 
related. The earliest appeared about 1800 B.C. in Peru and Mexico. 
__ During later periods, some of the platforms were enlarged by building ol 
walls around them, and then raised in increments of 3 m-5 m (10 ft-16 ft). 
# he platform became a series of terraces or a step pyramid, with very steep 7 
sides, and a level top which was the focus of religious rites. AR 
Still later, the platforms were augmented by faces of mites adobe brick, 
: and later of trimmed stone, battered at slopes of 60°-70° with the horizontal. ; 
These truncated pyramids first appeared about 1800 B.C. in Peru. Their construc- 
= tion climaxed with the Pyramid of the Sun at Pachamac, about 30 km (18 miles) | 7 
. oe of Lima, between about 200 A.D. and 1000 A.D. The pyramid is of Bi 
stones and adobe brick, about 220 m (700 ft) square, flat topped, and — 
7 _ The largest earth-stone pyramids were constructed in the Valley of Mexico. 
_ At Cholula, about 50 km (30 miles) east of Mexico City is the Pyramid of 
Tepenapa. It is possibly the largest structure built before the 20th century, 
constructed between 350 A.D. and 600 A.D. The base is 402 m ? (132 sq ft) 


LATFORMS AND PYRAMIDS ry 
| 
4 
it 1S Mm ign. Ougn it iS NOt as as 


th flat ™ an was replaced by a Roman Catholic Church which is still 4 


pyramid is constructed of light- -colored tuff boulders about (0. m 


. the heart of the pyramid have found five (or possibly seven) smaller concentric 
. e pyramids with plastered faces. A painting of a giant grasshopper on one face 
Deep in the jungles of northern Guatemala, along the humid Gulf Coast of 
southern Mexico and in the tangled underbrush of the Yucatan Peninsula are 
_the most magnificant structures of pre-Columbian North America, i.e., the Mayan — 
temple pyramids (14,20). Commencing about 500 B.C. and ending between 1200 
4 A.D. and 1400 A.D., complex religious centers resembling cities were constructed 
in the tropical lowlands. A wide variety of massive structures remain, many 
still choked by the jungle and awaiting archeologic study. Enough have = 
_ examined to demonstrate the evolution of advanced construction technology 


~ over a wide area (Fig. 6). Although each has its imposing structures, Tikal, 


GREAT PYRAMID: 2,600,000 


3,500,000 Yb? 
TEPANAPA: 3,800, 000 M 
M 
’ 


70 


5.—Comparison « of of Tepanape, Cholula, Mexico with Pyramid of 
Guatemala represents one of the earliest and longest sived, of t these areas (1,2, 2 q 
“a rs 1, Fig. 7 was completed about 700 A. D. near r the | climax of Tikal’s 
of its nine terraces 65°. It has a narrow truncated top, 30 m m (100 ft) 8 pag 2 
the Plaza. On top is a small temple with a soaring crest or comb rising another 
14 m (45 ft). Although other temple pyramids at Tikal are higher, and many 
in Yucatan, such as at Uxmal and Chichen Itza, are broader, Temple I remains a 
an impressive example of the best Mayan construction. A Stair leads up to 
‘the temple, so. ‘steep that one climbs it today with the aid - of a chain. The si 7 
pyramid surface i is of well-fitted soft limestone blocks laid in courses. __ an 7 
_ Surface repair of the Tikal pyramids reveals that they were built of boulder-sized 
soft limestone rubble in a matrix of limcstone chips and cemented with lime 
= mortar. Tunnels into the pyramids reveal smaller, older pyramids within (1,2). 


Large carvings on the older pyramid surfaces, hidden by the larger, a 
— construction, can be ‘seen in | exploration | adits, excavated by ar 

_ The quarries from which the soft limestone blocks were cut and the rubble m a 


‘obtained are ‘still i in une, furnishing fresh fresh blocks for restoration of the structur 


3,770,000 m’ or 5,000,000 cu yd is approximately 1.5 times greater 
q 
| 
7 
4 
i 
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‘Fig. 9). Deep, narrow were cut, using tools of chert" or flint. T 

; _ rock was then broken apart by prying or g or wedging i in the slots. caer 

_ Adobe is a universal construction material used throughout the Americas as 
_ as in the Old World. Technically it is unburned soil (clay) brick, but the 
has been applied to other forms of dried soil construction. ‘bes 


— Sun- -dried brick is a form of artificial stone used in arid regions where local 
a” 4 rock i is difficult to shape or is absent. Based on data from sites in the western 


— 
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YUCATAN PENINSULA } 


PITTED LIMESTONE SURFACE 
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_ 
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a GUATEMALA city . 
requirements are sufficient cohesion to provide dry and enough 
inert filler such as sand, fine gravel, and even straw, grass, or twigs to minimize _ 
shrinkage. Based on unpublished tests by the writer, unconfined compressive | 
strengths range from 100 kN /m? “nearly 700 kN/m/? (15 Ib/sq in.-nearly 1000 
in.). The tensile strength is typically 1/6 the unconfined compressive 
“a strength. Adobe resembles weak, plain concrete in its strength and ability to | 
absorb repeated strain and tension. It resists static loads in walls up to 15 _ 
(50 ft) high but crumbles under the dynamic loads of earthquakes. 
4 _ Adobe balls and cones were used as masonry units in coastal Peru as early — 
as 1000 B. C. Rectangular adobe blocks appeared shortly afterwards sometimes "7 


interlayered with stones. Th he adobe brick ‘structure in the ‘Americas: 


— 
> 

| 
| 
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The Pueblo Indians of the United States ‘southwest used both adobe and 
"stone for construction depending on material availability. The Casa Grande of 
Arizona was built between 1350 A.D. and 1400 A.D. as a part of a village 
> compound of the Hohokam, whose descendents may be the Pima Indians (9, 25). 
Lenticular pancakes, typically 0.35 (1.3 cu ft) diam of firm silty were 
aid in very irregular noncontinuous courses and level tiers, 1.2 m (4 ft) wide 
and 0.5 m-0.75 m (1.5 ft-2.5 ft) high. Each tier was shaped to make flat vertical — 
i surfaces and then allowed to dry before adding the next, producing cold joints | 
in which hand prints are still visible. Although all the soils in the walls were 
~ from the site, there are significant variations between tiers,  feflecting — 
: the different calcium carbonate percentages in the soil horizons excavated. The 


insides of the walls (and possibly the Lege were plastered with adobe and 


ace 


> 


* 7.—Temple | Tikal, Guatemala 


/ 50,000 man-hours of labor, based on the excavation, mixing, and placing of 
each soil pancake. A village of 20 families could have built it in a year or 
two after allowing for their farming and hunting. Similar adobe — 
has been found in other southwestern United States pueblos. beet Sy 

The temple platforms and pyramids are massive, requiring much labor, but 
limited ingenuity. Walled structures, which also evolved from simple heaps <= 
stones, required advanced technology to provide stability against overturning 
by wind, earthquakes, and security against enemy attack. 

_ The first walls in wed Americas were of cobbles, stone slabs, or adobe ponlll 


| 
| 
= 
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= about 2 2500 B.C. Similar cobble and mud walls w were built along the i celtic 


cliffs o of the ‘southwestern United States to enhance the shelter provided by i 


the sandstone overhangs. 
‘Wall construction evolved rapidly in the southwestern United States » using 
3 sandstone flags about 30 cm? (1 sq ft) by 3 cm-5 cm thick, laid in courses, 
with the stones in successive courses overlapping to form a bond. The earliest 
used flags. for uniform size. In later construction, 


FIG. 8. —Stone Face 2.5 m (8 | in Tunnel into Pyramid 


outer face of uniform thickness. A good example is Pueblo Bonito, Chaco Canyon, - 
_N.M., built about 1000 A.D. The higher walls consisted of outer casings of 
bonded and faced flat stones with cores } Of Tough flat stones crudely interlayered 


of the wall is reflected in n the greater width at the base, Im n (3. 3 ft) pontarsed 
a3 to 0.3 m (1 4 at a height of about 8 m (25 ft). The eenan nee laid with 
: 


tid” 
j 
2 
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stability was provided by cross walls; the rooms were typically 2 m-4 m wide 


x Coursed } masonry walls were used n more sparingly i in Meso- American construc 


masonry ‘of tuff blocks, 25 « cm x 50 cm, on yn the ‘outer rt face. There a are a | few 
“Mayan (700 A.D.-1200 A.D.) walls of soft limestone, generally laid A al 
mortar. The blocks were of varying sizes, from 20 cm to about 50 cm high, me 

i. but uniform in each course. © The — ae to width ratio of the blocks _ 


Masonry walls were constructed extensively in the Peru-Bolivia 


FIG. 9 Quarry in Soft Limestone, Tia, 6 
= 
of unequal height, 10. notching of the corners 
to fit successively higher courses s together. Th The fitting was very precise and 
most of the walls were laid without mortar. At Tihuanaco, Boliva, some of 7 
the stones were held together by I-shaped copper cramps, forced into dumbbell 4 
Shaped notches cut between adjoining stones. The use of the wall was reflected 
in the quality of construction in the later Inca masonry, 1200 A.D.-1500 A.D. e 
Terrace walls in rural areas were built of roughly trimmed rectangular blocks 
_ varying sizes, laid in courses that were not necessarily I level nor of F equal a, 
height. Military structures such as the Fort of Sacsahuaman near Cuzco o (Fig. 
10) were made of huge, more or less rectangular blocks, laid in irregular courses, 
with blocks of unequal notched to fit The fitting of these umesione 
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all. The : surfa ace texture of the stones suggests that they were shaped by poundieg 


™= small hard stones, , followed by surface ae stones and sand. 


4 upper blocks. The construction of Sacsahuaman, completed about 1480 A.D., 
required 50 yr utilizing tax labor, about 120 days/yr/person, and a legendary 
labor force of more than 20,000. One large block, the ‘‘Weary Stone,’” wel 
never trimmed or lifted in place. According to Inca legend (10), it was so =£ 
_ after being dragged to the fort site that it wept blood from two holes in one 
corner. More likely it was the laborers who wept—3,000 were killed when 7 
ee fell from a slope during hauling. ’ 


be 


, 
Phy: 


¥ FIG. 10.—inca Wall of Blocks: Fortress of Saceshuaman, Cusco, 


Peru 
_ Because of the irregular coursed, notched corners, and close fitting, the walls 
have demonstrated unusual resistance to earthquakes compared with the uniform 
- coursed, mortared walls built in the same areas by the Spanish conquerors. 
Some historians have speculated that the notched form represents a conscious 
effort for greater stability. However, if this were true, it would be used on 


the most important structures, i.¢., religious and governmental buildings. — 


the most sacred buildings in Cuzco,’ such | ‘as the Convent. of the Virgins o of 

‘the Sun and the Temple of the Sun. The writer believes that the notched masonry * i 
a. a compromise between the ordered appearance of uniform coursed j 
“Masonry and the labor required to produce uniformity from rock that did not 


- - dragged from the nearby quarry by hundreds of men utilizing twisted vine ropes _ ' 
either wood skids or rollers Earth ramps were constricted to place th 
| 
| 


Mortar. —Mortar i is to bed and cement t the units of a 
_ and to plaster over a rough wall to provide a smooth surface, often for decorating. 

A low to moderate plasticity clay was the most widely used material throughout 
the Americas. It is durable if protected from wetting, 
the masonry units were so irregular that thick mortar was required, 

_ shrinkage became a problem. Pebbles and stone chips were forced in the mortar 4 
= _ joint to steady the stones, reduce the mortar volume, and minimize shrinkage. 
_ This practice was widespread in the United States southwest and in Meso-America. 
Ps - In the United States southwest, caliche, a calcaerous silt or clay, is abundant 

a ‘Rear the ground surface and has been employed for mortar probably because 

availability. It has the € property of hardening somewhat, in alternating 


-_damp-and- -dry conditions, and g slightly water 


_ The Mayans discovered ‘that burned limestone would reharden slowly into 
soft rock upon wetting and subsequent exposure to air. Lime was produced 
limestone on firewood | racks (20). The mortar was: a mixture 


‘not appear in the Old World until a few hundred yenss B.C. The Mayans - 
‘& possibly as early as 500 B.C. for masonry construction as well as plaster. __ ~ 
Composite Construction.—The Mayans made extensive use | of massive com- 
_ posite wall construction. The interior of the wall was coarse, rubble masonry, 


cemented with clay or lime mortar and chips of stone into a homogeneous 
mass. A veneer of stone masonry, precisely fitted and coursed, covered the - 


q 
— 
“Sullace, Providing te WUSION Of maSOMry OIOCK Wall (rig. 11). practice 
fod can be seen in early structures at Tikal, about 200 B.C. and some of the latest [im 


a Chichen Itza and Uxmal, 1200 A.D. to 1400 A.D. 


wood-soil-plaster composite. The isa of wood poles Scm-10 
em (2 in.-4 in.) diam spaced 1 m-2 m (3 ft-7 ft) apart. Between them is 0 
- ‘matting or wattle made of parallel (but sometimes agp woven) a 


clay mortar is applied to both ‘sides to a wall 8 
cm or (3 in.-6 in.) thick. Such a contemporary house in Yucatan is seen in 
Fig. 12. The deteriorating clay exposes the wattle . The 


“Nunnery” (although its use is not known) of Uxmal. ¥ 
-‘Similar b houses can be seen today | throughout the ‘Mayan areas, in other pect: 


vo 


12. Wattle tle and Daub House, ne House, near Ch Cc ichen Itza, Yucatan at 


of Meso-America (and in Kenya and adjacent areas of East Africa). eae 
_ Retaining Walls.—Retaining walls served two purposes in the pre-Columbian ~ aa 
"Americas: (1) Agricultural terraces; and (2) support of steep slopes. The ,4 
7 - civilizations of the Andes Highlands of Bolivia and Peru built extensive terrace — 2 

_ walls to make level land on steep slopes and to conserve rur runoff. Most were — 

of untrimmed, small boulders laid in crude courses (Fig. 13). rem ates thaw? 

- to support steep slopes were a part of the structures of Monte Alban ~ 
in Mexico and Mayan Copan in Honduras. They appear to be incidental to © 
the structures to which they are attached and not built specifically for retaining. — 7 
a Pueblo Bonito, Chaco Canyon, N.M.., is built at the toe of a sandstone cliff. 7 


é Soft shaley layers near the base of the cliff were weathering and eroding, allowing 


the massive sandstone above to tilt outward from a vertical joint crack (Fig. ; 
Fenner The pueblo cae constructed a wall at the toe of the cliff to ouenail 
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weathering rock and to farther about 1000 

_ A.D. For more than 800 yr it functioned. In 1941 the cliff collapsed demolishing _ 

the wall and a portion of the pueblo. 


"Foundations in the Americas were usually simple. Although the bottom of 


a high masonry wall was sometimes wider than the top, there was no enlarging - 


of the base in contact with the ground to make a footing. At the sites of the 
larger Mayan and Mexican structures, rock was usually shallow; therefore, _ 
foundation bearing was no problem. The same is probably true of the — 
South American structures; however, foundations are not re in archeo- _ 


FIG. 13. 400 y yr ye Ago but Still Intact, 


logical reports. Some of the precisely fitted Inca masonry walls rest on crudely 7 
laid rubble below the ancient ground level. This suggests that a trench was 
dug, filled with stones, and the wall placed on top. The trench width and depth, 
however, are not significantly greater than the wall thickness. 


2 Wood columns in most North American construction were placed in holes _ 


“about 1 m deep, and backfilled with soil to provide lateral support during 


_ construction and (probably without intent) some vertical support for the eventual — - 

~ Joad. Stone columns of the Mayans and Mixtecs usually sat on the stone platform _ 
of the temple complex; separate foundations were not required. = 
‘The Pueblo Indians in the vicinity of Chaco Canyon, N.M. used prefabricated z ‘ 


“stone. disks at the bases of the columns for the Great Kivas. Three or four, | 
~ about 1 m diam and 4 cm thick, were stacked like pancakes under each | column. : 
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Wood piles served as both foundations and ‘unsupported Cbladens in areas 
os wet soil. The Spanish explorers found Indians in Lake Maricaibo on the 
‘north coast of South America living in pile-supported houses over the water. — 
= named the area Venezuela: ‘‘Little Venice.’’ There are similar pile supported © 
g in the area today; the only change is the substitution of wood-plank | 
- siding and corrugated iron roofs for the indigeneous pole ; and thatch. Pam oer 
Similar wood pile supports were used in the Valley of Mexico where ‘the ; 
ye built Tenochtitlan, the first Mexico City, in a marshy arm of Lake ff 
- Texcoco (14). Excavations for modern construction have uncovered wood pile — 
- foundations for houses, and to retain soil fill for roadways. Although there | 
no Tecords of how the Piles were driven, the writer observed m men driving» 
They 
“raised a hammer 25 cm in) diam and | 60 cm ft) long, by 


FIG. 14. —Re Erosion Wall Pueblo Bonito, Chaco Canyon, New 


uaa and it 60 cm on top of the 10 cm diam wood pole 
while standing on a wood fastened to driven piles. 


"Transportation Routes 

Transport of people developed late in the pre- -~Columbian 
Americas. There is evidence of loosely organized trade in the finding of native - 

~ copper in mounds in in Georgia, which probably came from the only known source 

_ Similar trade in mica and flint between different areas of the eastern United am 
; ‘States i is inferred from grave artifacts. The trade routes, which were also used 
: _by the early European explorers, were water courses and loose networks of _ 
foot paths, 0.6 m wide, which probably evolved from hunting trails. There 

, was little need for more. The only beasts. of burden were men, women, and — 


occasionally dogs. The largest boats were only canoes or rafts. In the Andes = 
of South America, the llama, a mini-camel, still carries loads of about 50 kg ‘ 
Ib) but has never pulled load. Coastal ocean commerce developed wail 


— 
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= in the Yucatan area and to a small extent along t 


Coast of the United States and Canada. woe 
p _ Surprisingly, systems of wider ‘‘highways”’ are found in many areas. Because 
these routes focus on major settlements and religious centers, archeologists 
~— speculate that they were required for pr processions or perhaps for armies moving - 
in compact units. ply 16? doen pow 
_ The best known are the Inca roads in Peru and Ecuador (22,23). One 5, 00 
km or 3,200 miles long followed the high valleys in the Andes, from the 
administrative capital at Cuzco, Peru northward into Ecuador and southward — 
into Chile. A second road followed the coast linking the irrigated valley _ ; 
civilizations across the sandy deserts. They were connected by east-west roads. 4q 
In the plains they were merely trampled yma, marked by rows of stones 3 


15.—Pile Driving by Hand, Acapulco, Mexico, 1949 x. 
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m6 m apart, and swept me of sand regularly. The roads were —— 


retaining walls, paving, stone stairways up steep aan (Fig. 16) and parallel 
_ marker stones. Suspension bridges spanned narrow canyons and tunnels pierced — 
_ steep obstructions. At intervals of 7.5 km (4.5 miles) there were distance markers. 
Relay stations were built 2.5 km (1.5 miles) apart for messengers who ran the | 
7 interval, passed messages and ‘perishable food to the next runner, _and then 


" ‘miles), there were rest | stations for official travelers port supply mage ll 
P 4 yer small military parties. Possibly the litter parties that carried the nobility — 
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only occasional llama pack trains. Parts of the Inca 
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7 in use, while aun: lie in disuse because of the steepness | of the pathways 
or because they are lost under the desert sands. al 
‘The Mayans likewise had an extensive , road system called ‘‘Sacboebs”’ 
- connecting many y of their cities and religious centers (Fig. 6). Many can be 
traced today as lines of different growth in the jungles of Guatemala and the — 
scrubby forests of Yucatan. The ground surface in much of the Mayan ting 
Ba extremely rough, deeply pitted from the solution of the soft limestone, making — 


constructed with level surfaces and in straight lines. 
_ The causeways were built of small limestone boulders, with limestone 
chips and limestone sand to form a smooth surface, anticipating the ideas of 


McAdam, the ‘Scottish engineer, centuries later and an ocean stone 


to a a roed suggested a a ‘roller to Mexican geotechnical engineers (18). They found 
, that 15 men could easily push it across freshly placed stones and compact 
the surface. (Some archeologists are skeptical of this explanation of the stone’s 
use.) Most of the ‘causeways have become overgrown with vegetation after 


of the United States National Service, studying 2 air 
graphs of Chaco Canyon, N.M., for prehistoric irrigation channels, were surprised 

by faint linear features radiating from some of the ancient pueblos (5). The 

lines go ae across topographic features but are best defined on the uplands 


| 
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= 
enturies Of disuse. A lew Seco ave Dee Ored DY the Mexican | 
‘pathways 6 m-9 m (22 ft-30 ft) wide. Across bed rock, they were 


or “roads’’ are still evident (Fig. 18). Their use is not 
; archeologists. speculate that they were paths from the cultivated fields and springs if : 


the rock. ‘Some terminated in staircases down the cliffs. Across sandy ground 
the pathway v was a shallow excavation. ‘Despite | seven centuries of neglect, these — 
a Although much of the transport of fms and commodities into the Aztec 
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} 
_ marshes connecting it with the mainland (Fig. 19). Closely spaced pole piles — 
feemes parallel retaining walls (Fig. - 20). Possibly reeds or canes were interwoven 
with the piles to produce a tighter wall with fewer piles. Earth fill from the 
adjacent marsh was placed between the pile walls forming an earth causeway 
with a canal for canoes on each side. There were occasional openings in the 
& causeways to allow water circulation and passageway for boats. According 7 


for military defense (4). "Cortes, his soldiers, and Indian allies were nearly rapped 
by the Aztecs who opened the bridges during his first sortie into Ge city. 
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THERE! WERE GIANTS 


Water and soil for agriculture are twin essentials to community erm, 
& many regions of the Americas, highly seasonal rainfall was unfavorable to { 
_ continuous settlement and agriculture. In coastal Peru, whees the sea was a 
source of food, lack of any rain made life nearly impossible. In mountainous 
areas, steep slopes eroded and slumped when the permanent vegetation was 
= by farming. Despite shortcomings < of climate, some of the - most 
advanced civilizations in the Americas in areas. | Water and 
_ Terracing, in steep terrain, with the aid of bos rubble a retained both 
“soil and runoff. There are. of such terraces: in the south- 


and 


FIG. 18. —Chaco Canyon Road or Path Sandstone 


southeast of Mexico City, terraces have been found in cultures dating jon 

_ In South America, terracing in the Andean highlands and in coastal settlements 
probably began about 500 B.C. (11,23). Because of the continuity of commonity ; 
life in the Andean highlands, the terrace systems were improved and enlarged — 4 
until the end of the Inca empire in the sixteenth century. There are many ry 
that remain in use today. Many others are in good physical condition despite ; 
four centuries of neglect, such as at Machu Pichu (Fig. 13). Others have crumbled; 
their remains are a testimony to better water and soil management than is practiced i) 
today in the same region. The terraces in some areas are as high as 2 m or : 


ft “constructed on slopes: age -20°-30° or 2. 7(H)-1(V)-1. 7(H)-1(V) but 
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with soil, underdrain of sand over gravel to prevent 
“soil saturation (and also to prevent hydrostatic pressure from overturning the _ 
~ wall) (15). ‘Unfortunately, the technical knowledge at and the political s system that ie 
s = made it possible to build and maintain these terraces were destroyed by the i a 
_ Water Diversion Canals.—Diversion of water from streams and springs to 
farm lands, and to terraces by canals or ditches, represents a logical step ra 
terracing. Unlined ditches have been found in the southwest United States bait 
by the -Hohokum n 2,000 years ago (9). The builders diverted streams 
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19.—Causeways and Dike, Aztec Tenochtitlan (Mexico City), 


d using temporary dams of woven fiber matting supported by short poles driven 4 
into the stream bed. In Chaco Canyon, N.M., there are stone rubble structures 
= in dry gullies to divert thunderstorm runoff into cultivated fields (24). i. q 
> % In Coastal Peru, where there is no rainfall, agriculture 2,500 years ago depended = 
on diverting water from the rivers which flowed from the mountains into the 
flat terraces alongside the arid flood plain. In the Andean highlands with seasonal 2 
rainfall, multiple use was made of the occasional dry weather springs. The 


water wv was collected i in n masonry channels, diverted first to fountains, and baths, 


| 

| 

| 
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trenenes leading tO terraced j 


_-- is a combination terrace and canal system (Fig. 22) (7). Warm water, ve 
-Tich in magnesium and calcium bicarbonates, flows from several small springs _ 
ea a steep hillside. The carbon dioxide is released in the air, precipitating the be. . 
minerals as travertine, a soft limestone. Ordinarily, waters cannot be 


The Hierve el Agua 20 east Milta in Mexico’ s Oaxaca 


the soft deposits were recut forming new channels on top of the old. Today, 

_ ‘more than 2,000 years after the system was begun, the minicanals are perched © 
g travertine ridges 1/2 m-1 m high and 1/4 m-1/2 m wide. These ridges — 
baw form terraces on \ the steep | hillside At intervals along the 
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“SECTION A-A ROAD CAUSEWAY SECTION B-B DIKE 


20. 0.—Cross Section of Road Causeway and | Dikes of (4) 


_of Belize and Guatemala and the Yucatan Peninsula, there is 5 little s stream flow 
despite abundant rainfall. The runoff disappears into cavities in the soft, porous rf 
limestone. The occasional large open sinkholes, or cenotes, were the only natural — 

- source of water. These became the focal points for many Mayan Temple cities. 

‘The cenote of Chichen- Itza, ‘Yucatan was a ‘sacred well. When the 
were sacrificed into itto bring the water back, 
x ‘At other Mayan Temple cities artifical reservoirs were constructed in the 

_ limestone to store surface runoff. Some were probably old solution pits in the 

rock, others were carved in the soft limestone. The pores of the rock were 
sealed with | lime- dust mortar and occasionally clay. Some at Uxmal 


masonry reservoir in a runoff-collecting system on the sandstone a 
top. Storm water was diverted ‘into | it by | ditch, 


builders constructed an elaborate system of minicanals using the precipitating 
travertine (Fig. 23). The lengthy network prolonged the air exposure, enhanced [ie 
_ precinitation and ed ae he minera onten Vhen he ine hii 
| 
4 
: 7 es water was handdipped from these to water the plants on the terraces (7). 
_ i” Other small canals with small cisterns for dipping water are found throughout 
= = — ancient Mexico; some were derived from springs or small streams but a few 
Ff 
| 
: /_catchments (Fig. 25). At Tikal, there were a number of open reservoirs with _ 


of the temple city y complexes of Meso-America and i in n the Andean es 


stone-slab culvert systems, draining to nearby gullies and streams (13). Grooved 
stone bars were used as runoff channels at Copan. Similar grooved stone water | 
conduits can be seen on the Andean highlands. It is possible that forse - 
v were joined to form closed | water pipes. 


used for drainage human waste. If: so, , they ‘accummulated decompos- 
ing solids during periods of low rainfall because there are no springs on the 
mountain top and rainfall is highly seasonal. ae 


FIG. 21. 


—Spring, Baths, and Water Diversion System at Tambumachy, Cuzco, Peru 


7 At Dainzu, Southeast of Oaxaca, Mexico, there are extensive conduit systems nae 
of grooved stones and stone slab culverts. There are also burned clay sad 
about 10 cm OD (Fig. 27). They brought spring water into the temple complex — 
A ribald sense of humor can be seen in roof drains for a temple in “‘Old 7 
Chichen,” a part of Chichen Itza. The roof downspouts are in the the form of 


when it rains, they discharge in a realistic manner. 


Valley of Mexico.—The valley of Mexico (Fig. 19) was once occupied by — 
a shallow lake between the surrounding volcanos. The early inhabitants excavated 


canals in the marshy lake margins and piled the excavated peaty, silty clays 
to form low rectangular artifical island called chinimpas (14). Several 
bos per year were possible in the warm sunshine | with the organic soils and 


ree adjacent source of water. Human wastes were collected re 


| 
4, 
= | 


F settlements to enhance the the natural fertility. The present “Floating Gardens of . 
_ Xochimilco’’ (the rem remnant of an ancient south arm of the lake) still produce © 
flowers and vegetables, although they are now mostly a tourist attraction. —__ 
a The Aztec capital of Tenochtitlan (now Mexico City) in the western marshes — 
of Lake Texcoco was an adaptation of the chinimpas (plus some older islands) 
- urban _ development. The island city was connected to ‘the mainland by 


large fresh water inflow from springs and streams compared with its surface 
area from which water evaporated. By limiting the area of the cross-canal — 
L: to the more saline larger eastern ml of t the lake, the A Aztecs were able le to 


_ maintain freshwater around the city. 


: had no. outlet _ The causeway isolated. the western part of the lake, with its 


aa 
___ Windstorms generated sufficient waves and water pileup on the main lake 


: that the causeways were overtopped and parts of the Aztec City flooded on > 
a Sow occasions. Shortly before Cortez’s conquest, the Aztecs constructed nll 


to a shallow saltpond.; 
a The pre-Columbian civil works » with their unusual variety and size, are a 
tribute to those ancient engineers who conceived their designs , and oe: i 
their construction. The civilizations are dormant and the names of those a 


of the old causeway water control system, and Lake Te excoco has shrunk 
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2 - that can be seen in the Old World (although this conclusion may be partially 


wi The: first test of the value of engineering i: is $ its contribution to society’ s ty yon 
physical needs, i.e., shelter, water, food, communication, and defense. Enough | 
examples of engineering works that meet each need remain unswallowed by 
jungle or unmelted by weather to demonstrate that they rival comparable structures 
of ae and as Middle East of 2, 500 or 000 yt ago. / As will be considered Sof 


Structures that expose man’s religious and spiritual feelings also fulfill = 
4 fundamental need. In size and grandeur the temples and pyramids of the ¢ Americans ~ 
compare with those of the fps ys Babylonians, an and Assyrians. In construction, 


ee 


23. of Hierva el Agua Canal ‘and Terrace System (7) 

_ A third test of technology is the adaptation to available materials. In this 


_ the pre-Columbian Americans excelled. Because long, heavy timber was scarce 
inthe Andean highlands, the Valley of Mexico, Yucatan, and the arid southwestern 


_ United States, local stones and soil were exploited. As today, the soil was 


dirt cheap—a versatile, strong, durable material if handled correctly, a pile 
of mud if not. Even a child knows this instinctively and empirically; some 

Soil was made into walls, and was used for the mortar for walls of stone, 
and the plaster that smoothed the surfaces. Reinforced with sticks and reeds, 
‘it became floors, roofs, and the walls of huts in areas where soil was too 
scarce for adobe block construction, or where moisture softened the blocks. ll 


responsible for these a 
| 
| 
| | 
| = © geotechnical engineer Can ting a Gceper mcaning in IS USC O 
gg earthy materials. Many religions teach that man becomes a part of the earth | 
7 __ when he dies, therefore, he must be made of earth. The ceremonial underground . 


it or kiva of the ‘Southwest Indian, is thought to represent the womb of the a 

a from which man springs. The civilized child sometimes finds comfort 
in a close space among rocks, a shallow pit in the ground, or a dark on 
of a room. The earth, or rock mound that supports a temple, reflects man’ 4 
(and all life’s) dependence on the earth. The pit grave and the burial mound 

signify his return on death. ys Pad) a 

q Thus, the geotechnical engineer can feel that he is working with something _ 
more fundamental than a mere construction material—he builds with life itself! 


— Detail Raised Min Minicanal Constructed of Travertine 


_ Unfortunately to some modern engineers, soil has a dirty connotation. However, 
to the Geotechnical Specialist, soil is a primary material for construction and 
dirt is something swept off the floor, or which he washes from his child’s — 
‘The | final test of technology is its ingenuity in transmitting an idea into a_ 
reality. The civilization of the Western world was born in Asia and Africa 
_ and matured in Europe and the post-Columbian Americas. The initial development 


was slow, with ideas emer; in widely 
Writing fostered the long distance communication of ideas and transmission in 
of knowledge from one generation to the next. After the development of writing, = 
technology developed at a rapidly increasing rate through the my ar of 
Development in the Americas was different. 
compared to Africa and Asia, as recently as 12, 000 3 years ‘(or perhaps as long oe 
as 50,000 years ago). The settlements were dispersed and their life was dependent | a 
on hunting and food gathering. Primitive technologies arose spontaneously - 
widely separate areas, about 3,000 years ago, much later than in the Old — 


they at as. as those | in OW during 


were other in the. in Americas. 
_ The lack of a logic system and mathematical computation inhibited an orderly _ 
_ understanding of the mechanics of technology. The lack of written — 
inhibited the synergism between geographically isolated areas as well as the 
transmission of technology to succeeding generations. 
_ With these differences between the Old and New Worlds, many historians. 
have speculated on the origin of the technologies that produced the ingenious 
works that are being uncovered. 2 
Old World Contacts.—There were legends of both the Incas in Peru and ss. 
the Aztecs in Mexico that describe bearded, light-skinned gods landing from _ 
te large ships. However, the stories were embellished with such fanciful details — 
that they could be entirely imaginary. There are a few scattered finds of 
pre-Columbian European artifacts. Reo Recent sm small-craft a of the Atlantic © 


> 
| 
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demonstrate that skilled seamen, each as the ‘could have 


to the Americas, probably because of bad navigation. The chances of a safe _ 
return were small, which would explain no records of such voyages. There 
are fragmentary stories of Norse expeditions t to ‘‘Vineland’’ but conflicting 
evidence of its location. Certainly there are no records of major expeditions — 
that might have brought Old World skills to the Americas. 
7 ~The evidence of comparative technology is also weak. There are some _ 
_ similarities, e.g., weaving, burned clay pottery, adobe, and pyramids. However, 
* = i a detailed comparison shows that the similarities are superficial. The American _ 


stone pyramids are far steeper than those of the Old World d. Most were built 


box-in-box and of rubble with only a veneer of faced stone. Their ‘primary 
use was in religious ceremonies; burial was secondary. 
_ Some of the most important Old World advances are absent such as the 
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FIG. 26(a). —Plan and Cross Section, Lake Reservoir, Verde, 


Colorado 


rossings). Most important, there was a nor written 
“communication. By way of contrast, _ there are some New World ideas that _ 
: m , lime mortar and composite masonry 
using ¢ a 4 rubble c core e with only : a a thin facing of finished stone work. ao lt nai —_ 
_ The absence of Old World diseases, such as smallpox, i in the Americas before 


the indigenous people were wiped out by Old World infections after contact — 
_ The available evidence suggests that although contact was possible,  . 
was no significant transfer of technology from the Old World to the New. | 
Creatures from Space.—The idea of extraterrestrial communication is intriguing, 

but entirely lacking in evidence. The Chariots of the Gods strains at comparing 4 
the weird stylized cr costumes of the Mayan gods with the modern space suit 


7 
: 
Wheel and metallurgy (except for Western Sout merica, remote from Atlantic 
3 | 
uthwest 


and suggesting that the Nazca lines were landing fields for space craft We | ie 


can understand the ancient Hebrews explaining still older megalith construction 
(lavstones) by stating, “There were giants in the earth in those days.”’ — - 


is little excuse for such a viewpoint today (21). 

The extraterrestial hypothesis as voiced by Von Daniken is taiedaget by 

4 logic of elimination: (1) the huge structures built today are the product of 

- modern technology; (2) the ancient people possessed no such technology; (3) 

they did build some extraordinary structures; and (4) because they had no modern — 

technology, they must have received supernatural help. yin 

a This conclusion is based on on the unstated assumption that the ‘ancients could | 

not have contrived an innovative, admittedly primitive technology to produce — 
their structures. It is a form of intellectual snobbery on cannot admit the i 
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—Mummy Lake Dike, Mesa Verde, Southwest Coloredo 
‘potential for achievement without the tools and gadgets which we now possess. 
Indigenous Ingenuity. —The alternative to the hypotheses of Old World influence 
and extraterrestial intervention for the New World construction achievements 
is indigenous ingenuity. The word ‘‘engineer’’ is derived from ingenious, the 
"capacity to contrive something new. It is related to genie, the mythical being © 
_ with magic powers to create. Confronted with a problem, the engineer (the 
ingenious one) contrives a plan to transform the available materials into a solution. 
' The plan becomes a reality through the organization of labor. St 
_ The available evidence supports this hypothesis. Each of the widely ‘separated — 
- ~pre- -Columbian civilizations exhibited a logical technical evolution within its 
geographic nucleus. The problems were common among the different civilizations. q 7 
However, each evolution was somewhat different, reflecting different materials * 4 
and different waned e.g., the Maya built their temple bases of limestone Tubble- 


| 
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aa volcanic stones in a soil-vol 


___ Although there were rapid technological developments between 2, 500 and 
500 years ago, many were forgotten before the Old World conquest. Moreover, ast 8 
_ the technology of one nuclear civilization did not appear to influence that of bei 
— its near neighbors. Each invented its own, more or less independently. When — - 
a society deteriorated, its technology was largely forgotten. The lack of an 
effective means of recording and transmitting ideas was a major deficiency 


of the pre-Columbian Americas. In my opinion, this loss confirms the indigenous 


4 For Today.- —A Sait of the past is of greatest value in its message to il 


t and To" the geotechnical engineer, our pre-Columbian 
af 


their mother or grandmother, the Judeo-Christian tradition maintains that we 
were created from it, our childish instinct is to play in it, and as engineers, — 
we build with it. Thus, our work with the earth as araw material anda foundation — 


_ The arrested development of of p pre- -Columbian technology is a warning that 


‘remain within a nuclear society eventually die with it. We continue to waste 
effort in “‘reinventing the wheel.” 
‘these achievements ‘Temind us ‘the most important factor in 


he essence of | engineering. Modern 
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27.—Drain Pipe, 10 cm OD, Dainzu, Oaxaca, Mexico 
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technology has us | to 0 build bigger, faster, and even n better . However, 


_ in love with them instead of their us use. The past r remains us ; that r man nan has produced ei 
7 magnificently with little except his ingenuity and labor, £ 
- The ancient Hebrews were closer to the truth than they realized when they 
conclude, ‘‘There were Giants on the Earth in Those Days.” There were giants 
4 Karl Terzaghi was such a giant. He generated new ideas 50 years ago that 
have transformed engineering. We can express this i ingenuity if we aren’t seduced a 
by our own tools. We too can be giants, and society will be the better, - 
have contributed observations, data, and ideas from which this ‘paper 
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a EFrEct OF WETTING AND DRYING ecloya 


By Me y Mehter } Mohamed | Alam! and Asuri Sridharan? 
ied The specinuns, thos te 
_ _ The fore cover in tropical and semitropical regions undergoes large fluctuations 


both lead to degradation of the parent rock to form soil, and result in the 
aggregation of soil particles and in interparticle bonding which can impart | ans 
intrinsic: effective stress to the natural soil. The study reported herein was 
_ conceived to examine to what extent the shear strength properties of a disturbed 
natural soil are ee when it is subjected to repeated wetting and drying. 


‘The othe in large p parts of ‘the earth’s land surface are partially saturated 


_ to considerable depths and are known as desiccated soils. These soils are subjected 
- to repeated wetting and drying as a result of climatic changes. In arid areas 
— _ with deep water tables, soil humidities are controlled by atmospheric humidities — 
- todepths of 20 ft (6 m) or more (17). Changes in color and arid soluble constituents — 
and the ‘gain and loss of other constituents can occur as a result of climatic — 
changes (14). The action of water on the soil consists of solution, hydrolysis, 
and decomposition of the soil compounds (6,8). It has been reported (1) that 
alternate wetting and drying, resulting in the buildup and breakdown of soil : 
granules, is partly responsible for a continuous increase in the percentage of _ 
water r stable aggregates. The extent to which discrete ¢ or elementary, ies 


a is a function of several factors such as the grain | size e distribution, amount of * 


‘2 "Asst. Prof., Dept. of Civ. Engrg., Islamiah Inst. of Tech., Bangalore, India. A i 
? Assoc. Prof., Dept. of Civ. Engrg. , Indian Inst. of Science, Bangalore 560012, India. 
ON a Note. —Discussion open until September 1, 1981. To extend the closing date one month, 
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ASCE. Manuscript was submitted for review for possible publication on January 10, 
1980. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 
of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT4, April, 1981. 
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lay, present, nature of the saturating cation 

= and presence and solubility of excess | entities as well as organic decomposition =] 

Several ways, ys, including w wetting and "drying, have also. been reported. (13) 
by which individual soil particles can become more or less 

5 together and in effect grow into large particles. Chemical attack by suitable a 

im reagents on the majority of soil components is capable of producing covalently — A 

bonded , macrolatticed compounds possessing considerable potential strength, 

. es and under suitable conditions this can be achieved by alternate Wetting and = 
drying cycles (10). ‘Repeated wetting and drying of ‘some si soils can result 
increased strength due to aluminum oxide bonds being generated in the initial — k 

_ states, and also due to the recrystallization of ferric oxide in the later stages — 

(8). In certain cemented clays, Al, Fe, Mg, and Ca  gomponads act as cementing ~ = 

agents and influence their compressibility (11, 15). Removal of these bonds by 
chemical treatment results in loss of the apparent preconsolidation pressure. — 


Chemical Amount present, as apercentage 


Ferric Oxide (Fe,O 
Titanium Oxide (Fe, Oo, 
Oxide (K ,O) | 


Sodium Oxide (Na,0) 


4 


Magnesium Oxide (MgO) <a 
Calcium Oxide (CaO) 


on 


‘ «i i been also reported that the effect of desiccation is ‘similar to that of — 


- From the available literature it can be _ concluded that repeated wetting and 
d 


drying of a soil can generate bonding between particles which can influence © A : 


‘its shear strength properties. This may be a source of the bonding present © = te 

in desiccated soils (2, 23). This paper investigates the effect of — oe ‘s 
and drying on the shear strength behavior of 

For this a disturbed sample of des desiccated was “taken the 

yard of the Soil Mechanics Laboratory, Indian Institute of Science, Bengeiore, _ 

_ India. This soil is locally described as Red Earth and its chemical composition a 

is presented in Table 1. It can be classified as a nonlateritic soil on the basis 

of the criteria given by Martin and Doyne (4,6), at % silica to alumina ratio: 

- = 3.0. Its grain size distribution i is shown in Fig. 1, and it is composed of a 
about 0% sand, 25% silt, and 5% clay. Kaolinite is dominant clay mineral | 


i 
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TING AND DRYING 

_ eeeent, and the x-ray ‘diffractio n pat attern for this soil is pre ented it in Fig. 2 > 
& index properties are: liquid limit = 32.1%; plastic limit = 18. 1%; shrinkage Z 


“4 


i: a _ The oven dried (110° C) soil was sieved through BS Sieve No. 7 and statically 


aa compacted into cylindrical triaxial test specimens of 1-1/2 in. (38 mm) diam — 
a : 3 in. (76 mm) height in an apparatus (16) having a water content of 16% 


anda dry density of 1.79 g/cm’. The water content of the specimens “Te 
fully saturated is 17.5% and is equal to the | shrinkage limit of the soil. This — 
¥ dry density was selected so that drying would result in negligible shrinkage, 
ss The specimens thus prepared were allowed to first dry in air under controlled _ 
= “conditions, till there was no change | in weight, before being dried i in an oven 


ey constituted the first drying cycle, the sp specimens ' were cooled in air and arranged 
_ in trays lined with blotting paper and completely covered with fine sand which 
was then saturated with distilled water. The specimens, thus subjected to their 
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% FIG. 1 —Grain Size Distribution of Red Earth ee’ ; 
first rewetting cycle, absorbed moisture pons the saturated sand surrounding z 
them with little erosion taking place. After rewetting had taken place for 3 
days, ‘the spe specimens in the trays were placed in an oven and dried for 3 ays 
7 at a temperature of 110° C. After cooling, rewetting was repeated in the manner 
re 4 previously described. Successive wetting and drying cycles were performed until 
ime specimens had undergone as many as 59 rewetting cycles but had suffered — 


negligible erosion and cracking. at todd acd bas 
. Specimens were selected at random from the trays at the end of their Ist, — 


25th, 40th, SOth, and 60th drying cycle and were subjected to isotropically a 


- consolidated undrained shear tests in the conventional triaxial shear apparatus — 
= after being saturated. The saturation of the specimens was done in the triaxial 
apparatus itself in the following manner and constituted the Ist, 25th, 40th, 
50th, and 60th rewetting cycle. The € dry | Specimens were set up in the triaxia 
apparatus under ar an initial all-round pressure of 5 psi (34.5 kPa) and were allowed | 


4 imbibe deaired water from the top and bottom drainage lines for 24 h. Next, a ; 


distilled deaired water was circulated upwards through the speciinens for a 


48 h via the 
ut air in the : specimens. Measurement of water content tat 
“the end of the shear tests indicated that this procedure resulted in a degree 
of saturation as high as 99.5% in the specimens. The specimens were then 
consolidated isotropically to the desired pressures and sheared with an axial _ 
strain rate of 0.08% per min, which was Sound adequate for the equalization 
of pore pressures mobilized during shearing, == 
w Some specimens were treated with 250 ml of 0. 2N EDTA (disodium salt 
of ethylene-diamine tetra acetic acid) solution after having “undergone 50 and ~ 
60 rewetting cycles. The EDTA is a chemical that would remove carbonates, oa ; 


pp 


fatal 


 « 


gypsum, and iron compounds, that is, it would remove materials that could . — 
act as particle | cementing agents (11). The EDTA) was circulated ‘upwards through a ; 
4 


the sp 
ze the top drainage exit, was chemically analyzed. The specimens were re next 

with distilled deaired water prior to consolidation to the required pressures, 

and were sheared in the § 
_ For comparison purposes, some field saturated undisturbed triaxial test 
specimens of the same soil were extracted and _ subjected to isotropically 


‘consolidated undrained shear tones in the laboratory. Remolded specimens © were 
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WETTING AND DRYING 


WETTING AND DRYIR 
void ratio present in the saturated undisturbed specimens, and were were subjected — 


: In this study it has been assumed that no significant volume changes occur — 
in the soil specimens due to cyclic wetting: and For ‘or this reason 


CYCLES 
#so cycles 


@60 CYCLES 
t 4 SO CYCLES EDT. 
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3 —Typical Deviator Stress-Axial Strain and Pore Wetting 
and Drying Specimens | beans >: 
_ specimens were remolded at water contents near the shrinkage limit. Further, 4 
as was confirmed by inspection at the end of the wetting and drying cycles, dl 
_ cracking and erosion due to the drying and wetting process would be minimal. i j 
.- has been assumed that whatever minute volume changes were likely to occur ; 
took place a at t the | end of the first — cycle, and at a stage all the specimens — 
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_ and gain the properties of undisturbed desiccated soils, and i tens tend to be 

analogous to them in behavior. As the s specimens are tested only in the saturated 

state, the influence of suction stresses is considered to be eliminated and any 
cohesion intercept obtained is attributable only to the intrinsic effective stress 
present. This intrinsic effective stress is anticipated to be imparted by the 7, 
desiccation bonds generated. ot oh. at sub line q 
‘The effective stre stress, “originally defined as the excess of the total stress over — 
the neutral pore pressure in the ‘soil (24), has: been redefined mean. the 
intergranular stress between the contacts or near contact points of individual “ 
soil particles in the modified effective stress concept (18,20,21,22). When 


presence of electrical charges on the su surfaces of clay sized soil particles and 


(2) | 3) | 


Co in pounds per 


, in degrees 

¢ (A - — R), in 
‘a pounds | per square | 


failure strain, 


_ other bonds are | accounted for, age stress is given by (21, 22) 


= intergranular or pry stress “between ‘soil ‘particles per unit 
area of soil; o = external normal stress; 7 = effective pore pressure; A = 
total interparticle electrical attractive and R = 


me Bos is ‘defined as the conventional effective stress recognized 


(2) 


in which o” = A — R= intrinsic effective stress. 
_ While a correct evaluation of C should lead to a zero cohesion intercept — 


_ on the Mohr-Coulomb diagram, the existence of a cohesion intercept is attributable 2 


tot the p presence : of an intrinsic effective stress is — cannot be ba 
q 


4 
7 pf Rewetting Cycles = 
50 4 50(EDTA) | 6O(EDTA) 
| | 10.73 | 12.54 | 12.80] 12.87 
25.4 | 24.3 | 23.6 | 25.5 | 26.6 
5.24] 9.79] 1.14] 13.49] 13.95] 49 | 498 
| 8.9 25.6 | 25.6 25.1 25.1 
9.49 | 20.43 | 23.25 | 28.79] 28.48 | 10.46 | 10.78 : 
10.5 | 82 | 59 | 28] 27 | 42 50 
a 
\ 
a 
| 
| 4 
| | 


WETTING AND DRYING 


easily evaluated dadepundentiy. It is this intrinsic effective stress that i is expected ~ 


to be affected by subjecting the soil to repeated wetting and drying. 


_ Stress-Strain Response. —Some typical deviator stress-axial strain and pore — 
water pressure-axial strain curves obtained from CIU shear tests on specimens 


UNDISTURBED 


FIG. 4.—Typical Deviator Stress-Axial Str and Pore Sinaia Curves for the Soil 


in its Undisturbed and Remolded Conditions 


which had undergone ‘1, 25, 40, 50, and 60 cycles of ‘Tewetting are shown in — 
‘Fig. a It is seen that the shear strength [defined as (o, — o;),.., ] progressively cg 
increases, the axial strain at failure decreases (T able 2), and the soil fabric ; 
rigidity i increases as the cyclic wetting and drying process repeats. The deviator Pe 
_ stress-strain curves show pronounced peaks indicating that failure tends towards _ 
the brittle type. Fig. 4 presents some typical deviator stress-axial strain curves 

from rom cIU shear tests on saturated undisturbed and remolded 
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fr PRIL1981 
of the same soil. $ Striking aiaiitiiien with the specimens subjected to cyclic 
wetting and drying is observed | if the first rewetting cycle is considered analogous ous 
aa was found that these deviator stress-strain curves could be pind 
into rectangular hyperbolas (12,19), and the values for the initial tangent modulus 
_ (ITM) obtained from these transformed curves have been plotted against consoli- 
S&S pressure (oc .) for different numbers of rewetting cycles undergone prior 


‘ 
425 
CYCLES 


® SO CYCLES 
° 60 CYCLES © 


| 
5. —ITM- for Wetting om and Drying | Specimens 
7 is as much as 9 times over the value for the first rewetting cycle when | jie 
_ soil has undergone 60 cycles of rewetting. A similar observation can be made 
_ when the undisturbed and remolded conditions are compared (Fig. 7). Fig. 8 | 
- presents th the variation of the ITM | Fatio (defined as the ratio of ITM for any 
‘rewetting cycle to that for the first 1 rewetting cycle) with consolidation pressure 
7 - for specimens sheared after undergoing 25, 40, 50, and 60 rewetting cycles. d 
The ITM ratio is found to decrease nonlinearly with increase in consolidation 
pressure for any rewetting cycle. This decrease is dependent on the slopes 
the ITM-o, curves for any rewetting cycle and 


tO SHOGIMIE (Pigs. anu 0). © tO MICICase Wil CONSONGAUO 
& pressure and the relationship is linear. The ITM for any consolidation ae 

isa 

ic 

"| 


on the value of ITM at o. Fig. 8 also presents the variation of 
ratio for undisturbed and remolded soils. The behaviors are similar. Fig. 3 also 
presents the stress-strain response of specimens treated with EDTA prior to 
_ Shearing after having undergone 50 and 60 cycles of rewetting. A decrease 
ee in shear strength, an increase in axial strain at failure (Table 2), and a decrease - q 
in soil fabric rigidity as compared with the corresponding untreated specimens’ = 
behavior i is noticed. Fig. 6 shows the ITM-o.. curves obtained for EDTA treated 
_ Specimens. The decrease in the ITM on treatment with EDTA is seen to be 
small compared with the increase taking place on subjecting the soil to repeated _ 


— 


6201 
4 SO CYCLES EDTA | 
2 SOCYCLES” 
3 GO CYCLES” 


“= 


=A 


FIG. 6.—ITM-Consolidation Pressure Curves for Wetting and d Drying S Specimens 


_ and drying. A similar observation is made in the case of undisturbed | 


From the foregoing it can be concluded that the increase of ITM with the 
number of rewetting cycles is due more to changes in soil fabric resulting from 
wetting and drying than to the generation of bonds in the specimens. = = 
Pore Pressure Parameters.—Table 2 presents the average values of B measured — : 
for specimens subjected to 1, 25, 40, 50, and 60 rewetting cycles. The low 
; B value of 0.3 for the first rewetting cycle can be possibly attributed to the 
fact that specimen preparation was by compaction. Further, it is seen that the - : 
< average B value decreases as the number of rewetting cycles i increases. ‘The 


7 

\ 
a 
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rE 
B value of a saturated soil is dependent on the sail skeleton’s compressibility, 
j _ which is itself dependent on the shearing resistance at particle contacts, and — 
hence on the intrinsic effective stress present. The decrease in B value with 
=“ increase in the number of rewetting cycles is a clear indication that repeated ~ ; 
- wetting and drying increases the intrinsic effective stress present. This is further 
: corroborated by the data for the undisturbed desiccated soil (Table 3), which | 
_ Show a smaller B value as compared with the remolded soil in which desiccation __ 
_ From the A,-o, curves presented in Fig. 9 it is seen that although there _ 
is no definite trend in A , Values as the number of rewetting cycles increases, 
- eamgies which have been rewetted to a number of cycles resulted in smaller = 
A, values than the one with | cycle at o. values greater than about 20 psi. 
The values of with EDTA after 50 60 
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FIG. Pressure Curves for Soll in its Undisturbed and Remolded 


“cycles a; | rewetting are more than the corresponding ones without treatment. 

- ‘Fig. 10 presents the A yo, curves obtained for the undisturbed and remolded _ 

soil, and also after treating it with EDTA. While EDTA removes the bonds 
i generated by wetting and drying and thus reduces the intrinsic effective stress — 

present, it does not alter the soil fabric. It has already been seen earlier that 

fabric also contributes to the ITM of the soil. The reduced intrinsic effective — 
stress results in reduced shearing resistance at particle contacts. Therefore larger 

_ amounts of the applied deviator stress (but smaller compared to those of the 

- first rewetting cycle specimens) are supported by the pore water pressures 
_ mobilized during shearing, and hence A , values greater than those for the untreated 
50 and 60 rewetting cycle specimens on obtained. The data for the undisturbed — 

and remolded soil further corroborate this conclusion 

@ Shear Strength. —Fig. 11 presents the shear strength-o o. curves for specimens 
q subjected to different of The 
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"strength with the number of rewetting cycles clearly shows that the latter 
contributes additional intrinsic effective stress to the soil by generating some 
bonds. Fig. 12 presents | the: corresponding data for the undisturbed and remolded 

- soil. The effect of treating undisturbed soil with EDTA is also shown. The 

_ teduction in shear strength on remolding may be attributed to the res 
of desiccation bonds in the soil. Fig. 11 also shows the shear strength obtained © 

on treating 50 and 60 rewetting cycle specimens with EDTA. The reduction 
in shear strength which occurs on treatment with EDTA indicates that the bonds # 


ee by wetting and drying are chemical in nature and therefore similar 
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5-UNDISTURBED SOIL 
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AG. 8.ITM Ratio- Consolidation Pr Pressure for 

to the desiccation bonds present in soils. Since treating 50 and 
60 rewetting cycle specimens with 250 ml of 0.2N EDTA solution met 
as much as | mg to 41 mg of Fe, 10 mg to 45 mg of Ca, 37 mg to 80 mg 

3 of Al, and 45 mg to 91 mg of Mg—similar quantities were detected from the 
chemical analysis of the spent EDTA solution collected from undisturbed soil— 7 
can bec concluded that at cementing agents of Ca, Fe, Me, are 


additional intrinsic effective stress occurring in the specimens. i ae 
; It is seen from Fig. 11 that the shear strength of the treated specimens is 
ee less than that of the specimens subjected to ata one rewetting — Undisturbed 
a 


- F relationship is linear in the range of a. used, and it is seen that the shear 
| 
| 
&g 
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| 
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TABLE 3. —Shear Strength Porameters for Undisturbed and Remolded Red Earth 
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Parameters condition condition 
Con , in pounds per 
, in degrees 
Cc’. in per square 


inch 


@” (A — R), in pounds per 
a square inch 
Average failure strain, 
asa 
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FIG. 9- Consolidation Pressure Curves for Wetting and Drying Specimens 
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has a similar when treated with EDTA Fie. This ‘suggests 
that the wetting and drying process produces microcracks in the specimens, 

and suggests the presence of weaker planes in the undisturbed soils. om) it 
Fig. 13 presents the sensitivity ratio (defined as the ratio of shear strength fe 
for r any ' Tewetting cycle to that for the first rewetting cycle) versus | consolidation , 
pressure curves for specimens sheared after undergoing 25, 40, 50, and 60 rewetting © 7 

2 cycles. It is seen that the sensitivity ratio increases with the number of rewetting 
: am. since the shear strength itself increases. Also, the sensitivity ratio for 
: any number of rewetting cycles decreases nonlinearly with increase in o.. This - 
_ is because the ratio depends on the slopes. of the shear strength-o a 
on the values for shear strength ato. = =0. Fig. 13 also presents the sensitivity 


ratio (ratio of shear strength undisturbed to remolded) for r the undisturbed soil. 
Its maximum value (2. 15) is somewhat larger than for the 60 rewetting cycle 
a 


ond: 


call 


FIG. 10.—A Consolidation Pressure Curves for the =e in | its Undisturbed one 


- _ The sensitivity ratio (Fig. 13) and the ITM ratio (Fig. 8) exhibit similar trends 
; for changes in o. .. Also, if the first rewetting cycle specimens are considered 
analogous to the temolded specimens, it is ‘seen that the decrease in | sensitivity 
 gatio is of a smaller order than the decrease in the ITM ratio with increase : 
Shear | Strength Par Parameters.—Table presents the shear parameters 
in terms of total and effective stresses obtained fromCIU shear tests on ania 
« _ subjected to different numbers of rewetting cycles. It is seen that the specimens 
a which | have undergone only on one rewetting cycle Possess significant cohesion i 
va intercepts in terms of both total and effective stresses. This can be attributed — 
to the stress history effects involved in the preparation of the specimens by 4 
static compaction, and also to the shrinkage stresses occurring in the first dryi ing 


/ 
q 


4 


- i ective stresses a as the number of rewetting ng cycles i increases, es, showing an increase 
- the intrinsic effective stress present in the specimens. This further confirms | 
that repeated wetting and drying induces bonds in the specimens. The angle _ 

2 of shearing resistance displays a slight decrease, and this can possibly be attributed 
to microcracks occurring in the specimens due to repeated wetting and drying 


we place (2). Table 3 presents the shear strength data for the undisturbed 


ples 
FIG. 11.— —Shear Strength- Consolidation Pressure Curves for and Drying 
and remolded ‘ili A decrease in the value of the total stress Parameters and 
the effective cohesion is seen to occur on remolding. __ to 
_ The 50 and 60 rewetting cycle specimens with were treated with EDTA prior | 
_ to shearing yield cohesion intercepts (in terms of both total and effective stresses) 
which are less than those for specimens subjected to one rewetting cycle only, 
while their angles of shearing resistance correspond more closely to those for bo 
the 50 and 60 rewetting cycle untreated specimens. A similar effect is “0 
to occur when the undisturbed soil was treated with EDTA (Table 3). This _ 


| 
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because the treatment with EDTA is very in removing vies the 
by and , while microcracks and and other defects remain 


S - UNDISTURBED SON 


FIG. 13.—Sensitivity ‘Ratio- ‘Pressure ‘Curves fer Wetting and Drying 


Specimens and for Undisturbed Soil 

unaltered. Since varying amounts of Fe, Mg, Ca and Al ions were removed 

by EDTA, it is concluded that repeated tie and drying induces some chemical _ 
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pie 
in soil. Since the deformation behavior of the 
7 undisturbed soil and that of the specimens subjected to repeated wetting and 
drying show good resemblance, it can be concluded that the bonds present 
both cases are possibly of the same nature. 
; _ Intrinsic Effective Stress, or Net Attractive and Repulsive Pressures.—From 
the cohesion intercepts C’ obtained for the various. samples (Tables 2 and 3) 
with the conventional effective stress concept (i.e., without taking into account _ 
the attractive and repulsive pressures), it is possible to estimate the contribution 
of intrinsic effective stress (i.e., the net values of A — R) to the effective 
stress (Eq. 1) provided the following assumptions are made (22): ce ae 
The modified effective stress defined in Eq. | is valid. — te a 
The effective stress-strength relationship is linear. 
—— The angle of shearing resistance >’ remains essentially the s same with or 
without consideration of the net values of attractive pressures. 
, 4. All shearing resistance in the soil is of a frictional nature. In other words, 7 


the true effective stress is known, then the intercept should become 


Inv view of these assumptions, the intrinsic effective stress = 


specified number of rewetting cycles, and in undisturbed and remolded na 

_ have been computed and are given in Tables 2 and 3. It can be seen that 
net (A R) value increases with gumber of -Tewetting cycles and could 
as as ad psi (198.4 


_ From the results of CIU shear tests conducted on specimens prepared in * 
“the laboratory by statically compacting a disturbed natural soil (Red Earth) 
and subjecting it to cyclic wetting and drying under controlled conditions _ 

_ the laboratory, and from the observed behavior of the undisturbed soil, the 
_ following conclusions cai can be drawn: ae ps 


oti times). ‘The of the ‘soil ‘also ii increases as the wetting ‘and drying 
repeats. 

Ss 2. Decreases in the compressibility of the soil fabric occur due to increases ae 
in the intrinsic effective stress brought about by repeated wetting and ‘drying. 
‘This i is shown wn by the decreases in the B and A, parameters. Dehighiedic ¢: 

a 3. An increase in the shear strength takes. place as the intrinsic effective | 
“stress i increases due to bonds generated by wetting and drying. 
? a 4. The cohesion intercept is found to increase as the wetting 2 and drying ~ ad 
Fi tn repeats, further confirming that the process generates bonds and hence _ 
intrinsic effective stress in the sol. 

a © The maximum value of the ITM ratio (9.0) exceeds that for the sensitivity | 
' ratio (2.15). A similar trend (3.8 and 1.64) is observed when these ratios are 


compared | for the ‘Soil i in its undisturbed and states. 


|. 
4 
4. [ 
| | — . Repeated wetting and drying increases the stiffness of the soil fabric, as = . 
| 


t= 


EDTA t treated specimens show the increase ‘in ratio due 


2 wetting and drying g to be due more to changes in soil fabric than to generation "a ’ 
Treatment with EDTA decreases the shear stre1 gth s significantly, indicating 
- the bonds ; generated by wetting and drying are chemical in nature. “ae . 
a Considering the first rewetting cycle specimens as analogous to the remolded | 
a state and specimens subjected to a larger number of rewetting cycles as analogous _ a 
a the undisturbed state of the desiccated soil, remarkable similarity in the 3 
7 _ strength and deformation behavior for the undisturbed state and when the soil * 
is subjected to repeated Wetting and drying is observed. This s suggests that 


“the desiccation bonds found in can be created in the also 
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Pc relative density of natural sand deposits depend b both on textural 
properties (grain ‘size, sorting, etc.) and on environmental conditions| during 

sedimentation. The extent to which each of these factors affects porosity is 

4 of considerable interest in civil engineering, hydrology, and geology. Information — . 
on the subject is limited however, mainly because of technical difficulties in __ 
obtaining and handling undisturbed sand samples. Studies of porosity in natural 2 


: , deposits which were conducted in the past were generally bas based on samples a 
_ obtained by driving a a sampler into the unconsolidated sand, a procedure which 4 


= caused considerable ene of the sand. Accordingly, results lacked 

In the present study a technique was developed in which undisturbed samples 
are obtained by hardening the sand prior to sampling. Porosity is measured | 

on small cubes cut from the sand Samples. Some preliminary results | cleaned 


in dune sand are presented and serve to illustrate potential uses of the technique. 
hy 
Undisturbed samples were otained by impregnating the sand deposit with 
a hardening material. Two types of organic hardeners were used: A solution | 
containing 5% cellulose acetate in acetone and a polyester resin. _ i 
: The method of impregnation is as follows: A thin-walled cylinder (a bottomless — 
tin can) is pushed a few centimeters into the sand deposit, the hardening solution | 
"Sr. Lect., Faculty of Civ. Engrg., Technion, Haifa 32000, Israel. .? 
Sr. Lect., Faculty of Civ. Engrg., Technion, Haifa 32000, Israel. on a 
- Note.—Discussion open until September 1, 1981. To extend the closing date. ‘io Lilt 
a written request must be filed with the Manager of Technical and Professional Publications, 7 
ASCE. Manuscript was submitted for review for possible publication on June 24, 1980. 
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that a ‘‘bulb’’ of impregnated sand is formed below the cylinder. After hardening, 
the bulb of impregnated sand is removed and transferred to the —poery. 


The size and the shape of the hardened sample depends on the amount of 


? 


Sand Sample C nto Cubes; (b) Corresponding “Pee 


is then poured into the cylinder and allowed to filter slowly into the sand so =—s-—sF™ 
AIG. 1.—(a) Hardened 


and the cylinder’s diameter. In the present study cylinders 


15 cm diam were used and about 600 cm’ of liquid was used» per ae ; 

samples obtained were about 20 cm wide and 15 cm deep. pawn, 
3 Impregnating cellulose acetate is relatively ‘simple. Hardening occurs 

through evaporation of the acetone which proceeds so slowly that even os 


field is limited however, because of the long drying period meieiood, and because — = 
= it fails to harden when the sand is soaked or Submerged. Good ‘Tesults: with 


fas: where hardening must be fast and where the high water content t prevents 
_ drying. For such conditions the quick hardening Peet resin must be used. — 
Cubes measuring 3 cm in each direction are easily cut the hardened 
‘sample (Fig. 1). Porosity measurements are performed on each cube. The | 
procedure consists of determination of the cube’s volume by weighing in a 


liquid of known density, removal of the organic cementing 1 ener by by.heating, 
and volume measurement grains in a pycnometer: 


L ik oun nis 


m 


FIG. 2.—Location of Samples and Cubes: Relative Densities are are Indicated, as a 


grained sand is compietely impregnate elore hardening sets in. 
| 4 


The ‘technique used for | measuring the ‘volume of the cubes requires walls pores 


in the field, however, are only partly filled with hardening material, in particular : 
where cellulose acetate is used. | The acetate fills only a fraction of the pore _ 
volume and gives the samples a porous and spongy appearance. Complete sealing 
of the intergranular space is achieved by reimpregnating the cube samples with % 
polyester resin. The cubes are submerged in resin and subjected to vacuum _ 
ensuring the complete filling of the pores. One difficulty encountered when 
reimpregnating sand cubes cemented with cellulose acetate is that they soften 
during the process and therefore may be disturbed. This is prevented by applying 
a thin protective coating of polyester which is allowed to soak into the cube’s iS 
4 faces and harden before the main phase of i impregnation under vacuum. ‘et a 
a The volume of the cubes is determined as follows. The cubes are weighed a 
with an ais balance to which saps Se are  smaee by means of a thin wire. 


Weighing is done first with the sample suspended in air and then with the 
sample submerged in benzene. The density of the benzene is accurately determined 
& a calibrated pycnometric bottle, in which the temperature is kept equal to 
that in the benzene bath, in which the sample is weighed. Determination of 
q volume is by division of the weight of the displaced volume of benzene (which 
equals the difference between the dry and the submerged weight) by its density. 
re After determination of i its volume, the organic cementing materials are removed 


F of the sand grains is measured in a pycnometer filled with water. Complete 
evacuation of air from the bottle is ensured by subjecting it to vacuum. Porosity ae 


. is calculated by on basis of the volume of the grains and the ame of the 


_ Loss of mineral matter due to heati ting shou 


| 
| 
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NATURAL SANC SAND DEPOSITS 
porosity. For the determination of loss sand samples were 
first heated to a (standard) temperature of 105°C and then to a temperature _ 


of 500° C. The loss of weight was 0.3%-0.35% for the tested beach and dune i 


sands (which a are mainly ‘composed o of quartz but contain up to 10% calcite). 
Thus the porosity determinations were modified by adding 0.3% to the measured 


The reliability of this procedure of porosity determination in the field was 


a 
bane ad) 


i 4 


; 
FIG. 4 —Slip upper Part o Sample 3A 
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 iaiian in a Proctor mold. The porosity measured in cubes cut from the 
hardened sand obtained from the mold was 37.3%-37.5%, after allowing for 
q loss of 0.3% of mineral matter. The porosity as determined in the Proctor 
One advantage associated with the use ‘ti ny acetate as an hardener 
in the field, is that it allows the preparation of replica’s or “peels” of the 
5 hardened samples, revealing in detail its internal structure. Preparation of peels 
of hardened sand allows the se of the sand to be related to its internal 
‘The method of peel preparation is aes similar to that employed f fc or obtaining — 
4 peels from. deposits in ‘the field (2). Since the peels of the hardened sand are 
prepared under tse conditions, a relatively high degree of perfection — 
_ The method oat for the preparation of peels is as follows. The — 
: of hardened sand is laid on a thin layer of polyester resin covering the potion : 
of a flat dish. The resin is allowed to soak “upwards into the hardened sand. 
_ After a few minutes the slice is removed from the dish and the polyester q ~ 
- allowed to harden. Subsequently that part of the slice which has not been 
_ penetrated by the resin is removed by washing the slice in acetone (cellulose — 
acetate is soluble in acetone, but polyester is not). The depth of penetration — 
_ into the sand is sensitive to textural differences such as grain s size and packing. 


— Slight differences in texture are reflected by the appearance of minute ridges — : 


6. Deposit 8, Part of Dune, Sample 6A 


i 
/ 
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5} “and gr grooves which parallel the sand layering. The internal structure of the sand — 
Porosity Measurements in Dune ai asleoqeb Sine nag 
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area: along the Haifa Bay (Israel). It formed part t of a more extensive investigation 
_ of natural sand deposits, including additional dunes in the same area and different _ 
: of sand deposits elsewhere. The present account is primarily intended 

as an illustration of the previously described technique and its potential applica- 

Seven samples were obtained after being impregnated in the field, using a 

cellulose acetate hardener. The samples were taken from the dune’s back, slip 2 

surface, and foot (Fig. 2). Textural properties of the sampled dune sand were 

found to vary only ‘Slightly, with mean size of 0.18 mm. The sand 


‘The procedures for obtaining t the sample, performance of porosity measure- a 


mer ments on cubes cut from the samples, and preparation of peels were as described 
"previously. Since the peels and the cubes were prepared from adjacent slices __ 
of hardened sand, the distribution of porosity values could be related to oe 
internal structure of the sand. In addition to values of porosity, the corresponding 
values of relative density were als also determined. These values w were iene ll 


n in 
natural porosity as determined for the 


in which D, ‘aaave density; n = 
= porosity of 


= porosity 0 of sand in its densest state; and 
[~ _ The value n,,;,, was obtained following ASTM D1557. The sand was compacted 
_ in a Proctor mold in five equal layers with 25 blows/layer. A 10-lb hammer _ 
7 and an 18-in. drop were used. The value n,,,, was obtained by pouring the 
e sand into the mold and controlling the intensity of deposition by holding the 


: funnel a few centimeters above the surface of the sand in the mold. As a Mg 


= these tests nin WAS determined as as 37% and aS 


The internal structure of the dune, as revealed by | the peels (Figs. 3- _ 


corresponds to that of a typical Barchan dune (1). On basis of the canal 


structure, two types of sand deposits may be identified: _ ixoeob ylauoiver 


1 1 . Accretion 1 deposits: show a pronounced layering which i is horizontal or slightly 
sg. 3). They were f formed through direct deposition from a a 


inclined (see Fig 

a -driven mass of grains saltating « on a flat surface. (alae bone bosie-muldom 
_ 2. Slip deposits show a much less pronounced layering which is steeply inclined _ 
and formed by a process of and of 


q 
4 
| 
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- Accretion ‘deposits are found 1 at the back of the dune and at the flat sur ‘surface 
ey at the foot of the dune. Slip deposits are found at the inclined slip slope, 
2 at the leeward side of the dune. Fig. 4 showing sample 3A contains accretion 
_ deposits in ius upper part and slip deposits in its lower part. The sample iy 
reflects the process of dune movement showing traces of slip surfaces along 
which the progressive sand slides have taken place. Figure 5 showing sample 
6A contains slip deposits in its upper part and accretion deposits in its lower { 
"part. It reflects the encroaching of the avalanching sand over an even surface 

Accretion deposits show low porosities (36. 8%-37. 7%) and high relative 
"densities (93%- 102%). These values may be attributed 
4 force. The process is s analogous to ‘the deposition of sand in the laboratory 4 
Ps by dropping individual grains from some height. Porosities which are obtained . 

in this is way are known to approach those produced by laboratory compaction a 
Slip deposits show high Porosities (39.9%-46.3%) and low relative densities 
a Z (2%-69%). These values are explained by the mechanism of deposition which 2 
; involves avalanching of large masses of sand grains. There exists a similarity 
between this mode of deposition and that employed for the production of maximum _ ; 
, porosity in a Proctor mold. The process limits the movement of individual grains al 
which are therefore locked into an open packing. Significantly, the largest 


q 
_ porosities were encountered in the lower part of the slip slope near the foot 


7 ; of the dune, where the process of deposition is purely by mass : avalanching. 
_ Lower porosities are encountered higher up on the slip slope where the sand > 
_ is formed partly by the accumulation of sand grains rolling and jumping over 


the crest of the dune and where it is less disturbed by sliding and avalanching. — 


Karol (6) iain a method of sampling unconsolidated sediments by means 
of grouting in boreholes. The writer claims that the relative - density can be 
evaluated indirectly from the uniaxial compression strength of the hardened 
specimens, however no supporting data were provided. 
_ A method for the measurement of porosity in sand which was , impregnated =) 
by an hardening material has been described by Windisch and Soulie (12). The 
_ _ measurement i is performed by microscopic inspection of thin sections prepared : 
from the hardened sand. The thin sections are scanned systematically - along 
- transverses, using a point-counting technique. Porosity was determined statisti- 
cally employing a procedure described by Kahn (5) for the determination of is 
In our own experience the following difficulties are 
previously described method: = 
OL . The problem of obtaining a sufficiently sharp > optical definition of the grain 
outlines in a standard thin section 0.03-mm thick (which is about 10% of a 
2. The tendency of grains to ) dislodge during preparation of thin | sections. 


“3. ‘The necessity to scan a very large auatber of transverses in ome to 
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- ‘The method used in the present study for the determinstion of porosity in 

hardened sand seems to be more accurate and reliable than the previously — 
& described method, and, in addition, it is technically easy to perform. One limitation 
inh inherent to the method is, however, that it can give accurate results only in in 
ok sands with a low content of minerals | which are affected by heating. In many 
sands this limitation does not apply. CRIA 
_ The study of hardened sand is especially suited to the determination of small — 

‘scale density variations in ‘natural sand masses, a subject about which present I 

knowledge is insufficient. One limitation in the application | of the technique 


: - so mat sand bodies may be impregnated and hardened at considerable _ 
depth, i.e., by means of boreholes. We believe that some of the commercially — 
available organic — used for grouting of of boreholes could be successfully a 
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EM SETTLEMENT P 


‘By Gregory B. - Baecher' and Thomas cok 


Foundation design “requires predictions of of total and differential settlement. 
However, as soil strata are spatially variable and exploration effort limited, | 
_ these predictions cannot be made with certainty. In deterministic — 
_ in which the subsurface is characterized as a uniform or trending continuum, 
" the only procedure for dealing with these uncertainties is sensitivity analysis. — 
‘This leads to an assessment of the effect of uncertainties 0 
but ‘cannot be used to assess the effect of spatial variability on differential 
y 
settlement. A soil model based on spatially variable soil properties would broaden 
the spectrum of predictions that might be made, but for deterministic analysis a) 
oe more site information than feasibly available. Therefore, a model is 
desired which more precisely encodes uncertainties about soil properties than . 
4 either random lumped parameter models or els or sensitivity « analyses do, yet which 
— the simplistic soil mechanics assumptions of current stochastic models. a 
Given the accelerating use of numerical models and the cost experience curve 
for computation, an obvious approach is to broaden finite element a 
(FEM) to incorporate stochastic variability in soil properties. This paper reports _ 
the results of two- dimensional second-moment settlement analyses using finite 
element it methods, and compares these results with one-dimensional uncertainty 


analyses a already i in the literature. 


Previous Stupies of Settrtement Uncertainty 


Previous studies of probabilistic settlement appearing in the literature can 


be grouped in three categories: random lumped parameter models; stochastic 
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parameter models propagate parameter uncertainty through deterministic equa- 
5 tions in which soil properties are assumed uniform but uncertain. ‘Gielie 
; _ models are distinct from lumped parameter models in that they assume spatial by 
variability in the | soil profile. The one-dimensional models address variation : 
4 along vertical lines in the subsurface. Finite element models address variation - 
_ One of the earliest lumped parameter models is that of Wu and Kraft (15), oa 
which propagates uncertainty in standard penetration test data through a Terzaghi i 
and Peck correlation for the load necessary to cause one inch of settlement. 
This is a full- distribution model, assuming the standard penetration test resistance z 
and model uncertainty to be normally distributed. Ramos (10) extended this 
_ model to include uncertainty in loads, in a second-moment formulation. _ 


+  Hilldale (7) developed a one- -dimensional settlement model in which modulus ; 


elastic stress field. Similar developed by Resendiz and Herrera 
and Diaz and Vanmarcke (5). The former model assumes 


latter assumes soil properties to be autocorrelated within layers. Both are 
Stochastic finite element techniques for two- (or higher) dimensional settlements | ‘ + 
were suggested by Cornell (3) in a general discussion of the applicability of . 
- second-moment techniques to linear systems, and developed by Ditlevsen (6) 
for solution by matrix techniques. Application of stochastic finite element methods 
_ to rock and soil mechanics has been made by Su, et al. (12), Paraseau (9) 
gad Cambou (2). _Dendrou and and Houstis (4) have used a similar model for ie 
flow. Su, ‘et al. and Paraseau use Monte Carlo simulations to ) generate element 
properties from specified distributions, repeating the calculations to form a 
sampling distribution. In both studies, element properties were assumed to be 
independent and identically distributed (iid). This yields low Variance all 
_ with marked sensitivity to element size, through an ‘‘averaging’’ of random 
fluctuations over larger | and larger “sample sizes,’’ as element nt dimensions are ag 
“reduced, Cambou applied -second-moment approximations. to linear solutions 
of the finite element method which include autocorrelation among properties. 
_ The procedure for selecting element properties, and ateeat the iid iid assumption 
was made, was not discussed. 


_ Existing probabilistic settlement models, including stochastic finite ay 
methods, are extensions of deterministic techniques in which input parameters 
_ are allowed to vary. Thus, these models suffer the same limitations as deterministic 


models and should be seen only as refinements of existing techniques. Most — = 
d 


both estimation error and model uncertainty. However, it in . principle, both latter Se 

sources of uncertainty can be incorporated with existing techniques. obey teri 

_ In the present work, soil variability is modeled as one realization of a ba 

second- order stationary random field. Clearly, soil properties are not random, 

however, they are ‘Spatially variable. Since observations can _ be mate at 
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TOCHASTIC | FEM 
a finite number of locations, some assumption m must t be made to infer properties — 
at unobserved locations. Following standard practice in geological modeling, 
the assumption is made that the statistical properties of the soil are spatially 
constant (i.e., stationary), and that spatial variability is introduced by random 
fluctuations about the mean or mean trend. This is purely nypenenets: “Tae 
: uncertainty, actually due to ignorance, is divided into two artificial components, ta. 
“‘random’”’ fluctuation and estimation error, to facilitate analysis. duidw 
The random fluctuation component of uncertainty i is ¢ characterized by a point — 
variance (“‘amplitude’ ”) and autocorrelation function (‘ ‘waviness”’). Both of these 
result from the way the modeling of spatial variability is made, and change _ : 
_ if the flexibility of the mean trend is changed (Fig. 1). By increasing the flexibility a 
_ of the mean trend, both the point variance and autocorrelation distance (i.e., 


the ‘distance at which the autocorrelation Teduces to a value of of "are decreased, 


| 


FIG. 1 1.—Stochastic tic Modeling aaa Soil saan by Mean Trend Plus Random Fluctua- 


uncertainty depends on both the random fluctuations and the estimation Al 
tainty, for practical purposes there may be an optimal level of flexibility in 
Box mean trend for a fixed number of-observations. The only way the total — 


uncertainty can be reduced, however, is byi increasing the number of observations : 


- integrate deformations along a vertical line due to an exogenously determined» 
nonrandom stress field acting through a spatially variable modulus 


BY 


which p = the total Zo = formation dey depth; = the vertical a 
4 stress at depth z; and M = the modulus at 3. The total 


= 
.. 
| 
— 
— 
— 
@ pecause more Of tne Spatial Variation iS Deing explained Dy Gelerminisuc 
a _ trend. On the other hand, more parameters must be estimated for the more _ 
“mu @ _ flexible trends, and this introduces more estimation uncertainty. Since the total 
. 
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in which C[M,,,M, 


[Par = 

in which C [M,, = = the | covariance of moduli. thet two a 
b. Typically, the profile is divided into layers and summations are used in 

4 place of integrations. If the point variance of modulus is large (i.e., coefficient - 


of variation greater than about 0.2 or 0.3) a better approximation of E i) 
a can be obtained by using the first three terms of the Taylor se series 3 expansion. ~ 


i 2 terms of Taylor series expansions to calculate means, variances, and Asie 
7 ery of nodal displacements. Because stresses area result of the analyses 
= _ and depend on the — random modulus, in the present model, stress is 


Nodal calculated i in the = by the 


of equat 
in which [Xx] = the stiffness matrix; {p} = a vector of nodal Pema rte 
{P} = = the load vector. The first-order approximation to expected 


ments (i. two terms in the series expansion) is 


which 


Be ‘ 


a in which (M), = the average modulus within el element i. Ascuming the ors 
vector to be independent of modulus 


Ne 


bths z, and z,. The 
surfaceis 
4 
: 
esses, 
mce require partial 
dilferenuation Of the displac espect to element moduli. 
: | 
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The solution of ‘Eq. 9 sooults i in an a xm ‘matrix of differentials of the form: 


a(M), a(M), 

in which n = the number of displacements; and m = the number of elements. = 

With this matrix of differentials the nxn covariance matrix | 

iso obtained from the first-order a approximation 
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models, in addition to random stresses. Two- dimensional 

uncertainties depend on uncertainties in a the properties of all the elements through 
: due to common dependence « on the realizations of particular random variables. ma 
This correlation would exist even if the element properties themselves were 

mutually independent. bes ansoca 


__ Calculation of the expected displacements i is the same as in the deterministic 
case, as - long as the coefficients of variation are relatively small (less | than 


becomes non- n-negligible, which i increases es computational requirements (see, e.g., 
Ref. 6). Computation of variances and covariences is not as simple as the mean, 


and requires large matrix operations. However, treating Poisson’s ratio as 
deterministic lessens the burden by ae the terms of the global stiffness — 


in which ed = = a term in . the global stiffness. matrix. [K]; w w= w = the number — 
of elements ‘ contributing to k,,; and a, and B, = coefficients depending on 


uncertainty. As shown by. Cambou_ (2), however, ‘uncertainty in "vertical 
displacement is relatively insensitive to Poisson’s ratio. Although this nah 
is to Combos’ s model, the relative insensitivity of to 
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= n’s ratio can also be derived by p physical rea reasoning. ' Se yo 
Differential settlement, defined as an absolute rather than relative - magnitude Sey 
is strongly and cannot be approximated in second-moment 
_ Therefore, to obtain ‘moment estimates, a distributional assumption must be 


introduced. Here, Ap ~ ~ is assumed, and the moments of differential 


settlement, as a function of relative 


> =(p,-—p) (15) 


The variance of relative ‘not based on any distributional 
— 


For the stochastic finite element analysis, soil properties are specified i 


_ element means and variances of element means, and as a matrix of covariances 

2 among the element means. These statistical descriptions derive from the underly- 
ing stochastic model, and are ‘functions of element size, shape, and — 
“with respect to the correlation function. Results | presented in the next section 

_ normalized in units of the autocorrelation distance Ry. 2 


Element means are obtained by integrating (i.e., averaging) the stochastic — 


_ process over the element. For a constant spatial mean, all the element means 
are equal. Element variances are obtained using the procedure outlined by 
Vanmarcke (13), in which a double integration of the autocovariance function 

_ is performed over the ‘element area. The autocovariance is simply the autocorrela- 
tion multiplied by the point variance. In the present case, mega integrations 
were performed numerically using a Monte Carlo algorithm (e.g., 8, 14). For 
simple geometries and autocorrelation functions, Vanmarcke presents tables of | 

_ variance reduction functions which may be used in evaluating element variances. 
Element covariances were obtained by double integration of the autocovariance 
function over the two elements in question 


in which C[M,,M,] = the autocovariance function of modulus between (x,,z,)_ 


> and (x,,z,). The dimensions x,,z,and x,,z, = the horizontal and vertical dimen- 
sions of the two elements, appropriate to their respective pennetes. The values 
_A,and A, = the respective element areas. of 


| a 
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is ‘to note that variances and covariances differ from 
element to element, depending on their geometry | y and area. Previous studies a 
which report strong sensitivity of solutions to element size have in some cases 

x used constant statistical properties over differing element geometries and areas; s 
_ this in itself could explain the sensitivity. With constant element variance, as =. 
the number umber of elements i increases, the ers in nodal ain penn a 

- Computing element covariances for grids with many differing element shapes, 
sizes, and orientations is costly. Approximations, as used in the present study, 3 
are discussed in Appendix]. 


Moves 
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The finite element 74 used in the present study is coarse to reduce computa- 


4. 
2 —Finite Element Mesh for Two- Model 


model and exact deste: solutions i is good (Fig. 3). By elastic inten i stress 
_ on boundaries is less than 15% | of the applied stress. +The footing is J 
& to be flexible, the base is a 1 fixed boundary, and lateral boundaries provide x 
horizontal restraint but allow vertical displacement. To facilitate comparison, 
the one-dimensional model used the same base geometry. __ re a | 
_ Two cases were considered: a homogeneous profile with spatially constant a 
mean, point variance, and autocovariance function; and a nonhomogeneous profile — 
with modulus increasing as the square root of depth, _ constant coefficient of 


--variation and constant autocorrelation function (i.e., variance and covariance 
increase with depth). Two autocorrelation functions were used: ts onl Py, 
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isotropic and anisotropic correlation were considered. 
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IG. 4- 4.—Total Settlement Uncertainty for One- and Two- Dimensional ‘Models 

Results for total settlement uncertainty for one- and cwediimenataial analysis {Jf 
for the simplest case of homogeneity, isotropy, and exponential autocovariance — 


are shown in Fig. 4. Results are normalized by the coefficient of variation — 
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of modulus and the autocorrelation distance The fo following trends 


For R,/B 0, COV [e] - 0, for Ry 00, COV [p] - 
[M]}. These correspond to perfectly uncorrelated at perfectly correlated values 
of modulus, and agree with intuition, 
_ 2. For other values of R,/B the one-dimensional model consistently yields 
larger total settlement uncertainty than the two-dimensional cn 
3 3. COVs for different points on the footing differ by only about 10%. 
Correlations among nodal displacements are shown in Fig. 5, and the implied > 
differential settlement i in Fig. 6. While the one-dimensional model yields consis- 


TWO-DIMENSIONAL MODEL 
ISOTROPIC CORRELATION 
= = ANISOTROPIC CORRELATION 


rams 


7 interesting result of a critical value of R, / B which maximizes expected dif! ferential 
settlement. In retrospect, this is intuitive. The “‘worst”’ autocorrelation distance — 
te differential settlement would seem to be i in the range of 0.75B-1.0B. Given | 
typical autocorrelation distances in the range of 100 ft-200 ft (30 m-60 m), 
this would imply large differential settlements for mat widths « of about the same 
D. Introducing anisotropic correlation produces the results also shown in n Figs. 
4-6. Here the horizontal autocorrelation distance R, was expanded to ten times _ 
the vertical, R,. Numerically, this was accomplished by defining a new distance = 
D’, obtained by the horizontal dimension ten | times in 


— | 
“4 FIG. 5.—Correlation of Total Settlement Among Footing Nodes 
tently higher predictions of differential settlkement—because functional depen 
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FIG. 7.—Total Settlement for One Two- Dimensional Models. for 
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FIG. 8.—Correlation of Total Settlement Among Footing Nodes for Nonhomogeneous 
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with the anisotropic model is more a matter of the than 
a real difference. When plotted as root mean square R,, for example, the 
differences are reduced. The interesting aspect of Fig. 4, however, is that one— 
would expect vertical autocorrelation entirely to control the variance of total 
settlement, and this would seem not to be the case. For the same vertical 
| Ro variance in total settlement is sensitive to horizontal correlation, at least 
in the two-dimensional case. This is not true in the one-dimensional — a 
which does not include mechanistic interdependencies. Differential settlement, _ 
which one would assume to be controlled by horizontal autocorrelation, seems 
to be so. The influence of anisotropic autocorrelation on differential settlement _ 
_ is more striking in the one-dimensional model than in the t two- dimensional model. ae 
a The influence of autocorrelation function : shape, at least between the ¢ exponen- 7 


é conclusions are similar to the homogeneous case with one important — 


_ The nonhomogeneous case yields larger variances in settlement, and smaller a 


correlations among settlements immediately under the foundation. This yields | 
much larger differential settlements. This is again caused by the functional 
_ dependence. Covariance depends more on soil elements near the top of a 
stratum where the differentials ap/a(M) are the greatest. The yo meng oe 


data from different locations in minimizing errors of predictions, whether or 
i not numerical modeling is used. It provides a way of accounting for _ the © 


2 
Conciusions 
Se 3 
econd-order stochastic finite eames techniques allow predictions to be made 
: in random field problems based on more realistic analyses than current one- 
_ dimensional techniques provide, and allow the incorporation of spatial variability 
of material properties which deterministic finite-element techniques do not. In 
‘principle, stochastic finite element methods can be applied to a broad spectrum — 
of problems, at computational costs which are reasonable. 
_ Comparison of one-dimensional and two-dimensional settlement analyses — 
indicates significant differences in predictions. These differences are primarily — 
due to randomness in the finite element method stress field and mechanistic 
dependence on commonly shared random variables. Disparities are particularly a 
‘distinct i in predictions of differential settlement. wt. 
4 Care must be taken in selecting statistical element properties to sudennate . 
reflect the underlying stochastic model. If care is not taken, the resulting 


reflect the element discretization and not the underlying 


| 
i 
even the correspondingly low-element variance does not entirely cancel. This 
7 _ isnoteworthy, for it provides a vehicle for deciding the relative value of exploration 7 i 
a 
| _ This work was funded under PROJECT SEAGRANT, through the Massachu- J 
Institute of Technology Office of SEAGRANT Research. 


€ as described earlier, the large number of element pairs (C3,) makes this 
a costly procedure. Therefore approximations were introduced. Eq. 19 is an 
. averaging of covariances among pairs of points, one from each element. Thus, 

for t widely separated elements (i. €., with respect to | the autocorrelation distance) . 

the autocorrelation function is approximately linear in the range e of interest, 

and the covariance at element centroids can be used in an approximation. — 

However, for elements that are not widely separated, the approximate linearity _ 
_ of the autocorrelation function is not a good assumption, and the use of centroid — 


covariance in an approximation i is rae into | question. For closely spaced 

I ] 190 


REDUCTION FACTOR 


o——o y*=Texpi-€/R,) 
SIMULATION (EQ.!9) 
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elements two approximations were considered, which seem to bound the elt 


of of element covariance as as calculated from Eq. 19: 


in which € = the vector separation distance between element centroids. The 
variance reduction | factor is defined as V[(M)]/V and ‘is introduced 
ps 
as a covariance reduction factor. Results of these approximations are shown 
in Fig. 10 for two possible geometries. For geometric reasons Eq. 22 provides. 
a better approximation than Eq. 23 when R(-) is approximately linear over 
q ‘ the range of separation distances between pairs of points, one in each of the 


two elements. The reverse seems to be the case when R(-) is not approximately — 
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~The following symbols ar ad used i in this pe paper: 


area of one i; 


B = 


2 


distance argument for autocorrelation function; 


expected (mean) parameter value within element; 
influence coefficient for total settlement; 
influence coefficient for differential settlement; 
element of stiffness matrix; 


autocorrelation distance [R(R,) = 1 


nud 


= number of elements contributing to ke hy 


vector distance between. element ko 
total vertical settlement, at node i i; 
nodal displacement vector; el 
differential settlement; 
telative settlement; 
vertical stress at depth z; 
= variance parameter; and at 
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organic matter includes plant and animal residues in various stages 


_ processes. Problems associated with decomposition in ouniie soils are often 
avoided | because of uncertainties in organic- -content Measurement 


| 
. can n be ignored. A decrease i in n strength o or + additional settlement during the s service 
ite of these materials may go unnoticed. In many cases the potential for problems — 
* is solved by organic-soil removal and replacement with mineral soils. A need 


degree of decomposition, and prediction of future engineering properties. wyatt 


E for improved methods suitable for evaluation of the organic content, 


; Measurement of soil organic content requires that mineral solids be separated 
y - from organic solids. Destruction of the organic material by ignition is the method 
, : more commonly used by engineers. The high temperatures needed for the ignition 
: if test require that a correction be made for loss of surface hydration water from 
the mineral solids. Also, the ignition test provides only partial information o 
the degree of decomposition based on weight of organic solids. Both 
ndecomposed organic matter and | microorganisms f. ormed during decomposition 
are destroyed by the high ignition temperatures. Information on the degree 
of decomposition and possible future decrease in = organic content would 
be useful on many engineering projects.§ 
Prof. of Civ. Engrg., Univ. of Detroit, Detroit, Mich. 48221. 
*Prof. of Civ. Engrg., Michigan State Univ., East Lansing, Mich. 48824. Me ht 
: a Note.—Discussion open until September 1, 1981. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, , 
ASCE. Manuscript was submitted for review for possible publication on April 22, 1980. 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings _ 
of the American Society of Civil Engineers, @ASCE, Vol. 107, _ No. _GT4, April, 1961. F 


3) 
&g 
| 
a 


4 


_ Nutrients added to the model soils, in proportions similar to those found ‘in 2 
an average bacterial cell (13), and storage of anaerobic samples at a temper: ure 

° i Wi 

& close to 35° C helped accelerate decomposition in the laboratory. Measurement 
of organic content of decomposing soils by use of the ignition test is reviewed — 

: and a correction factor which appears suitable for use with natural — 
Soils is presented with an error analysis. Al relationship between 1 microorganism — 
cell content, undecomposed organic matter, the mineral fraction, and the — 


on soil organic content measurement ‘using the ignition test. OMT = 


a EXPERIMENTAL Wom 


‘The experimental work it involved preparation of model organic soil samples, 


; materials, sample preparation, decomposition photos, and ignition test are 

a Materials Studied.—The kaolinite included particles with 100% passing 40 pm, 
93% passing 10 pm, and 42% passing 1 pm. The grain size distribution was 

relatively uniform. Consistency limits included a liquid: limit of 47.8%, plastic 

limit of 27. 5%, and plasticity index of 20.3%. 


_ Pulp fibers included a range of sizes with a weighted if ibis length of 1.6 


FIG. Appearance of Pulp Fi 
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2. —Partially Decomposed Pulp Fibers, 1000 Magnification (1): (1): (a) ) Cracks racks at 
Fiber | Surface Indicate Fiber Breakdown; (b) Advanced Stage of Decomposition 4 


_ 


— m?/g ‘of ‘ee fiber. A value of 1.54 was used for the specific gravity of the _ 
almost pure cellulose fibers. 
a Sample Preparation.—Freeze-dried pulp fiberboard was separated into a iefty 
mass. Selected proportions of dry kaolinite and dry fiber were mixed. Next, 
distilled water was | added in amounts needed to form a slurry. Nutrients and : 
‘seed microorganisms, in predetermined proportions, were added directly to the 
clay-fiber slurry to help accelerate anaerobic decomposition of selected laboratory : 
: samples. Nitrogen served as the base by which other nutrient quantities were ‘ 
bed computed. Approximate ratios giving optimal decomposition r rates included C/N 
- 1, P/N = 1:5, and (Mg, Fe, K, Ca, and 'Na)/N “8: 16. ‘Pulp fiber supplied | 
7 - carbon, NH, Cl supplied nitrogen, K, HPO, for phosphorous and potassium, 
MgSO, for magnesium, CaCl, for ‘ealiam, and FeCl, for iron. About 1% 
(dry-weight basis) of municipal sludge provided seed microorganisms. —__ 
a Decomposition Photos.—The extent of microbial attack on pulp fibers was 
using a scanning electron microscope. Details for sample preparation 
and mounting of samples on small aluminum stubs are given by Al-Khafaji 
§ A pulp fiber with no decomposition, shown in Fig. 1, has holes with dimensions 
roy from about 2 »m up to 15 ym. An advanced stage of ——— 
‘al in Fig. 2, indicates a general breakdown of the fiber’s structural integrity. 
The ability of microorganisms to maintain high concentrations in small cavities : 
may accelerated the ‘Process at certain locations as shown at the fiber 


- surfa ace in n Fig. %b). These photographs help one visualize changes i in the oll 

organic material due to decomposition. 
Ignition Test.—Soil organic material is wherens te mineral 1 

i constituent is incombustible and ash forming. Separation of the two soil compo- 


nents was accomplished by firing an oven- -dried (105° C) sample in a 

furnace at specified temperatures until the organic material was reduced to 

os ash. Test equipment included a platinum crucible with lid, an analytical 

balance having a sensitivity of 0.1 mg, and an electric muffle furnace capable | 

of maintaining selected temperatures within +25° C. Weight-reduction curves 

4 & the materials, pulp fiber and kaolinite, were determined using temperatures 

from 200°C to 900°C with oven-dried fiber sample ‘sizes sufficient 
to yield minimum ash (or residue) weights of 10 mg-20 mg. For samples containing 
mineral fractions, oven-dry weights up to 10 g were used. After ignition and 

Pe some cooling, the crucible with sample and lid was placed in a desiccator. 
When cooled to room temperature the weight was a 20 i The ash (or 


_ Weight of ash or 
The organic fraction X,, without corrections, was” taken to (l 
correction for loss of surface hydration water from the mineral fraction is 


mm, based on the test procedure given in Ref. 5, and typical diameters of 
about 20 were observed in the electron microscope. Surface area 
: 
| 
eveloped in later section. Urganic Tractions for the decomposing model sous 
were obtained in the same manner using small oven-dry — : 


ORGANIC. CONTENT MEASUREMENT Pe 


= Soil organic material fechades plant and animal residues in various stages: ay 
of decomposition together with soil microorganisms that are involved in these 
processes. . Background information, helpful i in understanding soil organic-content 
measurement, is used in developing a correction factor for loss of surface 
hydration water from the clay minerals in multicomponent and natural soils. __ 
= error analysis helps evaluate the accuracy of soil organic-content measure- ae 
ori adi, moti ine Yo 220! cee 
initially oven dried by 
at 105°C. for = hour 


Mantmortt 


Temperature (cy) 


FIG. ‘4.—Weight Reduction for Pulp Fiber with Ignition at at ‘Wemapivabied of 100°C- 
_ Background Information.—For engineering purposes the ignition test is the 
most common laboratory method used for measurement of soil organic content. 


The ignition test, described previously, can err by up to (15% due to loss of 
surface hydration water from the clay minerals (11). The American Society “sit 


factor. Pont with small mineral fractions would lose little hydration 


. { FIG. 3.—Weight Re on for 12 h) for Pulp 
Fiber, Kaolinite, and Montmoriionite (1) — 
i _ _______ HH 
' 5 r for Testing Materials (ASTM) standard method (Designation D 2974-71) of test 
for orcanic matter 20+ moaterioalc inuoaluec with no 
it 


ont pon indicate that matter nee dried at 
My 105° C) loses about 3.8% by weight when ignited at 550° C. Arman (2) r recommends — 
“use of oven-dried samples followed by combustion at 440° C with no correction — 
4 factor. Franklin et al. (6) found that burning at 400° C-450° C until constant — 
weight was attained gave reasonable organic contents. 
_ To provide accurate information on change in organic content for the de- 
composing model organic soils, weight-reduction curves for the model soil 
components, pulp fiber and kaolinite, were determined as shown in Fig. * 
These curves show that at a temperature of 400° C maintained for 12 h there a 

was a minimum loss of surface hydration water from the kaolinite. All the pulp 
fiber was ignited at 350°C. At higher temperatures use of the ignition test 
_ without a correction factor would give larger errors. Weight reduction v versas 
time curves are given for pulp | fiber in in Fig. 4 for kaolinite in Fig. 


J 


gs 


5.—Weight Reduction for Kaolinite with Ignition at Temperatures ot 900° 


two model si soil with oven-dry “weight can tee 


4 


_ Choosing an appropriate temperature T and duration D for ignition (Figs. 4 = 
and 5) tt the ‘sample weight W, after ignition c can be expressed in the form 
in which C, A, D) and C,(T, D) = the weight fractions remaining after r ignition 
for the pulp fiber and kaolinite, respectively. Eqs. 2 and 


7 ) 


W, 


define the organic fraction X,= = W, | W, and the mineral 


san 


| 

| 

>: | 


— 


"Selection of T 400°C D = h gives C, D) = 0.00 from Fig. 


at and C (T, D)= - 0.986 from Fig. 5. Substitution into Eq. 6 and rearrangement _ 


ay for the model organic soils. 


q soil. For special cases where the mineral type is known, this ates provides 
, 3" a tool for measurement of soil organic content without the need for complete — 
combustion of the organic matter. colle are re considered | in the 
‘next section. 


Combination Pulp fiber | Kaolinite | Sand 
Cy = 00 
= 0.10 
= 0.00 
ty 


Multicomponent Soil.— An extension of the two-component system is possible 
for soils with n materials with different thermal and or. characteristics. — 
Consider the preceeding analysis in a generalized fo form, thus — 

=D CMs j= 1,2, 


which i = the ith a given ‘soil; 


temperature T and duration D; C, = the weight fraction remaining of the ith 
soil component; and W, = the total weight of the soil sample at the jth combination 


of T and D. Fora three ae soil a = om summing Eq. 9 over i gives y 


‘Consider a three-component model soil wits contains fiber, kaolinite, 


: and sand. Select three combinations of T and D; (T,,D,) = (105°C, 24h); _ 
(T,,D,) = (500° ay 1. 5 h), and D,) = (900°C, 1. 5 oem 


\ 


= 
4 
A 
4 
q bmponent Model Soil 
ple weight for 
& * 
a | 
| 
| he jth Combination 
— > 
: _ example it was assumed that sand loses no surface hydration water. Using . 


1.316 W, + 13.158 W, 


W,= 


= 0.864 8686 W + 8.816 W, 
Note that only the oven-dry weight WwW, and the ignition weights W, and Ws 


needed to determine the weight fractions of pulp fiber xX x,= W, / kaolinite 
* Natural Organic Soils.—In the usual case mineral type(s) present in natural tof 
¥ 
_ organic soils are unknown and the methods described previously will be of 
E limited value. The immediate question concerned possible use of these methods 
_ for natural organic soils. Further development was encouraged by the following 


os 1. Dehydration weight losses for most minerals at temperatures around 400° c ; 
fall betw — that piaeageii by quartz and bentonite (dominant mineral compo- 
ah 


pied 


° 


ok 


is as in Fig. 6. give an upper 

_ 2, Composition of the soil organic matter is very complex and varies with» 

- and origin, however, the main organic constituent is generally assumed { 
- to be similar to cellulose (6). With this Lemaeaeed the Spectrum of of organic 


inverted matrix the weight of pulp fiber 
| 
= 
7 “ ignition must be selected so that a minimal loss of surface hydration water | d 
| 


occurs. Ball (3) indicated that critical dehydration for occurs between 
_ 4, Complete combustion of pulp fiber occurs at 400°C as shown in Fig. _ 
3, This fiber is similar in composition to many soil organic | materials is which — 7 
are fully combustible ile at temperatures below 400° ae (2,8). ods 
45 Examination of Figs. 3 and 6 and recalling the four items previously listed, 
wt attention is directed to two facts. Regardless of the mineral type the weight 
_ loss due to mineral dehydration appears to be limited to about 2% at a temperature F ; 
of 400° C provided that the soil samples were initially oven dried at 105°C. 
_ Skempton and Petley (14) observed that thermal dehydration curves for a number 
is 7 of common Soil minerals (excluding montmorillonite) showed less than 3% — 
; 3 loss due to dehydration at 400° =. Depending on sample size ignition at 400° C 


“without regard for mineral content or origin. ‘Consider Eq. 6inthe form 


pick C, (400° C, - ¢, ¢ 12h) 


in which Cc; (400° C Cc 12 h) = : 0 0 for ‘the organic matter and C, C,(400° C, RD bh) 


= 0.98 for montmorillonite (Fig. 3). Substitution of these values into Eq. 12 
; gives the organic fraction X, of a soil in terms (of the initial oven sad weight 


and sample wee W, at 400° C after 12h ignition, thus 


_ The mineral fraction is given by Eq. 7. The correction of 0.02 is smaller than : 
b é d the value 0.04 given by Skempton and Petley (14). Two possible reasons include 
- their higher ignition temperature (550° C versus 400° C) and their assumption 
that the organic carbon = 0.58 x organic content. Broadbent (4) stated that 
the constant (0.58) is based on very old work and the value may range from 
0. 5-0.58. Greater accuracy would be attained if dehydration information on_ 
- the mineral matter contained in a given natural organic soil were used in Eq. 
12. A discussion on possible errors is given i in the next section. — Ay 
Error Analysis. — Advantages of the suggested procedure and correction factor 
- become more apparent when errors associated with other ignition methods are 
_ evaluated. Three procedures considered include: (1) Samples air dried at 25° C 
- before ignition at 900° C for 75 min; (2) samples oven dried at 105° C before 
ignition at 900° C for 75 min; and (3) samples oven dried at 105°C before — 
ignition at 400° C for 12 h. Note that complete combustion of the. organic matter 
7 occurs in all three procedures. The total weight loss fraction X, experienced 


by the ignition soil samples can be written as 
in which C,, = the weight-loss fraction due to mineral dehydration; X,, 
mineral fraction; and = the true organic fraction. Note that C m iS 
temperature, duration of and type. Recall that (X,, + X; ) 


: 
| 7 
| 
q 
. 
— 
2) 
| 
j 


thus the fraction = (x, -Cc — C,)- Note 
the experimental organic fraction X, will always be gi greater than the true organic _ 
fraction X; since it includes loss “of mineral surface hydration water. As a 
result the weight loss fraction due to ignition, X, , with no corrections, is aetuiied 


‘mental and true “organic fractions. 


rer 4 


‘or MPD =|c,————_| 


Examination of Eq. 15 shows thet th the MPD i isa function of the loss 
fraction due to ignition. Note that the highest MPD that can be attained equals 


~ 


Oven dried @ 
05 |./ Ignition @ 400°C for 


FIG. 7.—Maximum Possible Difference (MPD) Between Experimental and True Organic 
Fractions for Three Ignition Test Procedures (1) 


C,, and this occurs only when C,, _ = X,. This means that the maximum hl 
in the ignition test will occur only when the weight-loss fraction equals that 2 
Consider first a soil sample that has been oven dried at 105°C followed | - 
: by ignition at 900° C for 75 min. The maximum anticipated error in the organic 
fraction will be about 0.15 (Ref. 11). This represents the maximum possible 
weight loss due to surface hydration water for any mineral. If weight loss fractions 
are less than 0.15(X, < 0. = this loss must be attributed to mineral dehydration =) 
_ and substitution of C,, = = X, into 0 Eq. 1S gives 


7 


: 
: Now determine the maximum possible diiference (MEI) ween the €xperi- ; 
&g 
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= 


xX tal lo | 


the weight- -loss fraction due to ‘ignition, is ‘greater than 0. it 


be that C,, 15 and Eq. 15 becomes 
MPD = | = (6. 176(1 — X,)| for = X,= 1.00 ba (17) 


These r relationships are » shown i in Fig. 7. Note that the MPD will be zero w when — _ 
_ X, = 1.00. In this case the mineral fraction must be zero. A similar = | 4 


E 4 was made for an air-dried sample and the new procedure as shown in aw 
_ 7. For the new procedure the 2% error is eliminated using the linear equation — 
solution a1 and the only error remaining is due to the | presence of sand | which | 
experiences about 1% rather than 2% weight loss in surface hydration water. 
Thus the maximum error of +1% may be expected for Eq. 13 using oven- dried 
‘samples (105° C) followed by ignition at 400° C for 12 h or until constant weight > 4 


new bacterial cells, and a decrease in the organic solids content. ” Completion 
of the process means that some organic material will remain, the non- gr: el 


By-products 


ra 
Microbial 
cells 
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organic 
matter 


FIG. 8. —Schematic Diagram Showing PPA So before and after or Partial 


his 2 
‘The: ignition test of soil content does 
By distinguish between types of organic material, hence a correction factor 


is needed when computing the degree of decomposition based on the initial 


___ Background Information. —Both aerobic- and anaerobic- decomposition proc- 
esses result in ‘substantial weight losses of organic matter. It is reasonable 


— 


“i ute these weight losses as a quantitative m measure of the degree of decomposition, 


— 


must ff 
| 
] 
q 
| 


‘in n which W,, = the initial weight of organic matter in the soil and W, 


a of endecompened organic matter after partial decomposition = shown 
in Fig. 8. However simple Eq. 18 may appear, "separation of the organic material ‘ 
from the —— fraction is very difficult. Define the initial organic fraction — 
X,. = W,,/ (Wyo + W,,) and the organic fraction after partial decomposition 
W,, (Wi, + W,,) where symbols are as described in Fig. 8. Solve _ 


4 


substitute into Eq. | . 18, and with rearrangement write « oF om 


of 


4 


esses involve conversion of organic solids into stable humuslike products. . The > 
substrate (material on which the organisms live) i is transformed into water, gases, a 
volatile acids, and new microorganism cells (13). A schematic diagram showing 

the changes for a typical organic soil undergoing decomposition is shown in 
Fig. 8. The microorganisms are lumped with undecomposed organic matter in "f 
the usual ignition test for determination of organic content. Consequently, the 
= of decompusition will always” be greater than the value given by Eq. <i 


19. The error would be corrected if the microorganism oe. content were - 
and subtracted from the current organic content. 


- defined the cell yield, y, of a decomposed substrate as the ratio of the Ceteees 
to the mass of the decomposed substrate. For aerobic conditions the cell yield a 
- ranges from 0.15 g/g-0.4 g/g of decomposed substrate. A range of 0.05 g/g-0.20 

_ g/g of decomposed substrate is predicted for samples undergoing anaerobic > 


These a starting in the absence of more 


in which initial organic matter in the soil; Wye 

- weight of organic matter converted into microbial cells; W,, = the weight of — 
| Sasegempeces organic matter; and W,, = the weight of decomposition byproducts. 
4 -™,) of 1. Sebe defined by Eq. 18 can be rearranged to the form 


= X aW, into Be. 20 20 and rearrangement gives 


| 
a 
an 
| 
The cell yield. y. in a given substrate can be expressed as 
(l-y 


in which. y= = cell yield; X,, fo and X,, = = the ‘initial and final organic fractions 
given by Eq. 13. The factor, (1 — y)~', is a biomass error term and it is 
dependent on environmental and microbial nutritional parameters related to the 


decomposition process. ‘When estimating the degree of decomposition ~ an 
organic soil deposit one can use an approximate 


the anaerobic case. The alternative solution would be to: assume ne y to 
e in Eq. 22. Values for X,, and X,, would be based on ignition test results. red 


Model org anic soi ils p repared from kaolinite. and pulp fibers have been | studied 

relative to use of the ignition test for soil organic content measurement and 

a evaluation of the degree of decomposition. These conclusions are intended to — 
reflect the findings of this investigation and are limited to the soil materials 


__ used and to the test procedures employed. hoe. 


oF _ 1. Use of dehydration curves for kaolinite and pulp fiber made possible 2 the 
fp elimination of major errors associated with the i ignition test for soil organic-content _ 
determination. Linear equations and coefficients permitted 
Cr organic-content measurements for various combinations of temperature and 
Ly 2. The degree of organic- solids decomposition, X,,, was defined in terms 
of the initial weight of organic matter, W,,, and the weight or undecomposed % 
organic matter after partial decomposition, W,,, in the form X = (W,, — 
The cell yield, y, takes into account the fact that microorganisms 
formed during decomposition are also destroyed in the ignition test. = 
eh. 3. Ignition test data taken during different stages of decomposition may be 
used for computation of | the degree” of decomposition. This approach can be — 
used to describe changes in the organic-solids content of — soll deposits — 
=. 4. Visual evidence of changes in a pulp fiber due to microbial activity were 
shown using a scanning electron microscope. Fiber breakdown included a 
reduction in length, diameter, and rougher surface features. The photographs 


_ illustrated some be the factors contributing to changes in soil properties d — 
— 
a This study was supported by National Science Foundation Grant No. ENG75- — 
13765 and by the Division « of Engineering Research, Michigan State University. 
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weight fraction remaining a after ignition; 
= weight- loss fraction due to mineral dehydration; 
possible difference between experimental and true organic 

- decomposed weight organic matter; 


weight byproducts (gases, H, volatile | 
w ight mineral material; 


a eight of soil component; wie 
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"CENTRIFUGE MODELING OF 
to « EMBANKMENTS" 


2 
By Mosaid M. Al- Hussaini," M. ASCE, Deborah J. Goodings, ndrew N OE 


Current practice for the design and analysis of coal refuse embankments 


7 facilities as outlined by the Department of Interior (2) involves the use of various 

analytical solutions. . Soil modeling utilizing centrifugal techniques to increase 

4 gravity forces offers an attractive method for verifying the available Sasa 

_ procedures. In this context a series of model embankments with various geometries _ 

i) and materials properties were constructed of coal waste, accelerated using a 

a 4-m radius centrifuge to 100-120 gravities, and subjected to various reservoir 

levels: until failure occurred. Companion triaxial tests on specimens at similar 

model densities were performed to obtain shear strength and permeability 
for comparing results of seepage and 


with those observed in the tested models. b= 


Centriruce Mopeune Principtes Pertinent TO SEEPAGE AND Srasury 


The basic concept of centrifuge model testing is to create a scale model 
similar in material properties and boundary conditions to ‘prototype, and 
to subject the scale model to an acceleration such that the increase in self-weight 
stresses at any point matches those at corresponding — in the petetyps. a 


‘= Prof. , Kuwait Univ.; formerly Research Civ. Engr., - United States Army Engr. 
Engr., TAMS, New York, N.Y. | 
Prof., Cambridge Univ., England. +. 
Prof, Univ. of Florida, Gainesville, Fla. 32611; formerly Research Civ. Engr., ‘United 
_ States Army Engr. Waterways Experiment Station, Vicksburg, Miss. 
-Note.—Discussion open until September 1, 1981. To extend the closing date one month, * ; 
a written request must be filed with the Manager of Technical and Proiessional Publications, _ 
ASCE. Manuscript was submitted for review for possible publication on October 22, 
1979. This paper is part of the Journal of the Geotechnical Engineering Division, ae 
5 of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT4, April, 1981. a 
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In ‘the case of modeling eubachnden stability, the pertinent ‘controlling factors 
are: (1) Height of embankment and strength of material; (2) seepage conditions; 7 
and (3) pore pressure dissipation. Assuming that steady-state seepage through a 
an embankment can be represented by a flow ‘net of N, flow channels and = 
Ny, -equipotential drops, then the rate of flow, rs at ‘earth’s gravity 

an unit width of the embankment with permeability, | c. is iis by: 


in which Ah = the head we ‘the “For a model of identical 
_ soil at 1/N scale in earth’s gravity, the rate of seepage flow/unit width through 
“the model would be 1/N times that of the prototype, as the pressure difference 
associated with Ah in the model and in the prototype is in the ratio of 1/N, 
€.g-, @, = *. x N. However, if this model were to be accelerated to N gs 
(in which g = the gravity acceleration), then the velocity of flow of ' water 
_ out of the model would increase by a factor of N due to the directly proportional “ 
increase in pressure head (assuming k is not affected by the acceleration). —_ ‘Ss 
the resultant relationships for a model accelerated to Ng are: = Ae 
_ For the centrifuge models in this investigation, the model width was 6 i : 
(154. 2 mm), so the measured rate of seepage flow, Q, in cubic feet per inde a 
est a model at 120 g would simply be multiplied by 240 (120 x 2) to 
give an estimate of the equivalent flow per foot width of the ae 
Considering failure of undrained clay slopes, a condition in which 
the nondimensional group, c/y,H, referred to as stability number appears in ; 
"stability calculations. In this case, c = the undrained cohesion of the clay; — 
| , = the unit-weight; and H = the height of the slope. Such a clay slope — 
on be modeled using identical clay with identical and and a slope 
_ height reduced to H/N, in which N = the scale of the model. By increasing 
forces N times, the stability number will be c/Ny,(H/N) = ¢/y,H 
as for the prototype. For example, a prototype clay slope 75-ft (22. 875- -m) high 
- ome correspond to models at 90 or 120 gravities with model heights, respectively, 
of (75 x 12)/90 = 10 in. (250 mm) at 90 8, OF (75 x x 2)/120 = 7.5 in. (187.5 
Test Equipment, Materiais, And Test Proceoures 
Centrifuge. —The Cambridge Geotechnical Centrifuge has a effective ‘radius 
of 13.1 ft (4 m) from the rotor center to mid-depth of the soil in a model. 
; om model was constructed inside a strong box, with embankment models being _ 


“ 


10 in. high, 30 in. long, and 6 in. deep (254 mm x 762 mm x 152 mm). The 
strong box and model are attached to a swinging platform mounted on the 
end of the centrifuge arm arm. As the: centrifuge is accelerated, the swinging platform a 
_ Swings into position and seats against the end of the centrifuge arm. In this — y 
; : position, the model lies in an horizontal plane with the centrifuge component — 
_ of acceleration directed radially outward. As the model lies in the horizontal 
plane | at the of the balanced rotor arm, observation of the section 


4 
- 


with the video camera ‘located at the rotor and a 45- -deg mirror mounted . 
the strong box (see Fig. 1); or (2) vertical flash photography synchronized as 
model rotated past a stationary camera through the centrifuge- “pit 


CCTV CAMERA 


AAS 


SILVERED 


R INGING 


we 
event Diagram Contrituge Model ‘and CCTV for Viewing 
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—Grain- Size Distribution of Coal Waste Used in Models 


“safety cover. Additional d details concerning the equipment have been presented — 
inRefs.3and7. j ed? baided stow 
_ Materials.—Fig. 2 presents the grain-size distribution of the coal waste materials — 


a in the models. These materials represent the minus 1.18-mm fraction of 2 


a coal w waste embankment material I sampled i in the western United States. — 


ofthe model can be made through a thick transparent face ofthe strong box. 
wie 
— 
3 
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TABLE 1 —Summary of Consolidated-Undrained Triaxial 


At Shear pecimen 


Maxi 


| "stress, 
| oy, 
content, Lie ity, k, in | kilograms 
ow, as pounds per | centimeters | per square | per 
‘ Test identification percentage cubic foot per second | centimeter | centimeter 


Material B; Gs = 2.43: 3: 


14-1 93. 1.67 x 107° 


4 |3.26x WwW 


89.8 | 3.72 x 107* 
30-1 2s | 923 206 


Bilinear Bilinear envelope use used in stability 1 Ib/cu ft - = 16.02 kg/m. 


variation was accomplished by adding various percentages of fines to o the ilies cay 
fraction. Table 1 lists the results of 35.6-mm consolidated undrained triaxial ~ 
7 compression tests (R) on these materials. These results show that permeability | a 
ranged from cm/s- cm/s. The results also show values 
ranged from 26.9-34.7 deg and c " values ranged from 0.025 kN/m? 0. 24 kN/m?. 


Construction 

The sequence of model construction consisted of five basic steps, which 

were: (1) Template placement; (2) material placement and rodding, (3) marking 


sequence as s follows: 
1. sheet templates the slope, crest, and the modeled 
embankment, and different sections such as the key and toe drain, were lightly 
“oiled, and fitted and wedged in the strong box. 
= Subsequently, the sieved mine waste material and coarse drain material — 
a were placed behind the forming template. Since the material had been wet ~ 
_ sieved and was stored under water, during placement it had a slushy consistency, _ a 
and it was rodded not for density control, but rather to ensure that no voids © 


remained behind the material was allowed to 


model; (4) template and (5) Strong, ‘bon detail, the 


f 
4 
— 
— 
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RIFUGE MODELING 
Compression Tests o on Coal Waste Model Material 
from thy foed A rise ig the 


MAXIMUM EFFECTIVE STRESSRATIO.c,/o; 


minor stress, 
| Principal 
kilograms 


per square stress ratio, Strain square o’, in per square 
centimeter percent centimeter | de 2s | centimeter 


10.19 


ws drain as the model and strong box were tilted to a 45-deg position. 1. After ebiiied 
the and model were laid excess material was screeded 


ar: The model surface | 
eads were placed on a 0.6 in. 1.2 in. mm x 30 mm) grid using a template 
to enhance photographic records of movements after testing. 
} - 4. The templates were then carefully removed from their positions and last — 
a minute patching operations undertaken prior to placing the Plexiglas window. -— 
5. The final assembly step consisted of placing the lightly oiled Plexiglas 
window over the model, with its frame, and securing it by bolts. Sealing between a 
the Plexiglas and box is accomplished by a sealant squirted on the edges prior 
The test procedure then involved a delay in which balance calculations were 
completed, a counterweight was S fitted at the other end of the balanced rotor 
arm, and the model package ' was secured to its swinging platform. 
_—s geometries of the models and corresponding prototype and post-test 
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"CENTRIFUGE MODELING 


Model MW-9. —This model was slowly accelerated t to 100 eh. water 
was supplied to the reservoir from the feed pipe. A rapid rise in the piezometric 

_ levels at the base of the embankment was observed in the manometers, and 

: _ the phreatic levels were also clearly observed to rise on the CCTV monitor. 4 - 
of As the reservoir level surged upward, the phreatic level | reached the — 
sketched in Fig. 3 and failure was observed to begin by seepage and erosion 

at the embankment toe. This initial erosion was followed by a succession all 
slip failures, which are indicated in Fig. 3. After these retrogressing failures, : 

i. debris lay as a beach choking the exit and preventing further regression, but 
a. allowing seepage to flow . Additional water was then added to the reservoir 
; which caused erosion of the slope and d breaching o of the ; e model embankment. _ 

Ina an analogous prototype experience, a ce, a coal waste embankment 73.6 ft (22. 

g high would be constructed with a slope of 1 Vertical:1.5 Horizontal, with 

4 _ the reservoir filled with coarse rocks behind the embankment. A small stream © 
= would be directed into the reservoir and water begins to flow out from the as 

embankment toe. Additional streamflow raising the reservoir level would then 
Whes the -eservolr War al Ge cquvelen! of 748-4; (72.91 


PHREATIC SURFACE 
JUST PRIOR FAILURE a) 


PROTOTYPE 


PROTOTYPE TIMES: 


3 HR TO FILL RESERVOIR 
b. 10MIN BETWEEN RESERVOIR 
FILLING AND SLIP A-A 
40 MIN BETWEEN A- 


cause erosion of the toe and diepages of the embankment face which lead lead 
to a stable beach with steady seepage underneath the embankment. A much 
larger streamflow would then be diverted into the reservoir, leading in rapid 
succession to a rise of the reservoir level, saturation of the coal waste embankment 
surface, water flow, erosion, and catastrophic breaching of the embankment. — 
5! Using a scaling factor of N, the time from occurrence of the initial toe erosion a 
; until practically breaching the embankment would be 140 min. inane 
Model MW-10. —This : model was higher than Model MW-9 (as shown in Table 
2) and was accelerated slowly to 120 g. At this stage the reservoir was filled 4 
in a more controlled fashion and attempts were made to determine seepage 
“flow rates. The final increment of flow was perhaps applied too quickly = 
“the phreatic line sketched in Fig. 4 probably does not reflect a steady-state “4 
Sg The full flow from the feed pipe was at a rate about 1/50 cu 7 ; 
— 00056 m°*) in about 26 sec through the 6-in. (150-mm) wide _ model. This is 
equivalent to a seepage of about 2 x 1/50 x 120/26 = 0.185 cu ed ; 
68 0172 m’*/s) of prototype embankment. _ The full flow from the feed pipe led 
to toe « erosion and retrogressive ) oe sketched in Fig. 4. The location of the 


| 
theres 
Thee - 


_ two major slip cutbiie BBs and cc closely repeated the result of the test 

model MW-9. bingy A ed? mor) si? 0} beliqags anw 
In terms of an analogous prototype experience, a coal waste 
is built to a height of 94 ft (28.66 m). Flow of water from a reservoir through 
Dues PROFILE JUST PRIOR TO FAILURE do hak £ a bad 


_- CALCULATED FAILURE SURFACE 


ant oni 
_ 
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60 FT AND HEA 
ne. 4.—Schematic Showing Location of Phreatic and and 
lated Failure Surfaces for Model Test MW-10, 


od? 


tew 

FIG. 5.—Model MW-11: (a) Prior to Centrifuge Test; (b) After Centrifuge Test 


the embankment leads to erosion at the tc tos and retrogressive slips when ‘the 
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age 110W IS aD0Ul 16.9 CU SEC 5) per TU Ti) OF Provolype 
_ be! embankment. The time between the initial observance of distress at the toe — | 


CENTRIFUGE MODELING 


| a Model MW- 11.—This model was constructed of material B (see Fig. 2) usin 

ie a diagonal (sloping upstream face) r reservoir filled with coarse waste meteors 
The model was accelerated to 120 g, which corresponds to a notional prototype 
95.3 ft (29.07 m) high. Fig. 5 shows the model before and after the test. ‘The 
line of the original slope template is clearly visible in Fig. 5(a) and the dark ve 
mass of material at the back of the box can be seen to have settled by about 


4 Ss by the time full acceleration was reached . Supply of water to the reservoir 


| 1. When the reservoir was at the equivalent of a 51-ft (15.56-m) head a flow 

fate to the model of 1/50 cu ft (0.00058 m’*) in 96 sec was measured. ‘This ; 
tate corresponds | to a flow through a notional prototype of about 


2x 


— x = 0.05 cu ft/sec/fe (0.00483 m?/s/m)............@Q) 


2. When the reservoir was at the equivalent c of 74.8- 74.8-ft ft (22. 814m) head reed 
‘flow rate to the model was measured as 1/50 cu ft (0.00058 m *) in 29 sec, 


= 0. 17 cu ft/sec/ft (0.016 m’ /s/m) 


this stage, toe erosion pa ted | to a series of successive retrogressions, 

; indicated in Fig. 6. After each slip, the accumulated debris formed a oer 

a slope of al about t 10 deg, with steady seepage from the ‘Teservoir onto the 

7 beach. It is significant to note that the steel surface of the base of the box — 
had not been roughened, and sliding may have occurred along this surface. 
The times of retrogression scaled by the factor N = 120 are also listed in 
Fig. 6. This factor N = 120 is based upon Malushitsky’ s (4) “‘modeling of 
models’’ and has not been confirmed. The rate of failure was accelerated S| 
raising the reservoir to full height, after which the Tetrogression of erosion 


and slip eventually led to a catastropic breach of the re reservoir wall and ‘a 


____ In analogous prototype experience, a coal waste embankment with . slope of 


MODEL PRIOR 


FAILURE 


PHREATIC SURFACE 
JUST PRIOR TO FAILURE 


FIG. 6.—Schematic Drawing Showing Location of Phreatic and Observed and Cal- 
culated Failure Surfaces for Model vem ‘1 what 


— 
| 
= 
| 
4 4 j 
ding 
| 
= 


7 = 1V22H is b 2H is built to a height of 95.3 95.3 ft (29.07 07 m). A steady dint of water from 
the reservoir through the base of the embankment leads to erosion at the toe — 
% when the reservoir is only three-fourths full. When the reservoir level is raised a 
to full height, “rapid re retrogression leads to a breach of the reservoir wall. The - 
time between initial observance of toe failure and catastrophic breaching of 7 
the embankment wasSdays, 8 = 
4 Model MW-12.—This was a model of material B, with a diagonal (sloping 
upstream face) reservoir, and an inclusion of ‘‘key’’ material of coarse — 
(2.36 mm-12.5 mm) at the center of the bank. Such a “‘key”’ feature had been 
included in the original waste embankment from which the model material had 
_ been obtained. The model was accelerated to 120 g in similarity with a notional 
; ? prototype 95.3 ft (29.07 m) high. Reservoir filling proceeded slowly over a period © 
equivalent to 2 days, by which stage the reservoir level was only at an equivalent . 
prototype height of 63 ft (19.22 m) when toe failure began. The reservoir supply 
was interrupted t to measure a flow w of 1/50 cu ft (0.00058 in ‘sec, 


= 0.067 cu ft/sec/ft (0.00644 m’° /sfm) .. 


m4 


(PROTOTYPE) 


— 


the key material allowed seepage failure to be initiated and to progress xl 
a lower reservoir head. The phreatic line just before failure, and ad initial 


slip surfaces are sketched in Fig. 7. 


In analogous prototype experience, a coal waste embankment 
be constructed to a height of 95.3 ft (29.07 m) with a 1V:2H slope. A rapid 

raising of the reservoir pool to 63 ft (19.22 m) in 2 days would cause seepage 

through the embankment, initiating a toe failure. . Thus, an with 


s comparable embankment without a key experiences erosion when three- tenets 
full. Inability to remove slope debris halts further erosion. = 
_ Having thus established that a notional prototype could not retain fulla reservoir _ 
_ of water, it was decided to explore the effects of a remedial toe drain to a ll 
Model MW-13.— —This model was constructed of material (see Fig. 2 


5 had a toe drain constructed with about the same grain size as the ‘ — "material 


q 
| 
| 


However, this toe por was to jadi 4 a more useful function; i.e., all 
high seepage to pass safely without damage to the toe of the slope. The model i 
was accelerated to 120 g. Water was subsequently supplied to the reservoir 
and the phreatic line rose higher than ever before; reaching an equivalent j prototype — 
height of 90.55 ft (27.62 m), at which it remained stable for the equivalent 
of 7 days without any distress at the toe of the bank. At this stage the seepage os 
flow of 1/50 cu ft (0.0 00056 in was toa prototype flow 


“a ‘a In an analogous | prototype experience, a a coal w waste ‘embankment containing 7 ; 
undrained ‘‘key’’ of coarse material, and downstream drainage control via 

a toe drain would be constructed on a 1V:2H slope to a height of 97.3 ft 

(29. 68 m), at which condition it is stable. Subsequently, reservoir filling could BY 
_ be initiated and full reservoir capacity achieved. With this full reservoir head, 

the embankment could successfully allow the | passage of seepage of 0.2 cu 

ft/sec/ft (0.0019 m*/s/m) through the toe drain and remain stable. 7 
_ Model MW-23.—Model MW-23 was constructed of material corresponding | 

“to gradation MW-23 in Fig. 2 and containing 8% fines. The model was accelerated 

to 100 g, and then the speed reduced to various levels, with different flow 

_Tates, to establish data to determine the permeability. At a throughflow rate 
1o~° cu ft/sec/ft or 112 mm*/s/mm of model width at 75 g, erosion 
occurred. Once initiated, the erosion was aggravated by the concentration of ; 

_ Water own to that point; as the centrifuge speed was reduced from 100 g 
4 66 g, » and the water supply ultimately increased to 3.54 x 10~° cu ft/sec/ft 
(329 mm’/s/mm) of model width, the erosion reached a clearly ‘unacceptable 
level, forming a uniform blanket of debris across the model. 
_In analogous prototype terms, this model simulates several notional | prototypes 
of different heights, permeabilities, and throughflows, but with constant strength Jf 
Parameters. A prototype embankment of height 60.7 ft (18.5 m), with permeability . 
- equal to 44 mm/s, impounding a reservoir 37.7 ft (11.5 m) deep, i.e., 62% 

) ae full capacity, and permitting a throughflow of 0.12 cu ft/sec/ft or 1.1 x 
10-7 m’/s/m of embankment width is shown to be stable. By simultaneously 

"changing height and permeability of the embankment for the same model = 

an embankment-dam passing a throughflow of 0.8 x 10-7 m°/s/m 

' of embankment width from a reservoir 29.9 ft (9.1 m) in height or at 66% 


4 


of full capacity is found to be unstable with respect to erosion. Once begun, 
the weakened point attracts more re water, and the final eroded profile occurred 
when the notional prototype was an embankment-dam 39.4 ft (12 m) in height, — 
with a prototype permeability of 9 mm / S, at full reservoir height passing a 
throughflow 0.2 ft/sec/ft or 2 x m*/s/m of embankment width. __ 
Model MW-24. 24. .—This model was prepared of material containing 34% fines 
(see Fig. 2) was accelerated to 9 g over 90 sec, and left at that speed for 
13.5 min. It was further accelerated, and immediately upon achieving 15 g, 
q failed in a conventional rotational slip circle re 8), and the intact wedge — 


of debris reclined at 16 to the 


— 
| 
| 
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In terms of analogous prototype experience, this model simulates a rate of 
= construction failure. A prototype embankment with permeability equal to 3.8 
mm/s is built to a height: of 5.6 ft (1.7 m) over a 30-min period and then 
left to consolidate for 18 h. (Since embankment failure resulted from excessive 
rate of construction creating high pore pressures, the governing time factor 
comparing prototype with model is N’ (3). Thus (13.5 min) x 97 + 60 = 18.2 4 
__h.) When, however, an attempt is made to increase the height of the embankment 
to 2.8 m, the slope immediately fails upon reaching that height. ae ete aad 
7 7 Model MW-25.—This model of material containing 9% fines (see Fig. 2) was 
_ prepared identically to MW-23 but with a much wider vertical upstream face 
~ due to the decreasing supply of model soil. The model was accelerated to 100 — 


g over 22 min, after which time various combinations of throughflows and 


g, with 
a a flow rate of 1.3 x 10~* cu u ft/sec/ft o or 1218 mm am?/s/mm. of embankment width, 
an slope was stable; for the same model flow rate, erosion began to occur 
2 : when the speed of the centrifuge was reduced to 30 g. The erosion was permitted — 
to continue at 30 g for 4 min. dade te fr 


FIG. 8. —Model MW-24: Rotational Slip 2 at 15; Relative to State at ato g 


In analogous prototype experience, as in model MW- 23, several different 
_ prototypes are simulated. An embankment of height 61.0 ft (18. 6 m), aan 


was 15.4 ft (4. 7 =) or at odes capacity, pan equivalent throughflow of 0.04 
cu ttn or 0. 4 x 10°? m */s/m of embankment width caused erosion at 


_ Six uniform prototype embankment sections were analyzed by conventional _ A 

_ techniques to locate the phreatic surface and estimate isotropic ne Dupuis # 
of the material. The Phreatic surface location: was calculated Dupuit’s 


! 
a 
= a 
cu ft/sec/ft or 1.2 x m’/s/m of embankment width, was found to 
. q _ Stable. However, simultaneously decreasing the embankment height and per- 4 


of ‘seepage surface or distance f ‘from exit point to toe; and a a= 


prototype embankment sections (models 9, MW-10, MW-23, and 
MW-25) have vertical upstream reservoir faces, while models MW-11 and MW- od 
100 = = RESERVOIR HEIGHT AS PERCENT a 


valvet 


_ have diagonal (sloping) upstream faces. The phreatic surfaces for each embank- 

_ ment section were calculated for 100%, 80%, 60%, and 40% of the maximum 


reservoir level or crest elevation, as shown in Figs. 9-14. Also presented in 
_ these figures are the phreatic surfaces as observed in the CCTV monitors or 
photographs of the centrifuge models. For embankments with a vertical le 


_ face, the observed phreatic surface closely resembles the calculated parabolic _ 


9 
9.—Comparison of Observed and Calculated Phreatic Surfaces for Centrifuge 
— 
— 
q 


surfaces, exceptions s that the modeled surfaces are more than 

_ Parabolic, their gradients are much steeper, and they intersect the downstream — 7 
‘Slope’ at lower elevations. _ However, in the case of the model embankments 
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STEADY SEEPAGE THROUGH 
FIG. 12. —Comparison of Observed and Calculated Surfaces for 
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— 100 = RESERVOIR HEIGHT AS PERCENT OF : 
EMBANKMENT 


FIG. 13. of Observed and Calculated Phreatic Surfaces for Contrituge 


cutoff, which is reflected as a free-fall tendency in the phreatic surface. For = 
: this reason, the observed centrifuge model phreatic surfaces deviate from the 
- calculated values as shown in Fig. 11 for model MW-11 and later in Fig. 12 
for model MW-13. In an attempt to understand the effect of the vertical cutoff 
on the | phreatic surface of models with diagonal reservoirs, a a two-dimensional 


FIG. 11.—Comparison aces for 
7 
| 


CENTRIFUGE MODELING 


finite: element method (FEM) seepage analyses ( (8) 1 was as performed. Mw- ll 
_ with a diagonal reservoir was selected as representative, while MW-23 with — 
a vertical reservoir was selected for comparison and control. The results of vr 
these FEM analyses are also presented in Figs. 11 and 13 and show that the 
FEM phreatic surface | essentially duplicates Dupuit’s (1) parabola for both 
ss: The coefficient of permeability per unit length, k, can be ‘ditecthente cd sing 
RS the flow net and the measured yin of seepage per unit length, 


in which q = flow rate /unit length the embankment; hN= “height 
of prototype 1 reservoir; and N,/N, = ratio of equipotential drops, Nu. a> to the 


_ number of flow channels, N,, as obtained from a flow net. PEE ra 4 
-100= R RESERVOIR HEIGHT AS PERCENT OF 
ad) 


hy a atthough only appropriate for clean coarse sands and gravels, Hazen’ s Equation 


in which k= permeability in cm/sec; and D,, = diameter, in centimeters of 


= (10% size. Accordingly, Table 3 presents estimates of permeability values based a : 
upon the aforementioned three methods for various models. A cursory examina- > 
: = of the data shows that for models MW-10-MW-12 permeability values were 
in the range of 10~°, with excellent agreement between flow net based estimations 
and phreatic surface based estimations. Likewise, Hazen’s Equation provided 
approximately the same values, which lends some confidence to the values 
By comparison, permeability \ values previously presented i in Table 1, obtained © F 
_ from triaxial compression tests, are in the range of 10~° for models of material 
2B, and 10~* for models 23 and 25. Thus, a fairly poor agreement was obtained ; 
: between laboratory and model k values. This lack lack of of agreement may be attri be er 


| 
a 
| 
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“a. 


TABLE 2 of Coefficient of Perme 


Flow rate, mea- 
sured in model, 
centimeters per centi- | height, height, 
millimeter per second | meter length 43 meters a 


1.02 x 107° 
1.02 x 107° 


os *k= q/(a sin a tan a). 


q/h-N x 


~ to leakage occurring along the base or sides of the model, resulting i ” eo 


It is important in cen not results with prototype 
events, but also with numerical forms of analysis. This section examines the 
factors of safety (FS) recommended for coal waste embankments (2), with those — 
re of the various models. In this fashion, the appropriateness of recommended — Pa 
a Stability analyses were performed utilizing the Corps of Engineers \ Modified — 


- ¥ - Swedish Computer Program, as described in Ref. 6. Strength parameters utilized 


in these analyses were previously in The 


results are listed in the — 


Observed 


4 
“These FS ‘albeit conservative, phenomena ‘observed in the 
centrifuge models. That is, models MW-9, MW-10, MW-11, and MW-12 were 


unstable and felled, while models MW-13, MW-23, and MW-25 were 7 


al The conservative nature of the triaxial strength parameters 


not considered by these analyses. 


is evident, and the recommended FS of 1.3-1.5 (2) is quite appropriate. es 


Test 
number 
gher 
ner 


“ability, k, tor Coal Waste Conttuge 


=a toe to seep- a k, in centi- aN ® Shape k, in centi- 


age exit, a Gravitational meters per 
in meters acceleration, N 


100° 0. ne 


_ to MW-11 with the critical arcs estimated by the stability ; analyses are shown — 
in Figs. 4-6. These critical arcs agree quite favorably with the failure surfaces 
_ observed in the models, which lends confidence to the FS derived by the analysis. ne 
« ‘The effect of the toe drain (model MW-13) intercepting the phreatic ian 
a, and preventing seepage on the downstream slope obviously increased stability 7 


of the embankment. With the reservoir at m the FS almost achieved 
recommended criteria. : 


All of the the previously examined centrifuge model tests were subjected to seepage q 
except model MW-24, which failed prematurely while accelerating the model 
to prototype acceleration. This type of failure corresponds to a rate of construction 
failure that occurs when high pore pressures are created by the fill placement 7 
being too rapid. The prototype dimensions at failure are presented in ~ 
2, while Table 1 presents the parameters used in the stability analyses. = 


. In analyzing stability for the end of construction case, the ne Corps of Engineers’ 


_ effective or total, is used, resulting in a bilinear envelope. Despitet this conservative | 
: approach, the results of the analyses provided a f a factor of f safety v value of 3. 60, 7 


indicating stability; yet failure occurred. te tine boot 
_ Based upon the centrifuge models and matociale tested, the following conclu- | 


- = . Estimates of the phreatic st surface based upon Dupuit’ s parabolic assumption — 
and FEM seepage analysis agreed | well with phreatic surfaces observed in models — 
% with vertical upstream slopes. However, in models with more realistic sloping i 


upstream surfaces agreement was less favorable. eee: ons both upstream > 


Be 


q 
7 
| 
i 
eis itty Ulal fie >, a> pu pelicia 
unknown for this case. However, as it is assumed that negative pore water — / 
| 
| 
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_ slope conditions, theoretical estimates of the phreatic —_ were higher (and 

thus more conservative) than surfaces observed in the models. It is possible J 

= steady-state seepage may not have been obtained for models. a ee 

_ 2. Permeability estimates based upon seepage quantities measured through ; 

es models were of the order of 10~* cm/s, while laboratory permeability values _ 
"were of order of the 10~° cm/s. This difference may be due to 

permeabilities of the material or leakage along the sides of the models. = a 
i a * . Stability analysis using the Corps of Engineers (modified Swedish) slope a 
- stability method d agreed well, albeit conservatively, with observed model perfor- 4 

mance. Factors of safety ranged from 0.7 to 0.9 at ‘the instant of failure | for 
_ models \ with retrogressing slips. By comparison, factors of safety for stable 
embankments ranged from 1.3 to 1.7. Thus, from the standpoint of stability, 
the recommended minimum factor of safety of 1.3 is sufficient. = = | 

4. For embankments of cohesionless materials, permeability can govern the 
: failure mode. In cases of low permeability, rate of construction may pose a 

problem, particularly when the initial water content is high, as in hydraulically 
an fill. In case of intermediate permeability, pore pressures due to seepage ~ 

can cause mass instability and retrogression. In cases of highly permeable material, 

or overtopping, erosion may pose a stability hazard. tw 
Exiting of the surface on the downstream slope generally leads 
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The following symbols « are used in in this paper: 


= length of seepage surface along Seem slope, 
c undrained, drained cohesive strength; 
= hydraulic gradient; 


= number of flow channels 
a = slope angle; 

= unit weight of wai water; 


unit weight of solids (soil); 


5 = deviator stress. 
wees view af shown i ie, 
by tie 
waa 


6-1/2 of the 
of 
inn 


| 
| 
q 
| 
— 


of safety 


4 


"ENGINE EERING 
NEAR 
AMATITLAN, GUATEMALA 


y HE Bolton Seed," F. ASCE, Ignacio Arango, Clarence K. Chan, 


Iwtropuction 


an liquefaction ‘which occurred at ‘the settlement of La Playa o on 1 
pos prog shore of Lake Amatitlan. The settlement is located on a deltaic deposit 
near the mouth of a small river, Rio Villalobos, and consists mainly of a resort 
area with several newly constructed brick and cement-block houses surrounded 
by grass lawns and trees. An overall view of the area is shown in Fig. A 


about away from ‘the lake. The “ground so some distance from 
the beach is devoted to agriculture. The area of extensive liquefaction extended . 
_ about 600 ft (180 m) behind the nearest shoreline and is marked by the shaded -_ 
 . in Fig. 2. A second zone, where evidence of liquefaction was readily — : 
apparent but its effects were less severe, extended about 1,200 ft (60 m) behind | 
An excellent description of the geology of the delta and the nature of - - 
_ damage has been presented by Hoose, et al. (2), and Krinitzsky and Bonis 
(4). As described by Hoose, et al., the Villalobos delta covers an area of about — 
16-1/2 km? of the northeastern portion of the Amatitlan basin, but only i: 
6% of the area was affected by liquefaction-induced ground failure during the 
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Zone of extensive liquefaction en. 
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1 Overs liew of Liquefaction Zone at La Playa 
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976 earthquake. The youngest, most active portions of the delta were ol 


most severely affected. The river has changed its course at least three times 

in the | past § 80 yr, and the delta has had a rapid rate of deposition so that = 
it now extends well out into Lake Amatitlan. The material deposited by the — " 
‘fiver is primarily pumice sand and gravel derived from extensive Quaternary 

Drs reser to the north of the lake. Detailed maps of the ground failures and 4 


site inspection an and aerial photographs. _ The zones | Shown ‘in n Fig. 2 ar are based 


evidence of liquefaction have been presented by Hoose, et al. (2), based on 


inspections by three writers of the present paper. 
‘Within the heavy damage zone (see Fig. 2) there was subsidence and flooding 
of beach areas, severe a nd cracking with cracks ranging in size ™ to oases 
ro 


4 


me FIG. 3 —Ground Cracking in Zone of Extensive Liquefaction | 
meters (Fig. 3), severe damage to houses where they were located in areas ~ 
of extensive ground cracking, as well as numerous sand boils. The ground cracking — » 
was generally more severe near the waterfront with crack widths decreasing _ 
progressively in size as the ground rose s¢ gradually i in elevation towards the inland i 
edge of the liquefaction area. 
_ Of the 32 houses at La Playa, 29 w were destroyed or 5 denen by differential 
lateral displacement, generally as a consequence of lateral spreading and subsi- — 
dence. Fortunately, most of the houses at La Playa were mainly vacation Te 
they were unoccupied at the time of the earthquake. Ov, Juods 
An! interesting feature of the | area was the expulsion, together with sand, oa 
of occasional pieces of pumice, ranging from small gravel to cobble sizes in 
some of the mounds formed around the sand boils. The pumice particles \ were 
«of wey light — Sa that they were readily carried up with the liquefied 4 


q 
= 1 
| AST | 
Bs, Oe — 
| 
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; sand through the vents and cracks through which water r moved upwards 4 
the ground surface. Areas of sand boils in the zone of moderate liquefaction 
+ In contrast to the extensive damage due to liquefaction in the La Playa area, 7 
ie damage to houses in nearby towns was slight. Even a cement-brick house located 
- near an area of intense sand boils and land spreading at La Playa was virtually 
- undamaged (4), and the shaking intensity in the area, based on this lack of 
damage, was assessed as VI on the modified Mercalli scale. 
_ In view of the scarcity of well-documented cases of soil liquefaction it was : 
considered important to determine the soil conditions within and immediately al 
adjacent to the liquefied zone in order to supplement the available empirical 
data base used for liquefaction evaluations of other sites. At the ‘s same time, 7 


Boils at st Ground 
¥ it was considered useful to investigate whether the liquefaction that ecvesed. 
at La Playa might have been anticipated the use of corrently used 


Guaremaca EartHauae or Fesruary 4, 1976 Ya 
Z The February 4, 1976 Guatemala earthquake occurred on the Motagua Fault ; 
and has been assigned a Richter magnitude of 7.5. Its epicenter (see Fig. 5) 
was located about 170 km northeast of Lake Amatitlan, but su surface ‘Tuptures 7 
along the fault were mapped over a length of at least 240 km and the western 
end of the major east-west fault break was only about 40 km due north of ~ 
_ the lake. In addition, some short length north-south offset faults were mapped a 


in and around km north of the site. The distribution 


yy 
4 


of modified ‘Mercalli intensities in the strong motion area of the cechats, 
as assigned by the United States Geological Survey (USGS) (1), is shown in 
| Fig. 6. Lake Amatitlan is located in Intensity Zone VI but very close | to the 
_isoseismal for Intensity VII. Based on the correlation of intensities with p peak — 
a ground acceleration developed by Trifunac and Brady (9), this would indicate 


°ALVADOR 


5.—Location of Lake Amatitlan in Relation to Motagua Fault Espinosa) 


Intensity Distribution 


° 


FIG. 6.—Distribution of Modified Mercalli Intensities from 1976 Guatemala Earthquake ee 
a peak ground acceleration in the vicinity of Lake Amatitlan of the order of 


_ Data on the damage to adobe houses during the earthquake has also been 
assembled by the USGS (1) and Presented i in the form of a damage distribution 


| 
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damage at La Playa would seem to confirm this ; ssiaaaiaane t and suggest that * 
the peak ground accelerations were nearer the Trifunac and Brady mean value : 
4 of 0.12 g rather than a significantly higher value which might be indicated by Py 
probability considerations within the available intensity versus peak acceleration - 
Peak ground accelerations generated by the magnitude 7.5 earthquake at Lake 
= can also be assessed using the correlations between peak acceleration, oe 
basic magnitude, distance from causative fault and local soil conditions 


Lake Amatition 

FIG. 7. —Distribution of to Ad Adobe Houses in in 1976 Guatemala (1) 


Typical attenuation curves — 
for M = 75 earthquokes 
(from Ref. 4) 
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FIG. 8 of Laws t to Estimate Peak Acceleration at Lake Amatitian” 
: proposed i in Refs. 5 and 8. The correlations would indicate a peak acceleration 
in rock near Lake Amatitlan of the order of 0.21 g, but for the softer soil 
conditions near the lake this value would be expected to be reduced to about — a 


0. 16 g. A similar value is indicated by the correlations "presented by Idriss o 


a. Unfortunately, there were no records of the grovad motions induced by the 
earthquake near Lake Amatitlan. However seismoscope records were cotained — 
in Guatemala City (about 25 kn km from the causative ae and an interpretation 
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of these by the USGS (1) indicates that the peak ground eaten neato, 
in this area was probably of the order of 0.25 g. At the increased distance 

_ (40 km) to Lake Amatitlan, the use of attenuation laws (see, e.g., Fig. 8) would 
indicate a peak ground surface acceleration in the vicinity of Lake Amatitlan 
The lack of damage to adobe houses in the vicinity of ‘the lake wo would seem 
to indicate that the actual peak acceleration in this area was nearer the lower _ 
bound of the aforementioned values, but since the range indicated by the dif! ferent 
procedures is compratively low, it seems reasonable to conclude that the peak 

4 ground surface acceleration developed near the lake was probably in the range 7 
of about 0.12 g-0.15 g and the behavior can be assessed based» 
on this range of of “4 


= 


The soil conditions i in the area wher 


ire 2. ‘One boring (No. 2) just outside the boundary of the liquefied zone; : this 
boring was especially carefully located, on the basis of field age epee. 
and mapping by Hoose, et al. (2) about 30 ft (9.15 nearest 


3. Two borings well inside in | the ‘nonliquefied zone (Nos. 3 and 9 
mS The locations of the borings are shown in Fig. 2. It would pan i been 
F desirable to have boring data in the liquefied zone before the earthquake. However, — 
_ it was reasoned that by making borings as described in the ‘preceding, it ji 
be possible to establish the soil conditions r representative of a marginal liquefaction 
3 - condition for the ground motions resulting from the 1976 earthquake. = 
a Borings were made with a hollow-stem auger using drilling mud to stabilize 
the sides and bottom of the hole. In each boring standard penetration tests 
_ were performed at intervals of 5 ft (1.53 m), and “‘undisturbed’’ samples were - 
exuacned | in thin-wall seamless steel tubes for cyclic load testing in the laboratory. a 
‘Standard penetration resistance values were determined using a rope and pulley & 
‘system (two turns of old rope around the pulley) to raise the 140-lb weight 4 
and measuring the number of blows per foot penetration at the bottom of a 
i” drill hole whose sides were supported by drilling mud. hk ee 
£ addition, two test pits were dug t to a | depth of about 5 ft (1.53 m) on 
the level ground about 40 ft (12 m) from ‘the beach. These showed the —— ; 
a4 ft (1.2 m) of soil to be pumice sand underlain by a thin layer of dense Gem 
¢ peats and organic silt. The soil exposed on the walls of the pits was a stratified 
brown fine sand with some silt. One pit was dug across four cracks with a 
. total width of 25 cm; the ground surface showed several centimeters of vertical © 
_ Subsidence towards tl the lake. The same vertical differential movement — 
@ bes seen in the walls of the pit and increased with depth. The second pit was 
dug across a 10-cm wide crack showing no vertical differential subsidence ~ / 


the crack was as partially filled with sand which had t been forced up from below, 
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wet density of 55. 5 Ib/cu ft-61.5 Ib/cu ft, a water content of about 56% oe 
a dry density of 35. — ft-39.4 Ib/cu ft. These tests were taken in the — = 
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FIG. 10.—Schematic Soil Profile » Through La Playa 


; ‘moist soil above the weter table. It may be noted that pumice sands were 
also present south of the La Playa area and there were other small areas of — 
liquefaction but these were not investigated in the present study. ‘Seah gab 
— ee to a einai of about 70 ft (22 m) is shown by the logs of 


| 
and had intersected the horizontal stratifications. Two field density tests in {Jf 
| 
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a a section extending, from: the shoreline to a distance about 600 ft (488 
ee _ It may be seen that the entire area is covered by a surficial layer of a _ 
- pumice sand varying in depth from about 5 ft-10 ft (1.5 m-3 m) underlain 
in the zone of liquefaction by medium coarse sand containing some a 
a fragments and in the nonliquefied zone by an intermediate layer of black clayey . 
silt above the same type of sand with pumice particles. The water table varies _ 
in depth from about 4 ft-5 ft (=1.5 m) near the shoreline to about 12 ft (~=4 
m) depth in the area where no liquefaction was observed. It is interesting to - 
gui i 
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gray to black sand 
with pumice particles 


aa 


very buitiveqs > 
chan: FIG. 11 —General Soil Characteristics at of Liguetestion 
e note that the very lightweight pumice sand at all boring locations was located ~ 
_above the water table and thus was not itself the zone in which liquefaction 
occurred. The source of liquefaction was presumably in the saturated sand 
4 containing pumice particles below depths of 6 ft-10 ft (~2 m-3 m). In the 
heavily liquefied zone this material had a penetration resistance value, N, of 
oe about 5 ft (~1.5 m) to a depth of 50 ft (~15 m) but presumably ‘this 
value is representative of the liquefied and restabilized sand. There is no a 
of establishing the penetration resistance of this material before the earthquake. 
> For this reason, the remaining borings were made in the nonliquefied zone, 
cal 


one very close to the boundary (Boring No. 2), and two well behind the liquefaction 
zone. With increasing distance the shoreline, the average 


which would tend to reduce the possibility of liquefaction of the sand deposit. 
“ examination of the soil profile in Fig. 10 indicates that the conditions 
at Boring No. 2 are probably representative of the limiting conditions at which — 
aa would just occur or just not occur for the ground motions and 


soil conditions in the La Playa area and this profile has therefore been used | 


for a more detailed study. The general soil characteristics are shown more 
_ clearly in Figs. 11 and 12. The latter two figures show the grain size distribution 

_ curves for nine samples taken from depths between 30 ft and 45 ft (9 m and — 

14 m) in Borings No. 3 and 4. Similar soil types were found below the water 

table in Borings 1 and 2 where liquefaction occurred. It may be seen that the 

sand ranges from fine to coarse, though most of the samples are medium coarse _ 

sand. However, because of the presence of pumice particles ranging in size ‘ 
f rom sand t to gravel, the sand tested 
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ne. 12.—Grain Size Distribution Curves tor aaeie of Sand Representative of 


‘tics: a average dry | density =58 pef (0.93 g/cc); av average gravity. 


_ of solid particles ~2.1; average water content of sand ~60%; and average density - 
of wet sand ~92.5 pcef (1.48 g/cc). For the generalized soil profile shown in 
Fig. 1, the saturated density was sim _ to a rounded value of 90 pet 
Ana ysis Basep on Empiricat “wale? 
Pe One of the methods currently used to assess the liquefaction potential of 
a sand deposit is a procedure based on an empirical correlation between the 
stress" ratio. induced by the e earthquake and the standard penetration 
_-Tesistance corrected to an effective overburden | pressure of l ton/sq ft (or 
kg/em (6). proposed empirical curves, based on analyses of other 
‘sites, which separate areas where liquefaction has occurred in previous earth- ; 


and the Genin Of the water tanie increaced tore 
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quakes | from which have n not are shown in in ‘Fig. 13. 

for determining the standard penetration resistance, N , , atan overburden pressure © 

of 1 ton/sq ft are shown in Fig. 14. The data for the Lake Amatitlan site 
‘may be used to ee the gunmen of the the curve for for magnitude 7-1/2 earthquakes a 
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test conditions showing liquefaction. 
Open points indicate sites where “ 
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FIG. 13.—Correlation Field Behavior of Sands for L Level Ground 
- Conditions and Penetration Resistance (Supplemented by Data from Large-Scale a 
ei For borderline conditions of liquefaction and no liquefaction at La Playa, 
a the soil conditions are shown in Fig. 11. From this it may be seen that liquefaction _ 
_ is most likely to have developed in the lower blow count zone, N = 8, extending = 
F from about 8 ft-22 ft (2.44 m-6.71 m) below the ground surface. . ad 


this depth, of = overburden pressure show 
= 4x 58 +4 x 70+ 7(90 62.4) 704 psf (or 0.35 kg/m” 
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pressure. For o,= = pef (0.35 ks the value of Cy = = 1, 6, and thus boil 


The c he cyclic Stress ratio induced by the earthquake may be estimated from the | 


in which q 
r, ata of 15 is about 0.95; and o, = 4 x 58 + 4x 70 +7 x 90 


a Since a,,, x was estimated to be in the range of 0.12 g-0.15 g, the probable - 


= value of is about 0. .12-0.15. ‘These values are plotted against the 
q corresponding value of N, 1 in Fig. 15, together with the boundary curve for 
a magnitude 7-1/2 earthquakes taken from Fig. 13. It may be seen that the data 
for the La Playa area are in excellent agreement with the proposed empirical _ 
curve. It is readily apparent that small changes in the values of unit weights — 
used in these computations could influence the positions of the plotted points 
but they would not t change significantly. Thus, the data obtained from the field 
ial “study would seem to provide further corroboration of the position of the empirical 4 
curve for potential in magnitude 7-1/2 


Lasoratory INVESTIGATION 


aaa by comparing the ‘earthquake- induced stress ratio with that determined 

to cause liquefaction or cyclic mobility in laboratory tests on undisturbed samples. _ 

_ During the course of the present investigation, eight cyclic triaxial compression a 
_ tests were performed on samples taken from depths ranging from 30 ft-45 ft — 
1S 15 m-13.73 m) in Borings Nos. 3 and 4. Samples were obtained in thin-wall 
seamless : steel tubes which were f pushed hydraulically into the soil at the base 
of the hole. In the laboratory, the samples were saturated prior to testing. 
a these zones the standard penetration resistance was about 13 to 16, corres- 
7 ponding to an N , value of about 18, so these samples are likely to show somewhat 


a _ higher res resistance to liquefaction than the soil in the depth range of 8 ft-20 
_ The results of these tests, showing the number of cycles required to produce ~ 

7 residual pore pressure ratio of 100% are plotted in Fig. 16. The scatter is 
7 due to the natural variability of the samples, both i in characteristics and density, 7: 


the point ir in the ‘test where the pore ‘pressure reached a value of 100%, — 
_ deformations so the data in 16 be 
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_ For a magnitude 7-1/2 earthquake, | the number of equivalent uniform stress. 
cycles is likely to be of the order of 15 and the field stress ratio, i 


: tests (6). Thus, from this test data the : stress ratio > required to « cause liquefaction :: 


induced stress ratio at a depth of 35 ft (10.68 m) mer! be 


é 


as am | 
loa 
FIG. 16.—Results of Cyclic Loading Triaxial Compression Tests on Undisturbed 
as described previously and shown to have have the f following val values: 
This would i badlente a factor of safety ranging between about rk 15 and | 1.5 which 
_ seems to be in accord with the observed absence of f liquefaction in this area. 
It should be noted that this result may well be the effect of a 
errors. The “‘undisturbed’’ samples were probably densified to some extent 
eat the sampling and handling process. At the same time they probably 
g some resistance which they had acquired due to the prior stress nee 4 a 


(sustained loading plus prior seismic loading). As a result, the measured strength 

-may have been not too different from the in- situ: strength but this would not > 

=" be the case for all such studies, especially those involving very 4g 


Since it is impossible to determine the pre- p-certhaneke | properties of of the sand es 
; in the zone of extensive liquefaction near the shoreline, it is not possible to 
make a detailed analysis of its potential behavior either by empirical methods 
or by analytical- experimental methods. However, some indication of its liquefac- 
tion potential may be obtained if it is assumed that the liquefaction characteristics ual 
; of the sand in this zone were no better than those of the sand in Boring No. 


2, just outside the zone of liquefaction. Conditions near the shoreline ne woud 


In the zone of extensive e's water was at an 
os depth of about 5 ft (1. 53 m) below the ground surface. At a depth of 15 ft ; 


Anes 


~ 6) 


| 


‘For vi values. of a,,,, Tanging from 0.12 g-0. 15 g this leads to values" of cyclic © 
_ If the sand at a depth of 15 ft (4.58 m) had the same penetration resistance 
= the sand in Boring No. 2, then the cyclic stress ratio required to induce ‘ 
liquefaction would be of the order of 0.12-0. 15, thus indicating a factor of y 
* = against liquefaction i in the near-shore area of only about 0.63-0.97. Under : 
these conditions the extensive liquefaction w which actually developed should not 
- be considered surprising. In fact, if the density or penetration resistance of 
the sand in this zone were also somewhat lower than that further away from 
the shore (near Borings 2, 3, and 4), the situation would be even more unfavorable. — 
alternative approach to assessing the resistance to liquefaction of the — . 
: i in the near- shore area is to use the laboratory test results i in conjunction 


— 


_ for the sand in Boring No. 2. Up to values of N, of about 30 or 35, the 
“yesistance to liquefaction of a sand seems to be essentially directly proportion 
to N, (see Fig. 13). Thus, since the stress ratio causing ngage in the 
laboratory tests corresponds roughly to sand for which N, = 18 while the 
value of N, ata depth of 15 ft (4.58 ‘m) in Boring No. 2 is about 13, the 
laboratory test data would indicate a Ryeetaction resistance for this sand 


| 
| 
7 | 
wi penetration tes ata tO assess stress ratio Causing lquelaction | 


corresponding toa cyclic s 

it may be seen ‘that this ifi it to the sand near the 


a 
would be much less than the earthquake-induced stress ratios estimated to range 4 


a between 0.155 and 0.19. In fact, it would indicate a factor of ‘safety in the | : 
= of 0.71-0.87. Again, if the sand nearer the shoreline were of somewhat A 
= characteristics, the situation would be even less favorable than that bes 


indicated by the preceding figures, 


- On the basis of these evaluations, it should not be considered surprising 


that liquefaction near the shoreline was so extensive in effects. 


en 
_ The low unit “weight of the pumiceous “oy at the La Playa site would = 
a - appear to have been a significant factor in producing the liquefaction that occurred. 4 
‘The primary reason for | this is the very | low effective stresses produced a ; 
sand at significant depths below the water table. 
a Consider, e.g., the conditions that would have developed at La Playa if the © . 
- sand had been a typical quartz sand with the same resistance to liquefaction ile 
the actual sand containing pumice, i. e., ‘Tequiring an induced 
For the quartz sand, the unit weight of the moist sand might have as 
about 105 Ib/cu ft (168.1 g/cc), and the total unit weight about 115 joie 
ft (184.1 g/cc). For a water table at a depth of 5 ft (1.53 m), as existed in 
the highly liquefied zone, the at a of 15 ft ™ 
while the total pressure w would be be 
= 5x 105 + 10 x 115 = 1,675 psf. 


Thus, the stress ratio induced the be about 
0.65 


For of Tanging from 0. 15 g, this leads to values of cyclic 


Thus, depending on the acceleration there i is a very 


that a Fe a Lin of quartz sand would = have liquefied at all, in contrast — 


0.98 


EARTHQUAKE- INDUCED = 517 

to the much lower factors « of safety indicated for the pumiceous us sand (0.63-0.97),  §& 

the evident high degree of liquefaction of the area. 
_ It may be concluded that very lightweight saturated cohesionless soils are fj 

; Fi much more vulnerable to earthquake-induced liquefaction than typical quartz: : 


‘oats and should be treated with special caution in seismically active regions. 
NGINEERING DIVISION, 


a tes Sesion presents the results of a field and ial investigation of 
_ the extensive area of liquefaction which occurred at La Playa on the shore + 
of Lake Amatitlan in the Guatemala earthquake of 1976. The investigation leads ~ ~l 


the following « conclusions: ong ni 7 gairod bish py 


The soil in which liquefaction occurred was a layer of sand containing 
- particles of pumice which occurred between depths of about 5 ft-70 ft (1.5 ‘4 
m-22 m) or more below the ground surface. It was covered by a surficial layer _ 
of lightweight pumice sand and because of the pumice particles i in the liquefied 
layer its total unit weight had the relatively low value of about 90 Ib/ cu ft 
ae In spite of the fact that the sand is somewhat lighter in weight in sand ; 
deposits which have liquefied in other earthquakes, its liquefaction characteristics 
are apparently i influenced by ‘the same factors as other sand 
overall behavior i is consistent with that exhibited by other sa sands. ig 
_ 3. The penetration resistance of the sand at the boundary between liquefied 
- and nonliquefied zones is in good accord with previously developed empirical 
correlations between liquefaction potential and the standardized penetration — 
€ resistance (N ,) at which liquefaction can just be expected to occur. Sst bag” 
4. The behavior of the sand in the liquefaction and nonliquefied zones was 3 
consistent with ‘experimental- -analytical pt predictions of liquefaction potential based 
on the averaged results of cyclic loading triaxial compression tests performed 
on undisturbed samples to evaluate the liquefaction characteristics of the sand | 
and conventional procedures used in a sa with this type of test data 
toevaluate liquefaction potential, 
5. The high degree of liquefaction at the La Playa site was probably due 
in large measure to the lightweight nature of the pumiceous sands. A typical 
- quartz sand at the same site with the same cyclic load characteristics as the 
sand containing pumice might well have remained stable in spite of the earthquake 
_ shaking. Consequently, lightweight cohesionless soils should be treated with 
_ special caution with regard to their liquefaction potential in seismically active 
6. The results of the investigation provide a an extremely useful case histor i 
in which field data on soil characteristics in an earthquake-liquefied zone and — a 
a nonliquefied zone can be correlated with field performance, thereby supple- — 
_ menting the limited number of available case studies of this type which -_— 
i be used as a basis for predicting probable behavior at other sites. The results 4 
also t tend | to corroborate currently used procedures for evaluating liquefaction 
: potential although in the case of the laboratory test data, this clearly depends : 
on the degree to which the in-situ properties of the soil are represented by 
the “‘undisturbed samples extracted from the deposit. 
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In spite of the inevitable limitations of any ‘this it is 

oe that the results of the investigation, used with appropriate judgment, will 


: and thereby to an improved ee, capability of this type of behavior in es 
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Acoustic Emission MontTorinG OF 


"By Robert M. Koerner,’ M. ASCE, W. Martin McCabe, pans: wot! 


ima a ai we previously published papers (1,2,3), the acoustic emission (AE) 
method for monitoring soil stability was described. Briefly, the technique is 
based on monitoring nonequilibrium conditions within a soil mass by using a 
pickup transducer (i.e., , accelerometer, hydrophone, or high frequency — AE 
_ sensor), which is connected to an AE device. This device amplifies, filters, 
and then counts the number of threshold crossings that the incoming AE pulses 7 
; _ create. These counts are directly related to the instability of the material producing — 
t ~ them; e.g., high counts are indicative of major instability, low counts indicate _ 
_ minor instability, and zero counts indicate stability. The technique has been 
a 7 used for earth dam stability monitoring and is being investigated for a number 
7 of other applications, one of which is the monitoring of water flowing — 
earth masses, e.g., seepage in earth dams. val OF 
¥. _ This technical note is presented with the goal i in mind of showing the technical | 
feasibility of the method and of quantifying certain aspects of AE monitoring © 
of seepage of water through soil. The work is a laboratory study in which 
a both clear water and turbid water were passed through a column of soil at J 
various flow velocities. Flow ‘Tates a and AE rates from these laboratory tests 
are then compared to a previously-conducted field study at the site of reservoir 
leakage beneath a small earth dam. 


The test apparatus ¢ consisted of a plastic cylinder 20 cm ‘in 
x) cm high. Water flows into the base of the cylinder, then upward through 
a perforated base plate, and then through the soil column. After passing through 
the soil, the water is collected in a reservoir, and under steady-state seepage 
conditions, a flow rate is determined. Fig. 1(a) shows the experimental details 


Prof., Dept. of Civ. Engrg., Drexel Univ., Philadelphia, Pa. 19104. 
?Research Assoc., Dept. of Civ. Engrg., Drexel Univ., Philadelphia, Pa. 19104. 
Grad. Student, Dept. of Civ. Engrg., Drexel Univ., Philadelphia, Pa. 19104. 
- _ Note.—Discussion open until September 1, 1981. To extend the closing date one month, - 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on June 24, 1980. 
¥ This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings _ 
ad of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT4, April, 1981. 4 
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- soil per liter of water. The experimental setup for this c case is shown in Fig. 
- 1(6). The soil through which the water (both clear and turbid cases) ae 
was a poorly-graded, angular, medium-sized gravel. It was placed in a medium - 
Pa dense state with a resulting dry density of 1.5 g/cm’. This type of coarse- aed 
soil was selected due to its high permeability which allowed for a wide = 
of flow rates to be tested. AE monitoring was performed by embedding a 
waterproof hydrophone at the midpoint of the soil during its placement. ‘The 
_ hydrophone, which is most sensitive perpendicular to its longitudinal axis, wa 
placed as shown in the sketches of Fig. 1. The cable from the a i as 
was then brought out of the test cylinder to the AE monitoring device. ofl 
_ system used has a voltage sensitivity of 10.5 V/psi and was used with high-pass 
and low-pass filters set at 200 Hz and 50 kHz, respectively. The threshold 
youre was adjusted so as to screen out the entrance water valve noise at 


Wi 


the highest flow rate, and was used at this ie for oral tests. ‘Thus, all equipment ve 
settings remained constant throughout, "which hopefully allows for relative 
comparisons to be made between test results. 


Pes allowed for an AE rate to ; be plotted against time from the start a 
- the test. Fig. 2 presents this information for flow rates from 38 ml/s-120 ml/s. 
For this and the two. curves to follow, no data points are shown due” toa 
7 ‘relatively high degree of scatter; thus the curves represent average trends. The 
significant region of these curves is the terminal long-term value of the AE 
-* under constant flow. These values are seen to gradually increase with 
increasing flow rates. Additional tests were conducted up to flow rates of 200 
m/s where the trend continued (see Fig. 3). Also of significance is the minimum 
value of flow which i is detectable using this test setup and monitoring equipment. 
‘In Fig. 2 this value is seen to be 45 ml/s. While the lower flow rate of 38 
ml /s did give an initial AE response, it decayed to lower than threshold value — 


ho: the case of clear water parate 
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- known but could have come from an initial ‘movement of the embedded pickup 


: - under each i increasing flow rate. It should be noted that this AE pulse gradually _ 

: decreases with higher flow rates and is entirely absent at the maximum flow 
rate of 200 ml/s. This highest flow rate value of 200 ml/s corresponds to 

a Reynold’s 's Number of 0.2; thus all flow values are well within the laminar 


tests is seen with ‘the following exceptions: 


The AE rate at an equivalent flow rate is much higher fo for ony water Pe 
_ than with clear water. For example, at 70 ml/s flow rate, the clear water AE = 
rate is 1.5 counts /sec whereas the turbid water AE rate is 2.6 counts/sec. a 
2. The minimum detectable flow rate is much lower for turbid water than 


for clear water. In this case, the minimum detectable value is 10 ml/s, “anes a 


3. The AE pulse at the initial portion of the test is relatively larger in these 4 
turbid water tests than it was with the clear water tests. This seems to hypsenendll 
that a somewhat greater location disturbance to the Lome hydrophone occurs 
turbid water than we, clear water. 


the voir supplying the turbid water to test cylinder. 


Conneanion WITH ‘Seerace Test 
was used to seepage axis of an dam, 
which was leaking with what appeared to be clear water (4,5). Under consideration 
was a 300-m section along the dam axis. in which the owner located ‘lh 
- gandomly- spaced, open-ended borings. He conducted flow rate tests in these _ 
_ borings with measured values from 0 ml/s-250 ml/s. The higher values — 
concentrated in a specific 62-m long location. Having access to the location - a 
- of seepage from within the cased borings, the AE activity was monitored at F 
the same ; locations using a | heavy steel | wire | as a wave ‘€ guide to to transmit the 


 a_iee systems were somewhat different between this field test and the 
_ laboratory tests described herein, the flow rates and resulting AE rates do compare 
2 within reason. (For the maximum clear water laboratory test flow rate of 200 — 

_mi/s, the AE rate was 29 counts /s sec. For direct comparison t between laboratory 

_ and field tests, it should be recognized that some attenuation probably occurred 

im the field as the signal traveled up the wave guide, bringing the two AE im 
count rates to a reasonable comparison.) Perhaps of the greatest significance _ 
was the fact that the field =, of maximum seepage location was eer 


: The parallel set of test results using turbid water, instead of clear water. 1 - 
4 
) 
CMUSSIONS the DOLUOM OF the DOTINE tO Ihe PICKUD | saucer. 
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s In examining the AE method for earth dam 1 stability monitoring, » the question 
While a modest field test had been previously completed, a more - controlled ' 

_ feasibility study, along with base-line data have been presented in this paper. _ 
Clear water seepage flow rates of approx 45 ml/s are required for AE detection 
using equipment and sensitivities described herein. For the case of turbid water 
seepage (as with an earth dam “undergoing some degree of piping), the minimum — 
detectable flow rate drops to approx 10 ml/s. The behavior of AE rates for 
flow velocities above these minimum values appears to behave as an exponentially 
function (see Fig. 5). it can be that high flow rates 
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A.E. RATE (COUNTS/ sec. ) 


Turbid Water Seepage Conditions 
proportionately many more emissions than do | the lower flow rates. A modified s 
7 experimental setup to test still higher flow rates would be necessary to evaluate te 
AE AE generation the flow into tl the turbulent 
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4 Faure in Benerrtr/Cost ANALYSES” 


Discussion m by Francis G. McLean,” ASCE 


_ The authors have proposed a thought provoking course of action. Their 
4 sasoakaion addresses both the economics and safety of a proposed project 7 
_ without, in the opinion of the writer, a proper separation of the two topics. 
_ A separation of the topics, in light of the actual planning, engineering / design, 
construction and d operation / maintenance sequence whi which r results i ina completed, 


, othe true safety of a project results from the quality of the engineering /design a 


_ (13,14,15) and construction sequence and is not impacted by the economic 


‘evaluation. Continued safety, of a Project depends on operation and maintenance 
procedures, ‘including performance - surveillance and continued technical assess-_ 
ments of criteria. A project could, for example, have the cosis of a hypothetical 
piobability of failure included in the benefit /cost analysis (BCR), show benefits — 
adequate to warrant construction, and still be grossly unsafe if designed, — 
constructed, operated, or maintained inadequately. In fact, if the costs associated 
with safe design, construction operation are unconservatively (under) 

estimated during planning, the BCR would | ‘appear more attractive. The real 
_ costs directly associated with achieving a safe project are currently incorporated | 


into the economic ovemation process in each of the areas delineated in the ¥ 


Site Investigations. —A thorough understanding of the site conditions, both — 
hydrologic and geologic, are necessary to achieving a safe project. Experienced 
personnel and consultants, coupled with adequately planned and executed site 
_ investigations and evaluations, are requisite to attain this understanding. The 
more difficult (less favorable) the site conditions, the larger the costs included 
_ in BCR calculations during the planning stage to accommodate a more extensive 
j program of investigation and 1 testing. In addition, | ‘more effort may be expended 


achieve more realistic BCR results. 
Design. —Builders of major dams are generally aware sof the hazard potential 


zones; or, potential for liquefication may be determined to exist, resulting in 


more plastic. cores provisions for s specially processed filter and = 
1980, Gregory Ba Baecher, and ‘Richard De Neufville 
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poor foundation materials or insitu densities being improved. 
hh each instance project costs may be increased and probability of failure reduced. 
_ Engineering of these projects is continued through the construction stage in 
order to safely adapt the design to conditions actually uncovered at the site. a 
4 Construction. —During this phase the precise lines of the designer’s plans 
turned into. Teality through the efforts of the contractor and owner’s 
_ representative, who supervises and inspects construction to assure safe adaption ~ 
_ of the design to the site. The experience and integrity of the contractor oT a 
_ inspection personnel account for the actual safety of the project, providing — 
that original design concepts were adequately formulated for the site. Attention 
to details, such as the techniques of foundation treatment and compaction of 
the core zone at the rock contact, hold the key to construction of a safe dam, 
and rigorous inspection to assure that proper materials and methods are employed 
is necessary. These costs are reflected in the estimates. The safe adaption = 
_ a design to actual conditions often results in cost overruns, to the betterment 
of safety and detriment of the BCR. 
_ Operations and Maintenance. —The profession is cognizant of the record of 
- failures during the early phase of a project life. Therefore, close monitoring 
of first fillings and project performance for the first 5 yr-10 yr is the usual 
7 practice on projects. Instrumentation is provided to aid in this monitoring. Regular 
inspections and evaluations of instrumentation data are required, and timely ¥ 
maintenance of the project assure a continued understanding of 


_ The BCR is the yardstick which is used to measure project * “‘*profitability”” 

to determine whether a Project should be built, or to select between alternatives. 

‘During planning site specific project feasibility and associated costs are established 

so that realistic BCRs can be estimated for purposes of decision making. This — 

BCR reflects the cost of project safety as indicated previously, and is nal 

Beste on ‘the foregoing it is not apparent that the inclusion of a discounted 

cost of failure in BCR calculations is ‘‘. . . more appropriate for the protection 

of the public than an implicit zero . probebility of failure’ as stated by the 
_ authors. The only apparent protection offered to the public would bethe significant 
- reduction of the BCR calculated for all projects, whether marginal or much 

needed, without an actual site specific assessment of safety. To the writer's” Ss 
_ knowledge n no recent dam failure, or near failure, (Baldwin Hills, Walter Bouldin, co , 
Teton, Lower San Fernando, Kelly Barnes) would have been prevented, or the 4 
_ public better protected, by the author’s proposed course of action. In addition, 
_ if this approach is to be taken, a uniformity of application would i require that 
_ the cost of failure due to > earthquakes, hurricanes, tornadoes, etc., , be included 
for all major construction n projects which pose a public hazard, i.e., buildings, cs 
liquified natural gas (LNG) facilities, refineries, etc. 


a The writer has avoided discussion of the data presented by the authors in 


an attempt to provide separation and perspective of the topics of safety and 
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economics. Probability methods wot be applied to the technology of 
dam design and construction to improve safety; or, a site specific design could 
. completed and assessed for the probability of failure to see if the resulting — 
3 risk is acceptable to the public interest. However, beneficial or extensive 
i in these areas appear to be somewhere in the future. Temitdics _ 
: _ ‘The discussion offered herein is based on the experience and p professional ; 
associations of the writer and is not offered as, nor constitute, the position 
of the Corps of Engineers. A debt of gratitude is due to the many professionals — 
who have engaged in the design and construction of safe dams through the 
years, the fruits of whose technical and beneficial legacy we all Geer ot 
13. Embankment—Dam Engineering, Puaitieale Volume, R. Hirschfeld and S. Poulos, 
ed., John Wiley and Sons, New York, N.Y.,1973. 
«14, Peck, R. B., “‘Let’s Get It Straight About Those Dams,” The Military Engineer, ue 


. ‘Federal Guidelines for Dam Safety,”’ Federal Coordinating C Council for Science, 


n. M. 


_ The studies summarized in this brief note represent a substantial step forward. 
a It would be difficult to disagree that BCR should include recognition of the 
_ probability of failure, or t that the ‘Suggested default value of 10" “/dam yr i is 

“ represents 

_ The authors indicate that reduction of this default value, or elimination “a : 

_ the need for relying on such a value, will depend on the success of research 

leading to analytically derived probabilities of failure. Although not intended 

_ by the authors, this might be construed to: suggest that reduction of the default : 


value by probabilistic analysis ‘could somehow | increase the safety “of dams. 
Such a view would be untenable, as the safety of a dam is determined by 
= _ The failure of Teton Dam dispelled any complacency regarding the inherent 
safety of major dams designed by experienced organizations. This failure would a 
presumably be assigned by the authors to “unexpected “causes,” but both 
investigative panels concluded that the failure should not have been unexpected 
- a and that the design was deficient in a number of respects. These deficiencies if y 
could have been lumped together as a failure to provide defenses in depth a 
b= possible weaknesses in the first line of defense (the impervious core: 
existed; it was not used. of 
/_Whether we like it or not, there are aspects a 


= Prof. of Foundation Engrg., , Emeritus, Univ. of Illinois at Urbana- -Champaign, 
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AP RIL 
dam design, that are not now may never be to quentiative 
analytical procedures. Most of these (except for overtopping, seldom a problem > a 
= in recently designed dams because of present conservatism in spillway require- 


7 " ments) appear in the foundation or at the interface between foundation and 
| iterate Most, like those at Teton, would yield to the provision of nile 


lines of defense against qualitatively assessed Most, in other) words, 


we discount judgment and the potential its in 
defenses in depth, we may discourage the use of judgment and may never 
decrease the default value below 10~*/dam yr. If we encourage it and 
apply it in developing defenses in depth, the historical record within a few 
: decades may justify a a lower default vai value. This may be a more promising approach — 
than awaiting the development of analytically derived probabilities of failure. 
At best, these analytically derived probabilities could only improve the reliability 
of estimating the risk of failure; iad a powerless to alter that risk» 
= creating safer dams. wan) 


Closure Gregory Baecher,° A. M. Marie Pate,’ 

wad Richard de Neufville,* M. ASCE 


_ That the failure of Teton Dam, as Peck remarks, should dispel complacency — =) 
_ regarding the inherent safety of major dams is a point to be emphasized. The 
fact that nearly all major dams prove safe, a truly remarkable achievement, 
- must not obscure the complementary fact that some dams fail catastrophically. € 
It should be stressed in this connection that analytical | estimates of the 
" probabilities of these failures are difficult or impossible to obtain. As a practical 
matter, civil engineers go to great lengths to guard against every potential cause 
@ failure they can identify. The risks remaining thus principally come from “ 
- sources not analyzed. These may be due to poor physical understanding, simple : 
- omission, or gross error. For purposes of analysis | of the planning decisions, s 


the result is the same, however: the identification of these risks is beyond 


We thus entirely concur with Peck’s suggestions that it may be best to justify 
J Beis. the default value for the probability of failure upon the = favorable 


planning dam projects, it is our responsibility to ‘into the 
- eee the real possibility that any dam might fail. We do not know in — 
_advance which dams will fail and how these failures will occur. If we did, 
we re would defend against these eventualities or not build the project. ace arly 


* Asst. . Prof., (Civ. -Engrg., Massachusetts It Cambridge, Mass. 
Asst. Prof. - Dept. a Civ. Engrg., Massachusetts Inst. ‘Tech., Cambridge, Mass. 
"Prof. _ of Civ. Engrg. and Chmn., Tech. and Policy Program, Massachusetts Inst. . of be 
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projects, those posing larger hazards to life uae property must have their benefits — 


_ As Peck and Mclean remark, analysis of probabilities will not in itself increas 
_ the safety of dams. This safety ultimately depends on the ex experienced judgement — 
of the designers and the care taken to ‘‘provide defense in depth.’’ The inclusion 
of failure risks in BCR is not immediately aimed at decreasing failure rates, 


but at better allocating limited resources among public investments. Yet in the 


_ longer run, the resulting tendency to avoid dam projects posing higher risks = : 


to the public might enhance safety. om) bates! 
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(OF F EARTH Dam! MATERIALS 


by Marvin R. Pyles,’ A 

‘The ‘authors have taken a long step in a direction. Erosion 

Ph of soil particles through a continuous medium such as a filter zone in an earth 

_ dam has been studied extensively, but erosion of soil particles through a crack 

- 7 has not. Producing a crack by tensile deformation while applying a hydraulic 


gradient should be a . good technique for ‘studying this behavior, but — 


_ The apparent erodibility of any soil is a function of: (1) The characteristics 
of the soil; and (2) the environment in which the erosion takes place. Of particular 
note in the authors’ test apparatus is the size of the opening (8) between the 
: top and bottom halves of the cylinder. T. In the test series on both soil A (moraine) 
and soil B (silty- clayey fine sand), the maximum size opening reached, as indicated 
7 by Figs. 7, 8, 9, and 10, ranged between 2.7 mm and 5.8 mm. Since the method 
of measurement employed by the authors identifies an eroded soil particle o 
‘one which passes through the opening between the top and bottom halves of 
_ the cylinder, it is important to consider the size of that opening. gatbn. 3 oadw 
probability of one of many soil particles passing through an 
of given dimension would look something like the relationship shown in Fig. 
- 14 as a function of the ratio of opening size to particle size. The reason that 7 
- particles smaller than the opening do not have a probability equal to one © 
of passing through the opening is that there are many particles in the population, 7 


and they can interact, , bridging across the opening and preventing ‘€rosion. In 4 


_ the authors’ test results, the ratio of maximum opening size to maximum particle 

a ‘size for soil A ranges from 0.88 to 1.46. For soil B, the ratio ranges from 
1B. 3 to 39.3, an order of magnitude greater than that of soil A. It is easy 
to see from the schematic probability plot in Fig. 14 that soil B is much more 


na February, 1980, by Erik I. Hjeldnes, and Bhagwat V. K. Lavania (Proc. Paper 15220). ; 
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likely to erode through the opening between the top and bottom cylinders than 
= is soil A, giving soil A the appearance of being “self-healing.” *»* In the authors’ 
7 tests, it is not a property of the soil, but a constraint imposed by the test 
is _ apparatus that has been illustrated. To compare the soils properly, they should — 
be tested in a manner s such that ‘both a are equally affected by t the « constraints — # 
| _ The authors | do not make it clear in the description of the testing technique 


_ whether the samples were saturated after compaction and before testing, or — 
were tested at the compaction water content. If the samples were tested at _ 
S compaction water content, then capillary stresses could have played an 
important role in the tensile behavior of the soil, and thus the results of the 7 
erosion tests. The condition in the field that is of concern is cracking in the - 


e of an “embankment, and subsequent leakage. Depending on where and _ 
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14. 4.—Schematic Erosion | Probability Relationship 


Boe both conditions may be of interest, and show significantly different behavior. 
| ‘The tensile stress field — by the apparatus used by the authors is surely 


x _ when cracking occurs in an ‘embankment, the core could be saturated or not, 


deformation at the center of the sample. This will mt any a tensile stress 
_ distribution like that shown in Fig. 15 obtained from a linear elastic finite element 
analysis of the test apparatus. The magnitude of the stress in the soil immediately 
yy to the division between the top and bottom halves of the cylinder 
Pe is on the order of 10 times the mean tensile stress on the cross section. Since 
a soil i is not ‘elastic, the actual stress concentration will no doubt be something 
less than this, but still significant. The presence of stress concentrations does 
aie tae invalidate test peneits, for comparing one soil to another, but 
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tensile stress concentrations ns could have initiated csacking a and subsequent leakage 
at a mean stress deformation level that would.s not have produced 


a dam. This is ‘apparently beyond the of the authors’ stedy, 


be ignored i in the long run. on} 
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FINITE ELEMENT MESH STRESS DISTRIBUTION 


AG. 15.—Stress Analysis Results for « Approximate Maximum Mean Tensile Stress | 


The technique presented | by the authors may have promise for the comparative - 
indexing of ‘soils for use | in earth and rockfill dams, but much study needs 


be done to answer the q questions Presented here. 


Closure by Erik I. Hjeldnes* and Bhagwat V.K. Lavanias 


The writers would like to thank Pyles for his in hisi interest in our paper. His comments — 
- us an opportunity to explain more thoroughly the testing technique and = 
the results achieved. The following sections address specific points raised by _ 

Pyles. To some extent the closure also bears evidence of a valuable personal — ¥ 


“communication with G. A. Leonards. 
_ Laboratory Approach.— Apparatus and testing technique are very simple and 


pom of course many limitations. However, they served the main purpose of 
_ “Assoc. Prof., Univ. of Trondheim, The Norwegian Inst. of Tech., N-7034 Trondheim- 


*Executive Engrg., Tehri Roorkee, State of U.P. 


to the physical considerations of the test 
. ve procedure, there is a need to consider soil and water chemistry, and their effects 
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‘ao teste, iis e., to study the behavior of compacted d soils with 1 respect to leakage 
- and erosion under simultaneously acting tensile strain and hydraulic gradient. — 


In our opinion this is a very important requirement of the testing conditions. 


- Erosion Properties of Soils A and B.—The writers are | of the opinion that 
the opening (8 in Fig. 2) through which the leakage has to pass does not seriously a 


affect our conclusions concerning the erosion behavior of soils A and B. As 4 


the leakage, the two soils responded quite differently. 4 
_ In soil A the erosion started at a relatively small given deformation. As ite 
deformation was kept constant, the erosion increased, reached a maximum, 
and thereafter decreased to zero 4(b)]. This phenomenon was designated 
the v writers as as “self-healing,” 


= 


started, it continuously ‘caused blow off w without further 
- Consequently the deformations at start of erosion for soils A and B as on 
were observed in the tests [see also Fig. 11(b)] are not in accordance with 
“the discusser’ s Fig. 14. ‘From his diagram it is more likely that erosion starts ; 
a ‘first in soil B. This discrepancy shows the importance of taking into consideration _ 
the soils’ different grain size distribution and behavior under tensile straining. — x 
_ For soil A it is very interesting to compare, at a given deformation (8 < 
8 ), the observed erosion WY, ) with the volume (V,,) of soil particles smaller 
than 8 in the volume 48. Here A 8 is approximately equal to the soil volume — 
affected by the hydraulic gradient, and therefore V,, should be the upper 
- for the erosion. Simple calculations show that Vi<< V,,, a result that confirms 
. s “self. -healing’’ effect and shows that the outlet size is not the only restriction. 
gy _ Saturation of Samples.—The curves in Figs. 7-10 show that leakage (q,,) 


started for deformations (8) much smaller than deformation at blow off. It 

is therefore evident, because of the water flow through the sample, that the 

soil close to the fracture plane was highly saturated, and that the capillary — 
forces had no significant influence on the observed behavior. 

_ Stresses on Fracture Plane.— ae is obvious that the ‘stresses o on the fracture 
plane are not uniform, and that there probably exist stress concentrations. The — 
_ discusser has shown this (Fig. 15) for an elastic material in tension. It is, however, — 
appropriate to mention that for all tests with hydraulic gradient (Figs. 7-10) 
the total force (P) on the fracture plane is a compressive force, except for 

a possible small tension at the very start in some few of the tests. Hence _ ; 
average total stress P/A is a compressive stress. ‘(Tensile deformation 
does not necessarily imply tensile stress.) 
: The total force can be divided into a resultant pore pressure (U) and a resultant _ 
effective stress (P’), i.e., P = . + U. The results observed indicate that 

in these tests U is predominant and consequently P’ is small. In further and = 

more advanced tests of this kind one should try to measure the ‘distribution ny 
~ magnitude of stresses on the fracture plane, at least the pore pressures. = 
_ Chemical Effects.—As mentioned by the discusser, soil and water chemistry — 
could have some influence on the erosion behavior. Primarily one would expect hod 
such effects in fine- “grained soils. The writers are of the opinion that in most 


a erosion behavior of soils. Emphasis should preferentially be placed on 
large-s scale experiments, as “prototype- -like’’ as possible. The scope must be 
to improve the design criteria and the guidelines for choosing the most suitable — ’ 


materials for earth and rockfill dams. 


=  _ cases this will be of minor importance, but research aiming to clarify the question j _ 
s __ The writers hope these remarks have clarified the points raised in the aint. ’ 
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